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TO THE COMMISSIONER FOR PATENTS: 

I, Dr. Clements C. Lambeth, declare as follows: 

1. I am a co-inventor of the subject matter in the above-identified patent application 
and I am employed by Weyerhaeuser Company as a Tree Improvement Department Manager 
and Senior Scientific Advisor. 

2. A copy of my curriculum vitae is appended hereto as Attachment B. 

3. I have considered the Office Action dated October 18, 2004 and the decision by 
the Board of Patent Appeals and Interferences, dated June 6, 2006, issued in the above-identified 
application. It is my understanding that the Board has reversed the Examiner's rejections of 
claims in the application based on lack of written description and lack of enablement, and has 
affirmed the rejection of claims based on obviousness in view of Bridgwater (1992) in 
"Handbook of Quantitative Forest Genetics," Kluwer Academic Pub., Dordrect, The 
Netherlands, pp. 69-95, in view of El-Kassaby and Ritland (1992j Theor. Appl Genet, 83(6- 
7):752-8 and Stoehr et al. (1998) Can. J. For. Res. 25:187-95. 
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4. It is well known by those with knowledge in the field of tree breeding that it takes 
a very long time to validate a breeding method with tree species due to their long life span and 
the long time it takes to induce flowers. Therefore, in order to validate some of the benefits of 
the tree breeding method claimed in our pending patent application serial number 09/618,307, 
referred to as polymix breeding with parental analysis (PMX/WPA), my colleagues and I, with 
the aid of a contract laboratory Select Breed Services, designed and conducted a quantitative 
genetics simulation study in 2002 comparing PMX/WPA to a conventional full-sib breeding 
system. The simulation study, described in the research report attached hereto as Attachment C, 
assumed certain population statistics, heritability, inbreeding effects and selection through eleven 
generations. The simulation study demonstrated that PMX/WPA resulted in greater genetic gain 
over ten generations of recurrent selection under both restricted and unrestricted inbreeding 
scenarios than full-sib breedings, at a significantly reduced cost. In addition, the parental 
breeding values were more accurately estimated under the PMX/WPA system, implying that 
orchard roguing and family allocations based on progeny performance would be more accurate 
than with the conventional full-sib systems. 

5. My colleagues and I conducted a study in 2005 in order to test the contribution of 
males in a polymix breeding system with parental analysis (PMX/WPA). In this study, ten 
females were pollinated with a 35-male pollen mix and ~ 44 progeny per female were subjected 
to paternity analysis. Paternity assignment in this study was based on the use of 12 chloroplast 
and 3 nuclear microsatellite markers, which allowed for the unambiguous identification of 70% 
of the fathers and approximately 80% of the possible 350 crosses had at least a partial paternity 
assignment. The results of this study are presented in FIGURE 1, (appended hereto) and 
demonstrated that males contributed to fertilization in a near random fashion and the frequency 
distribution of the contribution of males to the progeny was excellent. Because of this random 
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distribution of paternal success, many effective crosses were created, with nearly 80% of the 
possible 350 crosses having at least a partial paternity assignment. In contrast, in a circular, full- 
sib mating system, a preferred method in the field, only 70 crosses would be created. The 
random distribution of paternal success demonstrated in this study resulted in more crosses than 
would be possible in a traditional breeding system, which would therefore produce more 
diversity than traditional breeding systems. The increase in diversity that is made possible with 
the PMX/WPA breeding method allows for more choices of unrelated or weakly related trees for 
inclusion in a breeding group for use in a next generation of tree breeding to allow for increased 
genetic gain through selection and better inbreeding control. 

The results of this study also demonstrated that crosses between closely related pollens 
resulted in a low frequency, thereby adding an additional advantage to the PMX/WPA method. 
In this study, some of the pollen in the polymix was related to the female (self, full-sib, parent). 
Based on a mathematical calculation assuming equal contribution of all males, the expectation in 
this study was 1.25 progeny per each of the 350 potential crosses. However, what was actually 
observed was that unrelated crosses averaged 1.36 progeny per cross while self and close relative 
(full-sibs, parent-offspring) crosses had only 0.37 and 0.61 progeny per cross. Therefore the 
results in this study showed that closely related pollens were at a fertilization disadvantage 
and/or their progeny were selected against in the early stages of life. It is well known by those 
with knowledge in tree breeding that mating of close relatives is not desirable in loblolly pine 
due to inbreeding depression in this highly out-crossed species. Therefore, the low frequency of 
inbred progeny observed in this study demonstrates an additional advantage of practicing the 
claimed method which further increases the potential for genetic gain that can be achieved using 
the claimed PMX/WPA tree breeding method. 
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6. The PMX/WPA tree breeding method satisfies a long-felt need in the field of tree 
breeding. Prior to the publication of our journal article describing the claimed PMX/WPA 
breeding method (Lambeth, C.C. et aL, Theor. Appl. Genet. 705:930-943 (2001)), there was a 
long standing need for a workable solution to the use of a polymix breeding scheme in isolation 
from other crossing schemes in tree breeding programs. The potential value of polymix breeding 
method has been long recognized by authors in the field of tree breeding. (See for example, 
Shelbourne C J.A., "Tree Breeding Methods," Tech. Pap. No. 55, New Zealand Forest Service, 
(1969), attached hereto as Attachment D; and Burdon R.D. and Shelbourne C.J.A., "Breeding 
populations for Recurrent Selection: Conflicts and Possible Solutions," NZJ. For. Sci. 1(2): 174- 
193 (1971), attached hereto as Attachment E). However, the limitations of polymix breeding 
have also long been recognized and include lack of full pedigree control, which limits genetic 
gain due to lack of breeding value on the male parent, and inability to control inbreeding levels. 
These limitations have led to the predominant view in the field of tree breeding (prior to the 
publication of the methods of the instant invention) that polymix breeding should only be used 
for estimation of breeding value of the parent and not as a basis for selection of the next 
generation of breeding (see for example, McKeand, S.E. and Bridgwater, F.E., "Third 
Generation Breeding Strategy for the North Carolina State University-Industry Cooperative Tree 
Improvement Program," In Proc. IUFRO Resolving Tropical Forest Resource Concerns Through 
Tree Improvement, Gene Conservation and Domestication of New Species, Cali, Columbia 
pp. 234-240 (1992), attached hereto as Attachment F; Burdon, R.D. and van Buijtenen, J.P., 
"Expected Efficiencies of Mating Designs for Reselection of Parents," Can. J. For. Res. 
20:1664-1671 (1990), attached hereto as Attachment G). 

For example, as described in McKeand and Bridgwater, (1992) (see Attachment F), 
university based cooperative tree breeding programs in the United States have used a dual 
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complementary breeding system in which polymix crossing is used to estimate breeding values 
and full-sib crossing is done to produce a population for selecting the next generation of trees. 
This dual complementary breeding system is known to be costly and time consuming. However, 
prior to the publication in 2001 (Lambeth et al, 2001) of the method of the claimed invention, 
there was no proposed solution to combine the polymix breeding system with pedigree analysis 
to produce the method as claimed in the present invention. 

7. The PMX/WPA breeding method as claimed in the pending claims 20-31 of 
application serial number 09/618,307 solves the limitations of the polymix breeding system and 
eliminates the need for the dual complementary breeding system described above in which 
polymix crossing is used to estimate breeding values and full-sib crossing is done to produce a 
population for selecting the next generation of trees. The PMX/WPA breeding method is also 
less expensive and has the potential to provide more genetic gain in a shorter time span than the 
dual complementary breeding system. 

In response to our publication in 2001 (Lambeth et al., 2001) describing the method of 
the claimed invention, there has been acceptance of our breeding method by those in the field of 
tree breeding as evidenced by favorable comments we have received from scientists in the field, 
requests for collaborations, and copying of our breeding methods by others, as described in more 
detail below. 

The reaction my colleagues and I have received from several scientists in the field is that 
the claimed invention is a very creative solution to the problem of tree breeding and one in which 
many of them are interested in due to one or more advantages of simplicity, speed, excellent 
estimates of both parental and individual progeny breeding values, greater opportunity to reduce 
inbreeding and greater genetic gain than conventional full-sib systems. The following 
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representative examples are provided as evidence of acceptance by scientists in the field of tree 
breeding, collaboration and/or copying of our PMX/WPA breeding method by others: 

• Dr Dudley Huber, Quantitative Geneticist and Professor at the University of 
Florida has described our method as a novel and effective approach to an age-old 
problem (personal communications). 

• Dr. Satish Kumar and his New Zealand and Australian colleagues at ENSIS are 
very interested in our method and mention our publication as a "leading paper" in 
an email communication dated July 20, 2006, attached hereto as Attachment H. 

• Dr Terry Stranger, geneticist with SAPPI Timber Company in South Africa and 
his colleagues from the University of Pretoria have discussed the PMX/WPA 
breeding method with me and have launched a research effort to test the efficacy 
of our breeding method in commercial eucalypt species, as described in a poster 
article by Minique H. et Castro et al., 2004, attached hereto as Attachment I. 

• Boise Cascade Corporation has contracted with co-inventor Nicholas Wheeler to 
conduct a study on the utility of our PMX/WPA method in popular breeding and a 
paper was published on the results, Wheeler N. et al., Tree Genomics and 
Genomes 2:53-60 (2006), attached hereto as Attachment J. 

• Dr. Dario Grattapaglia has cited our paper describing the PMX/WPA method 
(Lambeth et al. 2001), and has employed a similar strategy for the improvement 
of eucalypt seed orchards in Brazil (see Grattapaglia et al., Theor. Appl 
Genet. 7 OP: 192- 199 (2004), attached hereto as Attachment K. 

• Dr. Rowland Burdon, Senior Research Scientist at the Forestry Research Institute 
in New Zealand, has cited our PMX/WPA breeding method and mentions its 
value in the control of inbreeding in Burdon et al., IUFRO Division 2 
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Proceedings, "Forest Genetics and Tree Breeding in the Age of Genomics: 
Progress and Future," Charleston S.C November 1-5, 2004, pages 75-85, attached 
hereto as Attachment L. 



8. AJ1 statements made herein and of my own knowledge are true, and all statements 
made on information and belief are believed to be true; and further, these statements were made 
with the knowledge that willful, false statements and the like so made are punishable by fine or 
imprisonment, or both, under Section 1001 of Title 18 of the United States Code, and that such 
willful, false statements may jeopardize the validity of the above-identified application or any 
patent issued thereon. 



Respectfully submitted, 




Clements C. Lambeth, Ph.D. 
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Weyerhaeuser Company 
810 Whittington Avenue 
Hot Springs, AR 71902 



Phone: 501-624-8510 

Fax: 501-624-8419 

Home: 501-627-1941 

Email: clem.lambeth@weyerhaeuser.com 



dements C Lambeth, Senior Scientific Advisor 



CAREER OBJECTIVE 



To work in a goal-driven, forestry research organization that places strong emphasis 
on technology implementation, environmental stewardship and employee job 
satisfaction. 



Degrees 
Earned 



EDUCATION AND TRAINING 

=> 1972 - B.S. in Forestry, Stephen F. Austin University, Nacogdoches, 
TX (Outstanding Senior Award, Texas Forest Service Scholarship 
Recipient) 

1 974 - M.S. in Forest Genetics, Yale University, New Haven, CT 
(Research Fellow Recipient) 
^ 1979 - Ph.D. in Forest Genetics, N.C. State University, Raleigh, NC 
(Weyerhaeuser Graduate Studies Program Recipient) 



Other Training 



Management Skills 

Podium Power (Public Speaking) 

Total Quality Modules (Weyerhaeuser) 

Managing for Continuous Improvement (Tennessee Associates) 
Covey's 7 Habits of Successful People 
People Skills 

Advanced SAS (Statistical Analysis System) 
Biotechnology in Forestry Short Course - University of Florida 
Various Microsoft Products Including Microsoft Project, Excel, Word, 
PowerPoint, Access 

Spanish proficiency through intensive courses and on the job usage 

Managing Multiple Projects, Objectives and Deadlines 

Diversity Training 

Self-Directed Work Teams 

Safety Including Behavior, CPR, First-Aid, etc. 



Other Interests Travel > Canoeing, Reading, Coffee Roasting, Building 



EXPERIENCE HIGHLIGHTS 



Weyerhaeuser Hot Springs, Arkansas 
1 992-Present Tree Improvement Department Manager & 

Senior Scientific Advisor 

Roles and Responsibilities (Managed up to 14 scientists and 
technicians): 

=> Lead a team that develops genetically improved stock for 4.5 million 

acres in the southeastern U.S. 
=> Provide periodic presentations and reviews to Senior Timberlands 

leadership for decision-making regarding the development, production 

and deployment of genetically improved stock in the southeastern U.S. 
=> Lead the somatic embryogenesis clonal development and testing 

program 

=> Provide the needed information to Southern Timberlands area foresters 
for successful deployment of genetically improved stock 

=> Represent Weyerhaeuser on executive committees of 4 University- 
based tree improvement and research cooperatives 

=> Provide consultation to Uruguay and Saskatchewan timberlands 
leaders on tree improvement strategies and help develop long-term 
plans 

^ Conduct research to improve the efficiency and genetic gain in 
breeding, testing, selection and orchard processes 
Serve on a company-wide Clonal Forestry Implementation team that 
develops the strategy for implementation of cloning on all 
Weyerhaeuser land 

=> Serve on a Seed and Orchard Functional Team that studies 
alternatives and makes recommendations on changes in strategy and 
operational practices that improve growth, cost and/or quality of 
Weyerhaeuser's Southern Timberlands through tree improvement and 
seed orchards 



Key Accomplishments and Communications: 

=> Presented a Southern Tree Improvement Strategy that was approved 

by the Senior VP of Timberlands 
=> Wrote the Southern Tree Improvement Plan which details strategy as 

well as breeding, testing and selection activities and their timelines for 

Weyerhaeuser's 4 million acres in 9 states 
=> Led in the integration of Proctor and Gamble, MacMillan Bloedel, 

Cavenham and Willamette tree improvement programs 



=> Played a key role in the development of a widely used Landform Table 
which details the cost and benefit of various silvicultural and tree 
improvement options for all southeastern landforms and allows the user 
to test the economic benefits of single or multiple early rotation, stand 
management options 

=* Developed an electronic south-wide seed catalog that is used by 
Operations and Research personnel as a planning tool for family 
deployment, new orchards and breeding work 
Played a key role in development of a long-term orchard production 
strategy that was approved by the Southern Timberiands VP 

=> Led a team review with the Southern Timberiands VP and the lead 
team which prompted Southern Timberiands to make a potentially $50 
million decision to plant North Carolina stock in Ark/Okla & NLouis. 

=> Developed a rooted cuttings research and implementation strategy that 
was approved by the Timberiands Senior VP and Southern 
Timberiands VP 

=> Managed a team that developed improved outdoor rooting techniques 

and mechanization of steps in rooted cuttings 
=> Developed a new breeding concept (PMX/WPA) to increase genetic 

gain and efficiency in improvement programs using paternity analysis 

Led in the implementation of top grafting, a new technology for 
breeding, new orchards and crown replacement in operational orchards 

=> Led in the implementation of control mass pollination (CMP) in 
operational orchards including a strategy presentation that was 
approved by the Southern Timberiands VP 

=> Led a multi-department (forestry, analytical sciences, solidwood, pulp 
and paper) team effort in a study of the inheritance and gain potential of 
wood quality traits in loblolly pine that led to incorporating them into the 
tree improvement program 

=> Developed tree improvement strategies for Saskatchewan Timberiands 
and COLONVADE in Uruguay and made presentations to the forest 
management leaders in those areas 

=> Worked with a team of outside scientists to redefine seed source 
transfer guidelines for southern pines 

1989-1992 

Smurfit Group Jacksonville, Florida 
International Forestry R&D Manager 

Roles and Responsibilities (Functional responsibility reporting to 
Timberiands VP's in Colombia, Venezuela, Mexico and the U.S.): 

=> Coordinated a multi-disciplinary and multi-country research team aimed 
at improving operations and yield in Latin America and the U.S. 

=> Advised the presidents and vice-presidents of Colombia, Venezuela 



and Mexican Timberlands organizations on R&D direction and led 

projects needing international coordination 
=> Arranged for outside consultants as needed when internal skills were 

lacking in developing forestry programs 
=> Conducted tree improvement research for Container Corporation in 

Florida, Alabama and Georgia 

Key Accomplishments and Communications: 

=> Developed breeding strategies for Colombia and Venezuela for several 

hardwood and pine species 
=> Conducted on-the-ground initiation of the Gmelina and Eucalyptus 

clonal programs in Venezuela 
=> Conducted research in controlled mass pollination in seed orchards 
=> Initiated the testing of exotic pine species in Florida 
=> Coordinated the publication of the Smurfit Group Forestry Research 

Report series 

=> Organized an international conference of the IUFRO Tropical Tree 
Breeding and Testing group in Cartagena and Cali, Colombia 



1985-1989 

Smurfit Carton de Colombia/Container Corp. Cali, Colombia 
Director of Forestry R&D 

Roles and Responsibilities (Managed up to 10 scientists and 
technicians): 

=> Advised company leadership in Colombia and Venezuela on research 

and reforestation activities 
=> Led a multi-disciplinary team researching forest nutrition, competition 

control, hybrid and species testing, vegetative propagation, and tree 

improvement of pine and Eucalyptus species 

Benchmarked tropical pine and eucalypt research programs in Africa, 
South America and the Orient 
=> Coordinated the classification of new species and description of 
species diversity in the Bajo Calima rainforest on the Pacific Coast 



Key Accomplishments and Communications: 

=> Developed a long-term cloning strategy, sold it to upper management 
and led a team implementation effort in a highly successful cloning 
program in Eucalyptus with commendation from the president 

=> Conducted research and made recommendations that led to changes 



in species and provenances with large gains in growth rate and value in 
company plantations 
=> Launched the forest nutrition research program 

=> Established tree improvement programs and new seed orchards for 
Pinus kesiya, Pinus tecunumanii and Pinus Patula 

Began an advanced education program for Colombian foresters with 
management and research potential 



Weyerhaeuser Hot Springs, Arkansas 

Tree Improvement R&D Project Leader (starting as a 

technician in Centralis, Washington) 



Roles and Responsibilities: 

=> Led a team conducting research to improve breeding, testing, selection, 

and orchard processes 
=> Led in the allocation of genetic stock for all Weyerhaeuser holdings in 

the southeastern U.S. 
=> Helped set tree improvement strategies for the southeastern U.S. 



Key Accomplishments and Communications: 

=> Conducted groundbreaking research in early genetic testing including 
the development of a juvenile-mature correlation prediction model that 
has been used by several international researchers 

=> Conducted research on plot design leading to significant changes in 
test design 

=> Conducted research on accelerated breeding of small elite breeding 
populations that led to changes in breeding population structure and 
long term strategies 

=> Conducted research that defined differences in growth and drought 
tolerance of local and nonlocal seed sources and led a team Seed 
Source Movement Review with the Senior VP of Timberiands leading 
to provenance changes and high NPV potential for 1 .5 million acres in 
Arkansas/Oklahoma 

=> Developed a unique and effective system for drought resistance 
screening of genetic stock 



RESEARCH LEADERSHIP AND 



TECHNICAL EXPERTISE 



Honors and => Adjunct Professor (current), North Carolina State University 
Awards => Adjunct Professor, Oklahoma State University, 1 984-1 986 

=> Adjunct Professor, Mississippi State University, 1984-1986 
=> Honor Societies: Alphi Chi, Xi Sigma Pi (Forestry), Gamma Sigma 
Delta (Agriculture) 

CAMCORE Lifetime Honorary Advisory Board Member 
=> CAMCORE Executive Committee Member 2000-2004 

Invited speaker at 1 9 professional conferences 
=> Weyerhaeuser Outstanding Achievement Award, 1 982 
=> Presidential Commendation For Clonal Forestry Program at Smurfit 

Carton de Colombia 
=> Forest Science article "Juvenile-Mature Correlations In Pinaceae And 

Implications For Early Selection" cited over 150 times in other scientific 

articles 

=> Nomination for and election to position of "Senior Scientific Advisor" - a 
competitive and peer-reviewed position 

Recent Weyerhaeuser Special Recognition Awards 

=> 1 998 Weyerhaeuser Special Recognition Award - for contributions to 
the Southwide Early Rotation Functional Team 

=* 1 999 Weyerhaeuser Special Recognition Award - for leading the Seed 
Source Movement Review with positive results 

=> 2000 Weyerhaeuser Special Recognition Award - for contributions to 
development of long term tree improvement strategies for Southern 
Regeneration 

=> 2001 Weyerhaeuser Special Recognition Award - for contributions to 
integrating clonal testing and review of D-fir advanced generation 
strategies 

=> 2002 Weyerhaeuser Special Recognition Award - for development of 
the electronic Seed Catalog which has facilitated key regeneration 
decisions for Southern Timberlands 

=> 2003 Nominated for Scientist 7 (Senior Research Scientist), a 
competitive position 

=> 2003 Weyerhaeuser Special Recognition Award - for contributions to 



Other 

Professional 
Activities 



developing and presenting the results that led to business buyoff in 
developing high MOE genetic stock 



Served on graduate and undergraduate committees at 4 universities 
Chairman, North American Quantitative Forest Genetics Group 

CAMCORE Technical Committee member 

Chair, IUFRO Section On Breeding and Testing Tropical Species - 

1988 to 1993 (led the organization of the international conference in 

Colombia in 1992 with over 100 participants from around the world) 

Short Course Instructor (in Spanish) - CATIE, Costa Rica - 1991 

Short Course Instructor (in Spanish) - DIRENA , Nicaragua - 1992 

Southern Forest Tree Improvement Committee - 1991 to 2001 

Chair, IEG-40 (Forest Genetics) - 1995-1997 

Review manuscripts for several scientific journals such as Forest 

Science, Forest Genetics, Theoretical and Applied Genetics, Southern 

Journal of Applied Forestry, Canadian Journal of Forestry Research, 

Silvae Genetica, Journal of Tropical Forestry 

Book chapter reviews for two books by Dr. Bruce Zobel 

Book chapter review (Forest Genetics) for Dr. Tim White 

Occasional reviewer of Forest Service and University Research 

Programs 



External 



REFERENCES 

=> Dr. Bruce Zobel, Professor Emeritus, N.C. State University 

=> Dr. Bob Kellison, Director, N.C. Forest Biotechnology Center 

=> Dr. Bill Dvorak, Professor/CAMCORE Director, N.C. State University 

=> Rafael Diaz, President, Smurfit Carton de Venezuela 

=> Dr. Tim White, Professor/CFGRP Director, University of Florida 



Internal 

(Weyerhaeuser) 



Dr. Rex McCullough (retired), Vice-President Forestry R&D 
Dr. Howard Duzan, Director, Southern Forestry R&D 
Dr. Peter Famum, Vice-President Forestry R&D 
Dr. Christine Dean, Director, Western Forestry R&D 



PUBLICATIONS 



ReCGPlt Popular => Generations of Quality: The Tree Improvement Story (Weyerhaeuser Timberlands Special 
I C O IT! Da n V) Re port) 

* K / / Tree improvement ... The Best Value For Your Money (Weyerhaeuser - Tree Growing 

Articles Times) 

=> Ideal Trees That Are Infinitely Reproducible: A Clonal Program For Eucalyptus Grandis (La 
_ r , Revista - Smurfit Carton de Colombia) 

Refereed 

Journals 23 Total 
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Evaluation of Alternative Pine Breeding Strategies 



1. OBJECTIVES: To evaluate, using stochastic simulation, the predicted genetic gains that would arise from 
pine breeding programs using a range of breeding strategies. 



2. PINE BREEDING STRATEGIES 

The genetic gains from using three alternative breeding strategies were modelled using stochastic simulation 
techniques. These techniques take account of probability and sampling variation and thus generate information 
about the likely variance of outcomes 

The five pine breeding strategies that were evaluated were as follows: 

Strategy 1. Full-Sib Breeding System (Full-sib) 

40 parents mated at random to create 80 families (four crosses per parent in a circular mating system). 

□ Parameter estimates: 

□ GCA of 40 parents, 

-* SCA. 

□ Index selection: 

Q -» I = Pindividual + biSCAf a mily + b2GCAf ema ie + b 3 GCA m ale 

Strategy 2. Polymix(PMX) Breeding System (PMX) 

□ 40 parents mated with bulked pollen from the same 40 parents. 

□ Parameter estimates: 

□ -» GCA of 40 parents as females, 
-> No SCA. 

□ Index selection: 

□ I = Pindividual + bjGCAfemale 

Strategies 3 and 4. Poly mix and Limited Pedigree Assignment ( PMXLPn ) 

□ 40 parents mated with bulked pollen from the same 40 parents, 

Stage- 1: Top 100 (PMXJLP100) or 400 (PMX_LP400) individuals selected on 

I = Pindividual + biGCAf em ale 

Stage-2: Selected individuals are fingerprinted. 

Final individual selection on: I = Pindividual + biSCAf ami i y + b 2 GCAf ema ie + baGGA^e 
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□ Parameter estimates: 

□ -» GCA of 40 parents as females, 

SomeSCA. 

Strategy 5. Polymix and Full Pedigree Assignment (PMX_FPA) 

□ 40 parents mated with bulked pollen from the same 40 parents, 
All individuals are fingerprinted, 

Individual selection on I = Pmdividuai + biSCAf amjIy + t^GCAfemaie + b 3 GCAmaie 

□ Parameter estimates: 

-» GCA of 40 parents as females and males, 
- SCA. 



3. MATERIALS AND METHODS 

Genetic and phenotypic parameters for the trait analysed (Height Growth) are contained in Table 1 . 
Table 1: Genetic Parameters for Tree Height 



Genetic parameters 


Height-' ■ 


;^Grpwthv 


Individual heritability 

Population mean 

Phenotypic variance among 
individual trees (a 2 P ) 

Coefficient of variation 

GCA additive variance (a 2 A ) 

SCA additive variance (a 2 N A) 
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3.1. Basis for Simulation 

3.1.1 Genetic Parameters: A finite locus model was used for the genetic components of 
variance. In a finite locus model, gene effects are simulated as the sum of the effects of many 
small genes, instead of as from a continuous distribution (eg Fernando. Strieker and Elston, 
1994 ). Finite locus models are well suited to modelling dominance and epistatic effects, and 
hence can easily encompass models including inbreeding depression and combining ability. 

The genome was modelled as 200 loci, each having alleles 1 and 2, each with effect vector (-a, d, a) for genotypes 
(1,1), (1,2) and (2,2). Each allele was given a frequency of 0.5 in the base population. 

Inbreeding depression was applied using the formula d = 2F^d pq , where d is the inbreeding depression and 

p and q are the allele frequencies in the population (Falconer, 1989). This measure of inbreeding is absolute 
rather than the relative measure required by the rate of inbreeding depression of 4% for every 0.1 increase in F, 
however, with p and q both set to 0 0.5, ^d = 4.8 equates to the appropriate change when the mean is 6.0. 
For 200 loci the value of d is 0.024. 

Given that d = 0.024, for allele frequencies of p = q = 0.5, the dominance variance is (2 pqd) 2 = 0.000144 for 
each locus, giving a total <j] = 0.0288 . As the total genetic variance was specified to be 0.1 12, this leaves 
a\ = 0.0832 , or Ipqa 1 = 0.000416 for each locus, giving a = 0.0288. 

The vector of allele effects then is {-a, d, a) = (-0.0288, 0.024, 0.0288). 

GCA and SCA variance are functions of F, additive and non-additive genetic variance and the selection history of 
the population (eg cross between divergent lines). As the variances have been determined above, the appropriate 
levels of SCA and GCA were obtained by simulating breeding cycles between the base population and the 
experimental population. The choice of 200 loci was also made to produce a dominance variance conducive to 
obtaining the required levels of SCA and GCA. This process is not without error. Repeated simulations with 
constant genetic parameters and population structure can produce amounts of SCA and GCA, which are quite 
different. 

3.1.2 Population Structure: A base population of 250 individuals was generated, and a second generation of 200 
individuals bred from the base through random mating. From these 200 individuals, 40 were selected on the basis 
of phenotype to be the parents of the first generation of experimental progeny. It was assumed that female parent 
was known for these individuals, and no more than four progeny from any one female parent were allowed. 

The 40 selected parents were mated according to the specified strategies, and for each strategy a crop of 7800 
progeny were produced. EBV's were estimated for these as specified for each strategy, except that for indexes 
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including SCA the phenotypes and parents of all 7800 individuals were used in an animal model to produce 
BLUP estimates of EBV. This is equivalent to using an index of own phenotype, and mean phenotypes for SCA, 
GCAfemaie and GCA^e. To evaluate the accuracy of the EBV's ten further rounds of selection were performed, 
each with 7800 progeny from 40 individuals selected from the previous round. 

Two simulations were performed for each strategy, the first with no restriction at all on inbreeding, the second 
with limits on relationship for the 40 selected parents for each crop of 7800 progeny. The inbreeding restrictions 
were: for strategies 1, 3, 4 and 5, where both parents were known, no more than one full sib and no more than four 
progeny from any parent were eligible for selection; for case 2 where only female parent was known, no more 
than four progeny from any female were eligible for selection. 

4. RESULTS 

4.1. 

Validation : To check that variance components were as expected a series of populations were simulated and 
GCA variance and SCA variance estimated. A summary of these estimates is contained in table 2. The estimated 
SCA variance is close to that desired, while the estimated GCA variance is a little high. 

Table 2: Validation of rate of inbreeding depression, GCA variance and SCA variance. [Estimates are 
means from the 200 founder populations] 





CGA variance 

(sum of male and female) 


SCA variance 


Ratio GCA/SCA 


Desired 




. o;oo83;- ^ : 


v^r •••-;-.•:/ -3:53 , , -i ' 


Realised 


0.0341 (Q.0004) : ' : , 


0.0075 (0:0003) 5 


' 4.53 ' 



4.2. Response to selection 

For each selection strategy, for each of unrestricted and restricted inbreeding, 200 replicates of the simulation 
were performed. The results presented are means and standard errors of means for both response and inbreeding. 

Figure 1 illustrates the relative gains made through the use of the 5 strategies with and without restrictions on 
inbreeding. With no restriction on inbreeding PMX_FPA and PMX_LPA400 produce almost identical response, 
and are superior to the other strategies, with PMXJLPA100 similar to Full-Sib, which is in turn superior to PMX. 
The means and standard errors in Table 3 suggest that by cycle 10, the differences between the groupings 
identified above are significant. 



Table 3. Mean (se) of tree height after 10 rounds of selection. 
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Figure 1) Gains in Tree Height (metres) over 10 rounds of selection under Five Strategies 
a) No Restrictions on Inbreeding 
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Figure 2 contains graphs of levels of inbreeding coefficients for each strategy over time. Of interest is the fact 
that the Full-Sib strategy has the highest levels of inbreeding when there is no restriction put in place. By contrast 
when the restrictions are implemented the effect is to reduce levels under the Full-Sib strategy to the lowest of all. 

Figure 2) Levels of Inbreeding over 10 rounds of Selection under five Strategies 
a) No restrictions on Inbreeding 

I — # — Full-Sib 
-h»-PMX 



PMX_LPA_100 
^~PMX_FPA 




Generation 



b) Restriction on Level of Inbreeding 




Generation 




Strictly Proprietary (Red): Disclosure strictly limited to persons on a managed list. Contact 

Proprietary (Yellow): Disclosure limited to persons confidentially bound to Weyerhaeuser on a need to know basis. 

Non -Proprietary (Green): Disclosure unlimited. 

WTC/pmxwpa efficiency tr.DOT 2/6/03 



A 



RESEARCH REPORT 
Research & Development 



Project No. Error! Reference 
source not found. 

Page 11 of 12 



Table 4 clearly shows the cumulative effect of the inbreeding without restrictions by the 10 cycle and the clear 
effect of the restrictions, whilst also revealing slightly higher levels of inbreeding under the PMX_WPA 
strategies. 



Table 4. Mean (se) Inbreeding levels after 10 rounds of selection under the five selection 
strategies. 
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Accuracy of selection of individuals and parents in the initial and final cycles of selection are shown in Tables 5 
and 6, respectively. The strategies employing parentage assignment are clearly more accurate. It is also 
noticeable that the accuracies have declined with time. This is likely due to the way in which the BLUP's and 
Indices were calculated. We did not include full pedigree information, which would have correctly accounted for 
selection. That is, by using only records and parentage of trees under consideration in any cycle, we were 
assuming "base" variance components for what were actually highly selected populations. 

Table 5. Accuracy of Selection of Individuals (Correlation between True Breeding Value and 
Estimated Breeding Value) under the five selection strategies. 

FyiSib PMX _ _ PMXJ.PA100 PMX_LPA40Q _ PMX_FPA 

Initial Cycle ^60(0^)^ 
Final Cycle §$1^1^ 

Table 6. Accuracy of Selection of Parents (Correlation between True Breeding Value and 
Estimated Breeding Value) under the five selection strategies. 

F y!LS' b j_ PMX _ p MX_ L fA 1 90 ^^Jy^°° M>U£A 

Initial Cycle p705®0g|IIf^ 
Final Cycle g^Slffifc 



5. Conclusions 

On the basis of these simulations inbreeding restrictions clearly need to be applied in order to achieve optimal 
gains under all strategies. This is dependent on the rate of inbreeding depression being as specified. An added 
benefit of applying restrictions is that the general population has lower levels of inbreeding. 

The polymix strategies are clearly the best strategies in the unrestricted method and also when inbreeding 
restrictions are applied. Of the PMX with parentage assignment strategies, the FPA gives the highest response for 
the trait under consideration. There is also the potential to increase gains even further under this strategy by 
reducing the inbreeding even further using mate allocation strategies. 
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SUMMARY 



Tree breeding practices are reviewed and related to the methods used in 
plant breeding. These methods include: mass selection (collection of 
open-pollinated seed from the best trees in forest stands), mass selection 
with progeny testing (open-pollinated seedling seed orchards), simple 
recurrent selection (untested clonal seed orchards and seed stands), 
recurrent selection for general combining ability (clonal seed orchards 
with progeny testing, second stage orchards and control-pollinated seed- 
ling seed orchards), hybrid varieties (biclonal orchards, interspecific and 
interracial hybrids), reciprocal recurrent selection, clonal selection and 
propagation, and clonal selection and propagation with clonal testing. 

Expressions are given for predicting genetic gains from different tree 
breeding methods (mainly as developed by Namkoong et a!. (1966), with 
additional expressions for methods they did not cover). 

Actual data from a 14-year-old open-pollinated progeny test of Pinus 
radiata (height 67 ft) for bole straightness (heritability = .60), and dia- 
meter b.h. (heritability = .19) are used to calculate predicted genetic 
improvement in these characteristics for different breeding methods under 
a range of operationally realistic selection intensities. (Some additional 
assumptions are made about genetic variances beyond the information 
provided by the progeny test.) The time needed to produce the first 
improved planting stock, assuming a P. radiata time scale is also calculated. 

Highest gain for straightness is predicted for a second-stage clonal 
orchard but mass production of cuttings of the best tested clones show- 
higher gains in diameter, almost as high gains in straightness and produc- 
tion of improved planting stock in 11 years as compared with 28 years. 
Progeny tested clonal orchards, biclonal orchards and control-pollinated 
seedling orchards all gave lower gains that were of the same order. The 
main reason for the much higher gains for diameter among tested clones 
bulk propagated vegetatively was the large amount of non-additive genetic 
variance assumed for this trait. Untested clonal seed orchards and open- 
pollinated seedling seed orchards showed genetic improvement that was 
much lower than that of the clonal orchard with progeny testing (about \ 
for straightness and i for diameter). The use of seed stands gave genetic 
improvement about half that of the untested clonal orchard, while the 
collection of open-pollinated seed from selected trees gave considerably 
less improvement than the seed stands. The gains predicted for the latter 
two methods however, were still appreciable. 
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Tree breeders have a tendency to consider methods of breeding and 
improvement of tree species in isolation from the general methods of 
improvement of crop plants and domestic animals. This is largely because 
of the length of time needed for most tree species to reach sexual maturity 
and then economic maturity. Thus a single cycle of breeding from initial 
selection to selection in the next generation takes a very long time and it 
is rare indeed to find discussion in tree breeding literature beyond the 
first stages of the breeding cycle. In New Zealand, exotic forests have the 
highest growth rates of the temperate regions of the world and as high or 
higher than many tropical areas. Because of this, breeding cycles are 
much shorter and an examination of long term or multi-cycle breeding 
methods has more relevance than elsewhere. 

Tree improvement work has been active during the past 15 to 20 years 
in New Zealand. Initial emphasis was largely on selection and testing of 
species and provenances, and later an improvement programme for 
P. radiata was initiated with the aim of providing adequate supplies of 
improved seed from clonal seed orchards. The objective was to provide 
all future requirements of seed, which also would have the moderate 
degree of genetic improvement which can be expected from intensive 
field selection only. With the further development of breeding work in 
this species involving the realisation of extra genetic improvement by 
progeny testing, as well as the initiation of new breeding programmes 
with other species, it has become necessary to make a general evaluation 
of all possible breeding methods and their underlying theory so that the 
method most suitable for each situation can be adopted and so that future 
developments are soundly based. This pattern of development in tree 
improvement is found with minor variations in most countries where 
active tree improvement programmes exist. 

This review examines methods of crop breeding as described by Allard 
(1960), and relates them to methods of tree breeding currently practised 
or advocated. A review of the comparative theoretical study by Namkoong 
et aL (1966) indicates the composition of genetic gain for some tree 
breeding methods and expressions are developed for several others. 
Finally, actual data from 14-year-old P. radiata are then used to make 
numerical comparisons of genetic gains in two traits of contrasting 
heritability using the different methods under optimum practicable 
intensities of selection. 
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RELATION OF PLANT BREEDING METHODS, 
TO THOSE OF TREE BREEDING 



In a consideration of breeding methods currently in use with crop 
plants and their application to tree breeding, the framework used bv Allard 
(I960) is most suitable. From the point of view of the breeder species can 
be divided into two groups; those which are mainly self-pollinatin* and 
those mainly cross-pollinating. This distinction determines what breeding 
methods are available for the species, as those applicable to the self- 
po Iinators are distinct from those that can be used with cross-pollinators 
All cross-pollinating plants (unless they have been purposely inbred) are 
highly heterozygous and inbreeding of these will usually result in loss of 
vigour and other adverse effects. Self-pollinated plants consist of collections 
of independently existing, homozygous inbred lines which show no loss of 
vigour in this state. The goals of breeding programmes in the two tvpes 
involve the production of vigorous, fairly uniform, heterozvgous offspring 
of the cross-pollinators and homozygous pure lines of the self-pollinators' 
Where asexual reproduction such as by vegetative propagation is possible 
additional breeding procedures can be used. The reproductive biology of 
the species usually determines its particular breeding procedure." For 
instance, the ability to tolerate some generations of selfine by maize a 
a naturally out-crossing species, has made possible the use o~f inbred lines 
and hybridisation. 

The main- difference between the material tree breeders and croo 
breeders have to deal with is that in most cases crop plants have undergone 
countless generations of selection so that the species is quite different in 
its characteristics and population structure from the original wild popula- 
tion. Again, because crop plants are frequently annuals, the numbers of 
generations of selection that have been undergone in the past and that 
ran be undergone in the future are vastly greater than with tree species 
Tree breeding is also a new science that has only been widely applied in 
the last two decades. Usually, tree breeders are dealing with unimproved 
species, frequently as natural populations. When they arc grown as exotics 
as restricted portions of natural populations, thev sometimes become 
hybridised with other populations in the process of introduction. 

Among the species of interest for tree breeding, all conifers are out- 
crossing, essentially monoecious, wind-pollinated species with separate 
male and female flowers. Of other species of interest, most are out-crossin* 
frequently wind-pollinated, sometimes insect-pollinated. There are few 
species of tree breeding interest that are largely self-nollinated so this 
group will be omitted from further consideration. 

For out-crossing species of crop plants, every plant is likely to be 
heterozygous and this heterozygosity must be maintained during the 
breeding programme or restored at the end. In crop plants there is a wide 
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variety of agencies and mechanisms controlling cross-pollination and also 
a big variation in the ease with which this can be controlled by the breeder. 
With tree species, for conifers the collection of pollen and isolation of 
female flowers are in most cases quite simple operations. Controlled 
pollination in eucalypts involving emasculation of bisexual capsulate 
flowers is more difficult but nevertheless quite feasible. Perhaps the 
greatest practical problems of controlled pollination involve access to the 
crowns of mature trees, and the long periods between pollination and seed 
harvest, e.g. of over two years in P. radiata. However, as grafting or 
production of rooted cuttings enables the mature tree crown to be propa- 
gated nearer ground level, and the propagules normally have earlier sexual 
maturity than seedlings, a partial solution to this difficulty is available. 
Nevertheless, controlled mating of trees is an expensive, time-consuming 
and laborious procedure which, when it must be used, adds several years 
to the breeding cycle. 

The most important breeding methods used with cross-pollinating crop 
species, enumerated by Allard {Joe. c/7.), are: 

(i) Mass selection 

(ii) Backcross breeding 

(iii) Hybridisation of inbred lines to form hybrid varieties 

(iv) Recurrent selection 

Each of these will now be discussed and their application to tree breed- 
ing indicated. In addition vegetative propagation can be used to multiply 
sefected trees from natural populations, or to multiply any product ot 
other breeding methods, and its role in tree improvement will also be 
discussed. 

(i) Mass Selection 

With this method desirable individual plants are chosen and their 
open-pollinated seed collected and bulked without a progeny test to 
produce the following generation. Selection is thus based on female 
parent phenotype only, as there is no control over pollination and thus 
male parentage of the seed collected. This corresponds, in the forestry 
context, to collecting open-pollinated seed from the best trees in the forest 
and using this to establish further plantations from which, in turn, the 
best trees are selected and so on 1 . It should be pointed out that collection 
of seed from silviculturally thinned stands involves an element of recurrent 
selection (see below) improving male and female parentage. The effective- 
ness of this procedure is determined by the heritability of the trait involved 
and the intensity of selection. If heritability is high, accurate identification 
of the genotype of the female parent results and the offspring are appre- 
ciably improved. The purpose is to progressively increase the proportion 
of superior genotypes in the population over several generations and thus 
to raise the average population value for the trait. 

Mass selection has been effective in modifying traits in crop plants that 
are easily seen or measured. It has also been effective in developing varieties 

i Mnf* that certain tree breeders term the selection of superior trees and propagation of these in orchards 
Is procedure male and female parentage of the mulling seed .s select. 



for special purposes and in changing the adaptation of varieties to fit them 
for new sites. With trees one could expect that mass selection for mor- 
phological traits such as bole straightness and branching habit would be 
effective. Similarly the adaptation of an exotic species would be assisted 
by some generations of mass selection. In actual fact, mass selection for 
bole straightness, branch morphology and vigour has proved effective in 
Queensland with P. elliottii van elliottii (McWilliam and Florence, 1955). 
Selection of the best trees in the forest and collection of their open- 
pollinated seed resulted in an increase in acceptable stems per acre from 
112 from routine seed, to 206; there was a similar increase from 4 to 23 
"plus" stems per acre. There is evidence from provenance trials of several 
exotic species in New Zealand, e.g. P. rcdiata, that a few generations of 
natural selection in its new habitat with some mild silvicultural selection, 
is sufficient to improve the general performance and particularly survival 
and hardiness of the species considerably. 

Mass selection has proved rather ineffective, at least in the short term, in 
modifying yields of well-adapted varieties or species of crop plants. This 
has presumably been due in the main to the low heritability of yield which 
would result largely through a progressive reduction in additive genetic 
variation in the face of previous generations of selection, either natural or 
artificial. Other causes would include uncontrolled pollination so that 
male parents are random, not select, and the fact that intense selection 
can reduce population size to the point that inbreeding occurs with 
consequent loss of vigour. Tree size is a character for which competition 
has imposed countless generations of natural selection, reducing genetic 
variability. Variation between populations, provenances or species is 
likely to be much greater and should be exploited. 

Mass selection combined with progeny testing can help overcome some 
of the deficiencies of simple mass selection. In this case in crop plants 
open-pollinated progenies are grown from a portion of the seed of all 
individual selections and only the seed which gives rise to the best pro- 
genies is bulked and sown to form the following generation. This pro- 
cedure bears some resemblance to the practice in tree breeding of estab- 
lishing an "open-pollinated seedling seed orchard 1 ' (Wright and Bull, 
1963; Shelbourne, 1962). Here, open-pollinated progenies of select trees 
are grown in a progeny test and later, on the basis of the measurements in 
the test, only the best trees in the best families are retained and seed is 
then collected from them. Such a process involves family selection 
(removal of poorer families) and in addition some within-family selection 
involving removal of all but the best individuals within families. The 
latter procedure cannot normally be utilised in mass selection with progeny 
testing of crop plants. 

This procedure is also somewhat similar to the crop-breeding procedure 
of line breeding. Here, after several cycles of mass selection and progeny 
testing, the open-pollinated seed of a number of the best progeny-tested 
plants is bulked and planted in an isolated plot where random mating is 
allowed. Annual re-establishment of such a plot provides a continuing 
source of improved seed. Excessive inbreeding is avoided by ensuring 
that an adequate number of unrelated lines are involved. 
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It is proposed that a more practicable and economical method of 
establishing open-pollinated seedling orchards than has been previously 
suggested is to bulk open-pollinated seed of select trees, establish the 
orchard at close spacings and then thin selectively on the basis of pheno- 
type alone. Such a method will enable large acreages of orchard to be 
cheaply established and, as will be shown later, can be expected to give 
nearly as much improvement as the progeny test type of orchard. To avoid 
risks of inbreeding in offspring of the orchard,' large numbers (200 or 
more) of open-pollinated families should be used. 

(ii) Backcross Breeding 

The backcross method is particularly suitable in crop plants for trans- 
ferring specific genes to a good variety which is deficient in one or a few 
characteristics. In this method recurrent backcrosses to the more desirable 
parent are made while selection is practised for- the characters being 
transferred from the donor parent. It is important that sufficient plants 
of the heterozygous recurrent parent variety are used so that its character- 
istic features will be maintained in the backcross-improved products. 

The prospects for the use of backcross breeding are rather limited with 
tree species mainly because a single sexual generation takes periods 
varying from about 4 to 20 years. However, in species with early sexual 
maturity, provided the trait to be transferred can be assessed at a very 
young age, backcross breeding might be applicable. The traits most 
suitable for transfer are those with high heritability and high economic 
importance such as disease resistance. If disease resistance is easy to 
transfer it may be controlled by a single locus and in turn will be subject 
to mutations of the pathogen. 

(iii) Hybrid Varieties 

Hybrid varieties depend for their superiority on the hybrid vigour or 
heterosis associated sometimes with F t hybrids. Inbred lines, clones, races, 
varieties or even species can be crossed to achieve the heterotic effect! 
Once a suitable hybrid combination has been identified the parental 
genotypes must be maintained without change so that the same hybrid 
can be produced repeatedly. This has limited the use of hybrid varieties 
in crop breeding to those species that can be seJfed repeatedly to give 
homozygous inbred lines and also those species that can be propagated 
vegetatively. Homozygosity also gives rise to parental lines of consistent 
known combining ability, and because gametes are genetically quite 
uniform, gives hybrid offspring of considerable uniformity. 

The greatest development of hybrid varieties has been with maize (corn) 
whose biology is well suited to this method of breeding. Desirable parents 
are selected in open-pollinated populations, selfed through several 
generations to produce homozygous inbred lines and finally "these lines 
are crossed. In practice the yield of seed from a cross of two inbreds is 
so poor that usually the best single-crosses, as those are called, are them- 
selves crossed to form a double-cross (or four-way cross) which is the 
commercially produced hybrid. It is important to realise that if this 
procedure is followed, with no selection at any stage and random crossing 
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of inbreds, the hybridisation will merely restore the vigour lost by selfing. 
The improvements in grain yields of over 20% that have been achieved 
have resulted because of selection of the original parents, selection between 
and within inbred lines and finally intense selection between hybrids. The 
latter is the critical operation that is responsible for the big improvement 
in yield associated with corn hybrids. Of the approximately 100,000 inbred 
lines developed by 1951 (Kiesselbach 1951) only about 60 have proved 
suitable in certain combinations, from which to develop commercial 
hybrids. Selection for double-cross yields by single-cross performance 
has been found to be quite satisfactory. 

There is little possibility of using hybrid varieties of this type in tree 
breeding. Any adoption of the same techniques of crossing inbred lines 
used with corn would be impossibly slow, e.g. involving 42 years to take 
P. radiata through four generations of selfing, even supposing this is 
possible for more than a very few genotypes. In addition, the performance 
of hybrid varieties of this type is often quite specific to a given environ- 
ment ; and in forestry we must breed varieties adaptable to the uncontrolled 
conditions found in a forest. 

Heterosis, however, could be utilised in three other possible ways with 
tree species: interspecific hybridisation, interracial hybridisation and 
utilisation of specific combining ability between individuals within a 
population. 

Interspecific crossing is being used in a limited number of tree improve- 
ment programmes as a means of achieving genetic improvement not 
available by other methods. This usually takes the form of combining 
useful features of one species with those of another, and results in general 
intermediacy of characteristics for the hybrid, as well as hybrid vigour. 
In Korea the P. rigicla x P. taeda cross provides a vigorous hybrid that 
performs better than either parent species (Hyun 1963); in California the 
cross P. radiata x P. attenuata extends the range of the P. radiata parent 
on to higher altitude, colder sites while preserving its valuable features 
(Stockwell and Righter, 1946); in Europe the larch hybrid Larix decidua 
x L. leptolepis is probably the best example in forestry of a well-proven 
heterotic interspecific hybrid. As a means of achieving a commercially 
useful hybrid involving production of hybrid seed on a large scale, inter- 
specific hybridisation in tree breeding is not generally very useful because 
of low viable seed yields, and difficulties in obtaining good yields of 
exclusively hybrid seed. Hybridising orchards can be used to produce F, 
seed by planting selected clones of one species in a matrix of selected 
individuals of the other species, which act as pollen parents. The resulting 
seed will always contain some of the pure species and there are difficulties 
in matching flowering phenology of the clones of the two species. Inter- 
specific hybridisation is more promising if it is associated with vegetative 
propagation by rooted cuttings. By using this technique the best individuals 
from interspecific crosses can be bulk-propagated from hybrid offspring, 
as is commonly done in poplar breeding. 

Inter-racial crossing, that is, crossing of different races or populations 
of the same species, might be expected theoretically to provide hybrids 
with better balanced genetic constitutions than interspecific hybrids while 
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giving high yields of fertile seed. Provided the two races were genetically 
divergent, the resulting cross would be heterozygous and there would be 

a good chance of obtaining hybrid vigour. Few of such crosses are yet ? 

old enough to evaluate and experimentation with this method is con- jL 

tinuing. The races chosen for hybridisation should be well adapted to the | 

site on which the hybrid will be grown or alternatively, should be from | 

two opposite extremes of native environment to the intended site. A hybrid t 

of the latter type is likely to prove more physiologically versatile, though £ 

the former would be better adapted. | 

The use within a population of specific crosses between individuals which * 

show high specific combining ability for vigour or any other characteristic £ 

is a method of possible application in tree breeding, though at the present * 

stage of development this is unlikely to be adopted until after several |. 

generations of mass selection or recurrent selection (see below). As in I 

the use of inbred lines with corn, large numbers of "hybrid" or specific | 

combinations must be tested to identify the best crosses for the traits in jf- 

question and which also show other features of desirable quality. The ; 

success of such a programme will depend on the existence of appreciable \ 

amounts of dominance genetic variance for a characteristic in the popula- > 

tion in question. Very extensive crossing and testing programmes will be ; 

required to identify the heterotic combinations, and there is the further j 

problem of mass producing the seed of a particular cross. Two-clone [ 

seed orchards might be one solution but it is not known yet how serious * 

the effects of selfing would be under these conditions. Artificial pollination [ 

may be the only practicable answer. Because forestry demands genetic [ 

material that is quite variable as an insurance against catastrophe, several [ 

crosses would have to be produced and their seed bulked for planting. ; 

Specific crosses could not be used as a continued cyclic method of \ 
breeding as selection and crossing within the offspring of such two-clone 

orchards would result in sib mating and inbreeding. If a further generation f 

of two-clone orchards were initiated, it would be necessary to return to t 

an improved population derived from recurrent selection, mass selection, £ 

or recurrent selection for General Combining Ability. An alternative pro- i 

cedure which has the advantage of allowing a continued and cumulative £ 

breeding programme to exploit heterosis is discussed below (reciprocal | 
recurrent selection). 



(iv) Recurrent Selection 

There are four types of recurrent selection, depending on the way in 
which individuals with desirable traits are identified. These are simple 
recurrent selection (identification on parent phenotype), recurrent selection 
for general combining ability (on average offspring performance from 
several testers), recurrent selection for specific combining ability (on off- 
spring performance after mating with a single tester) and reciprocal 
recurrent selection. 

Historically, the ideas and theory of recurrent selection were developed 
in response to the realisation that in corn breeding selecting within and 
between inbred lines and finally between different hybrid combinations 
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to find the best heterotic hybrids was an inefficient procedure. Raising the 
frequency of desirable gene combinations in the parent population couid 
instead be done cumulatively, at the same time maintaining genetic 
variability and allowing favourable gene frequencies to be increased by 
genetic recombination. 

Simple recurrent selection 

Desirable genotypes are selected on the basis of their own phenotype 
and are intercrossed without restriction to produce populations for the 
second cycle of selection. This process is therefore different from mass 
selection in that both male and female parentage is select. Additional 
cycles of selection can be carried out as long as improvement continues. 
For given intensities of selection the genetic improvement expected will be 
double that for mass selection. 

In tree breeding, the selection of plus trees, their propagation in a 
clonal seed orchard and the harvest of the seed represent a single cycle of 
simple recurrent selection. The clonal propagules are allowed to inter- 
pollinate randomly in the seed orchard and the seed produced is used to 
establish plantations from which the next cycle of selection is made. As 
no test crosses are made, simple recurrent selection, like mass selection, 
is only efficient for traits with heritability high enough to make accurate 
selection for a superior genotype from phenotypic evaluation. The other 
types of recurrent selection are more suitable for selecting for low herit- 
ability traits as they involve test crossing to measure combining ability. 

Recurrent selection for general combining ability (RSIGCA) 

This involves mating phenotypically selected individuals with a hetero- 
zygous tester stock of broad genetic base. In tree breeding practice this 
may involve: — 

(i) crossing selections with a standard group of tester individuals 

(ii) polycrossing select individuals with a pollen mix of several select 
individuals 

(Hi) collecting seed from open pollination of the select individual with 
the unselected population 

(iv) using vegetative ramets of the select individual. (This provides 
adequate screening of individuals for GCA provided non-additive 
genetic effects are not very large.) 

Measurement of the offspring when planted in a suitable field design 
will provide half-sib, full-sib, or clonal family averages which give esti- 
mates of general combining ability. The individuals with the highest GCA 
are then propagated (in corn by selfed seed or in trees by grafts or cuttings) 
and intercrossed. The seed from this intercrossing is bulked and used to 
form the population for the second cycle of selection. 

This procedure corresponds to the use of a clonal seed orchard with 
progeny testing that is in general use in tree improvement programmes. 
Superior trees are selected on the basis of their own phenotype in the 
forest, and propagated vegetatively in an orchard or clonal archive. Seeds 
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are either collected from the select parent trees from open pollination, or 
alternatively from controlled crossing or polycrossing in the orchard to 
generate half sib families (polycrossing) or a group of full sib families, 
e.g. North Carolina Design IL, 4 tester design (Roberds et al. 1967) from 
which GCA can be estimated. "Families 1 ' of clonal ramets can also be 
used for this purpose (Libby, 1964) though this has rarely been attempted 
with trees. 

Where GCA is estimated from open pollination of the original parent 
trees in the forest, it must be assumed that the number of pollen parents 
is infinite, and that for different parent trees located in different parts of 
the forest the same pollen spectrum is being sampled. The larger the 
number of male parents that pollinate the female flowers on the parent 
tree, the better the chance that the pollen will be a good genetic sample 
' of the population. The use of open-pollinated seed is thus liable to induce 
bias in estimation of GCA because of differences in pollen spectra and 
numbers of male parents represented. This method cannot be safely used 
where appreciable genetic differentiation of local populations in natural 
stands has occurred. It is worth noting that collection of quite large 
numbers of fruits from different parts of the tree, and for serotinous pine 
species, from different years' pollination will alleviate this problem to 
some extent. Similarly, collection of open-pollinated seed from single 
trees not growing in forest stands should be avoided. The big advantage 
of using open-pollinated testing is that the seed may be collected in the 
same year that the tree is selected and the progeny test can be planted 
about the same time as the clonal orchard. 

Where polycrossing is used to generate half-sib families, a standard 
pollen mix of 10 or more clones can be used to pollinate all clones, or 
alternatively different pollen mixes can be used to pollinate different 
groups of clones. Where individuals in the pollen mix are being poly- 
crossed their pollen should be left out of the mix. 

Where full-sib test crossing is used, any number of tester males can be ' 
used though a minimum is four. This method provides information on 
specific combining ability as well but is the most costly in terms of con- 
trolled pollination and size of test (four or more families per select clone 
tested). Both polycrossing and test crossing usually involve the grafting 
of ramets of seiect clones and a wait of several years for these to produce 
male and female flowers and then control-pollinated seed. 

If groups of ramets of the select clones are used to screen for general 
combining ability, because of no genetic variation within clones or between 
ramets of the same clone, relatively fewer cuttings are needed to estimate 
the clonal means than are seedlings to estimate half-sib or full-sib family 
means. Also the clonal test can be established at the same time as the 
orchard or could even be combined with it. If topophytic or physiological 
age effects are important with the rooted cuttings however, clonal testing 
for GCA would not be appropriate, and the gain from clonal testing would 
be much reduced. 

The best portion of the original select clones can be retained in the 
orchard with the remainder being removed, or alternatively a new orchard 
can be established with these clones only. A point that is frequently over- 
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looked is that to achieve an appreciable increase in genetic improvement 
a rather large proportion of the clones must be rejected. The seed collected 
from the reconstituted orchard will produce stock for establishing planta- 
tions and these trees when sufficiently mature will form the populations 
for the next cycle of seleccion. 

Recurrent selection for general combining ability represents the main 
breeding procedure that tree breeders are adopting or at least working 
towards. The very slow progress in comparison with agricultural crops 
is due firstly to the long sexual generation and secondly to the long 
"assessment* generation^ the time taken before meaningful evaluations 
of progeny can be made, which is usually the longer of the two. Prospects 
for reducing both these generation times are quite good. Some success 
has been achieved in accelerating initiation of flowering in some species, 
and as data on juvenile-mature correlations are built up, earlier evaluations 
of progeny will be possible for some characteristics. 

Another procedure that comes under the designation "recurrent selec- 
tion for GCA" is the controlled pollinated seedling seed orchard. Pheno- 
typic selection is made in the forest in the normal way and controlled 
polycrossing is carried out either on the original parents, or on their . 
grafts growing in a clonal archive a few years later. The polycross progenies 
are then planted in a progeny test which is later converted to a seed 
orchard after assessment of the progenies, by removing all but the best 
families and the poorer trees from the best families. In practice for silvi- 
culture! reasons, selection intensity within the progeny test is limited to 
a total of about 5%, or about 25% for family selection and 20% for 
within-family selection. An extra stage of selection is thus involved 
beyond that of the clonal orchard with progeny testing. Because of the 
operational difficulties of this method, low family and within-family 
selection intensities, as weil as the possibility for considerable sib-mating 
in this orchard resulting in inbreeding, the relatively high gain expectation 
must be weighed carefully against the other disadvantages. Control 
polycross seed can be bulked, as was suggested for the open-pollinated 
seedling orchard and the test area thinned on the basis of individual 
phenotype alone. This again has many practical advantages and reduces 
gain only slightly, as will be shown later. 

A further" method of recurrent selection for GCA involves similar 
stages of selection to the control pollinated seedling orchard but these 
are more intense. Normal phenotypic selection is followed by clonal 
orchard establishment and matings are made to generate full-sib families 
in some design such as a partial diallel or 4 tester. Identification of the 
clones of highest GCA from the best crosses is made in the progeny test 
and the best individuals in these families are selected and propagated to 
form a "Second-Stage" clonal orchard as proposed by Roberds et ai 
(1968). The main requirement for satisfactory operation of this method 
is a large range of original phenotypic selections among which extensive 
crossing to generate a large number of unrelated families has been carried 
out. Improvement is then due to initial phenotypic selection, selection 
for general combining ability (family selection) and finally a quite intensive 
within-family selection. 
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Recurrent selection for specific combining ability (RS/SCA) 
^ This was proposed by Hull (1945). Test crosses were to be made to a 
single homozygous tester to develop complementary strains that would 
achieve maximum heterozygosity and thus heterosis in hybrid combination. 
Repeated cycles of family selection on the basis of test cross performance 
would mean however that if some defect became apparent in the original 
tester much effort would have been wasted. In addition, this method 
would produce improved populations of very narrow genetic base and it 
is not considered satisfactory for tree improvement. 

Reciprocal recurrent selection (RRS) 

This is the fourth and final type of recurrent selection and was proposed 
by Comstock, Robinson, and Harvey (1949) as a procedure useful in 
selecting simultaneously for both GCA and SCA. Two heterozygous 
source populations (A and B) are required and should be as genetically 
divergent as possible, though not sexually isolated, and should have 
generally desirable growth and technological characteristics. A group of 
individuals from population A is selected and mated to a similar sample 
from population B. In the same way a select group from population B is 
test-crossed to population A. The individuals in each population which 
show highest general combining ability with the other population, based 
on the performance of the test crosses, are later mated to each other 
within populations and the resulting offspring of each population are 
used to make the second cycle of selection. It is necessary to use enough 
plants in each generation for intense reselection after the progeny tests, 
and to leave a sufficient number after this to avoid serious inbreeding. 
The populations developed by several cycles of RRS could be crossed at 
any stage and, depending on the amount of divergence in gene frequencies 
and thus, heterozygosity achieved, the resulting hybrids will be more or 
less heterotic. This scheme will build up two populations that are high in 
GCA for the characteristics selected for, and, if inter-racial heterosis for 
characteristics such as yield docs exist, these populations will become 
progressively more divergent in their gene frequencies so that maximum 
.heterosis will result when they are crossed. In addition, the hybrid will 
be generally intermediate in characteristics between the two populations. 

Comstock et al (1949) compared theoretically the efficiency of this 
method with recurrent selection for general combining ability and recurrent 
selection for specific combining ability under three levels of dominance- 
partial, full, and overdominance. They concluded: 

(0 if dominance is partial, RS/GCA and RRS are about equally 
efficient and both are superior to RS/SCA 

(ii) if dominance is complete, the three methods are essentially equal 

(iii) if overdominance exists, RRS and RS/SCA are equal and both are 
superior to RS/GCA. 

These comparisons assumed no epistasis, two alleles per locus and linkage 
equilibrium. The presence of epistasis (non-allelic gene interactions) 
would favour RRS and RS/SCA and linkage disequilibrium would also 
favour RRS. This system of breeding shows some possible applications 
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in tree improvement. In programmes working with fast-grown exotics, 
two races (provenances) of a species each with equally good performance 
but from distant parts of the natural species range can be combined 
into a heterotic hybrid or a hybrid that will tolerate the contrasting 
environmental conditions of each race. Such a scheme would take no 
longer than a conventional one for production of progeny tested seed, 
but would increase the number of selections necessary as two populations 
must be maintained with an adequate genetic base. This would involve 
the following steps when applied to a tree species: 

1. Select 200 best phenotypes each from provenances A and B. Propagate 
by cuttings or grafts in randomly mixed clonal archive, with each 
population isolated from the other. Also propagate all 400 clones in 
mixed hybridising orchard at close spacings. 

2. Make controlled polycrosses on all clones using pollen mix of sample 
from population A on population B and vice versa. Collect seed, grow 
stock and plant in progeny test. 

3. After 10 years or so evaluate all progenies and remove the poorer 80% 
of clones from each population archive. Allow interpollination within 
each population and collect open-pollinated seed to establish com- 
mercial plantations of each population. Remove the poorer 80% of 
clones from the hybridising orchard, and allow interpollination and 
collect mixed hybrid seed. 

4. After 10 years or so make second cycle of selection of best phenotypes 
in commercial plantations of each population, as in 1 above and repeat 
test crossing, etc. Continue cycles of selection, test crossing, etc. 
Natural hybridisation of the two populations in the hybridising orchard 

can occur at any time. Alternatively, manual pollination of clones in the 
single population archives by pollen collected from the other population 
would result in production of pure hybrid seed. Furthermore, provided 
the species involved was commercially propagatable by cuttings, any 
superior hybrid at any point in the scheme could be multiplied by this 
method (see below). 

This scheme or variants of it could be initiated with any species provided 
two desirable but genetically divergent provenances existed. Little extra 
work need be involved as total selection intensities within each population 
could be adjusted to the same as would be used in a single population. 
The main prerequisites for such a scheme to be initiated would be the 
availability of stands of two genetically divergent populations (proven- 
ances) of sufficiently desirable nature to make improvement within these 
a satisfactory proposition. In addition, some data from preliminary 
crosses between such populations giving indications of heterosis for yield 
and satisfactory quality of the hybrids would be needed. 

(t) Vegetative Propagation and Clonal Selection 

The use of vegetative propagation to produce commercial planting stock 
is not strictly a method of breeding, though it is normally associated with 
selection and achieves a considerable improvement in the species con- 
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cerned. Phenotypic selection combined with clonal testing with vegetative 
propagules and ultimate release of tested clones is a procedure bv which 
considerable improvement of a species can be achieved within a single 
sexual generation. It is thus a useful adjunct to anv breeding procedure 
in that clonal testing will enable precise evaluation of the genotyoes to 
be made, and if propagation techniques are economic, will enable any 
superior genotype to be replicated indefinitely for commercial use 

Extensive use is made of this approach in horticulture and agriculture 
with high value crops like potatoes, strawberries and fruit trees which are 
all dependent on vegetatively propagated, single genotypes. Possible 
applications of this technique in forestry have been limited by the ability 
to root cuttings of tree species on a commerciallv economic scale and 
particularly to root cuttings from trees old enoush'amons which to make 
meaningful phenotypic selections. Forest tree plantina stock, unlike fruit 
trees, must be mass produced cheaply. An additional use for rooted 
cuttings is in establishing clonal tests for general combinins ability; even 
when rooted cuttings cannot be produced easily and cheaofv this method 
of testing has many advantages over using families, but some disadvantages 
as well Also, in conjunction with other genetic experiments, clonal tests 
can help in estimating non-additive genetic variance. 

Rather few forest tree species have been propagated commerciallv bv 
rooted cuttings. Poplars and poplar hybrids are an outstanding example 
of a tree genus that has probably reached a higher level of improvement 
than any other forest tree species. Rooted cuttincs of Crvptomeria japonica 
of certain selected clones have been widely planted in Japan, and Saiix spp 
are usually multiplied by cuttings. Recent work on P. radiata (Thulin and 
Faulds, 1968) indicates that commercial propagation of this species bv 
rooted cuttings is feasible, and in a further paper (Thulin, 1969) the 
operation of an integrated programme for the production of progeny 
tested seed and commercial production of cuttings of tested clones is 
described. 

Once commercially applicable techniques of rooting cuttings of forest 
tree species have been developed, phenotypic selection and clonal testin- 
of any superior individual, the product of any breedinz method be ft 
backcross breeding, interspecific hybridisation or even natural or induced 
mutation 1 can be multiplied and utilised by commercial forestry. Perhaps 
the greatest advantage of using clonal plantations at the present sta»e of 
development of forest tree breeding is the ability to effect more rapid" and 
greater improvement m genetic quality of planting stock than could be 
achieved by any other breeding method during the first cycle of a breeding 
programme. It should be emphasised, however, that provided the whole 
population has been effectively examined, the gain from cycles of selection 
is not cumulative as in simple recurrent selection; further gains can only 
be made by selection from an improved population derived from selection 
and genetic recombination. 

1 After an initial enthusiasm for mutation breedim fthe induction nf -u _ ■ 

artificially to produce grater phenotypic variation). r^^S^f^Su^SSSS^A^SSSS 
choosrae to make use of generally abundant existing variability rathe "than to £5 k 'I"'^ 
mutations, mostly deleterious, that were produced by Radiation. "* d °" na P hl » rd 
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EXPRESSIONS FOR PREDICTING GENETIC GAINS 



la the preceding review of plant breeding methods applicable to cross- 
pollinated species it was attempted to relate these to the tree breeding 
procedures at present in use. Comparisons between different procedures 
were of a qualitative nature. It is possible, provided estimates of certain 
population parameters are available, to make comparisons of predicted 
genetic gains from different breeding methods in quantitative genetic 
Terms. The best method of comparison would be by an economic analysis 
in terms of money invested in breeding on a per acre basis, and the 
expected increase in income per acre at the end of the rotation by the 
improvement in quantity and quality of wood. This type of analysis is 
well beyond our present knowledge. At best, in planning a breeding 
programme at present we can try to maximise our predicted genetic gain 
for a given expenditure and this involves adjusting the amount of work 
and thus expenditure involved in different phases of the programme in 
relation to the amount of genetic gain predicted, e.g. the costs of selection 
against the costs of progeny testing. 

Several authors have attempted to make theoretical comparisons of 
different tree breeding methods, e.g. Wright (1962), Wright and Bull 
(be. c//.), Stern and Hattemer (1964), Johnssen (1964), Libby (1964, 1965), 
Namkoong (1965), Namkoong et al. (1966) ? and Squillace et al. (1967). 
It is unfortunate that some of these attempts contain errors and in- 
consistencies, e.g. Wright, Wright and Bull, and Stern and Hattemer. By 
far the best and most comprehensive account is by Namkoong et al. (1966). 
They review heritability and gain concepts, and derive expressions for 
predicting genetic gains under different methods of breeding. Their 
expressions "for gain" prediction will be reviewed and expressions will be 
given for genetic gain prediction for additional breeding procedures not 
covered in their paper. 

To estimate genetic gain, provided a properly constructed and precisely 
estimated heritability is available for the trait in question, it is only 
necessary to determine the selection differential. Breeding programmes in 
forestry commonly involved several stages of selection, each having its own 
selection differential and heritability. Selection is usually on the basis of 
a fraction of the population saved, and so the selection differential is 
measured indirectly. If a normal distribution of values for a characteristic 
under selection, i.e. values from individual trees or family averages exists 
with a variance a\ y the difference between the mean of the select group 
and the population mean, i.e. the selection differential, can be derived to 
be i<x where i is the selection intensity. The values of i for various popu- 
lation sizes are given by Fisher and Yates (1953, Table 20), Falconer 
(I960, pp. 193-4), and Nanson (1967). Values of i for a population of 
greater than 50 are also shown in Appendix Table 1. 
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TABLE 1. SUMMARY OF PLANT BREEDING METHODS 
WITH EQUIVALENT TREE BREEDING PRACTICES 



NETIC GAINS 
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Plant Breeding 



Tree Breeding 



1 . Mass selection 



2. Mass selection with 
progeny testing 

3. Simple recurrent 
selection 

4. Recurrent selection for 
general combining 
ability 



5. Reciprocal recurrent 
selection 

6. Hybrid varieties 
(utilisation of 
heterosis) 



7. Selection and mass 
vegetative propagation 

3. Selection, and mass 
vegetative propagation 
of tested clones 



Collection of open-pollinated seed from selected trees, 
usually at relatively low selection intensities (0.1-5 % of 
stand selected). 
* Open-pollinated seedling seed orchard (also involves 
additional within half-sib family selection), 
(i) Clonal seed orchard, 
(ii) Seed stand (seed production area). 

(i) Clonal seed orchard with progeny or clonal testing 
followed by roguing or establishment of new 
orchard. 

(ii) Control-pollinated (polycross) seedling orchard 
(also involves additional within half-sib family 
selection). 

(iii) Second stage clonal seed orchard (selection of best 
individuals from the highest GCA clones' families 
in progeny test, with establishment of clonal 
orchard). 

Reciprocal recurrent selection — not yet in use but has 
possible application in combining RS (GCA) with 
interracial hybridisation. 

(i) No real equivalent currently in use in tree breeding 
(with inbred lines). 

(ii) Heterosis may be utilised by: 

(a) interspecific hybrids 

(b) inter-racial hybrids 

(c) specific combining ability — biclonal orchards 
Selection and afforestation with rooted cuttings — 
limited use so far with poplars, willows etc. 
Selection and afforestation with rooted cuttings of best 
tested clones. 



Predicted genetic gain is given by. the expression 
AG — ia u b 

where b = heritability for the particular selection system 1 



(1) 



fa* 



Thus: AG = 

Where *i = 
k 



additive genetic variance 

fraction of the total additive genetic variance in the covari- 
ance of additive values for the particular relatives in question, 
e.g. k = I for selection amongst half-sib families. 

In the preceding review, several breeding methods having application 
in forestry were discussed and a summary of these with the equivalent 
tree breeding practices is given in Table 1. The expressions for expected 
genetic gain for these methods will now be compared. 

l b has been referred to as the heritability for the particular selection system (Namkoong ft ai. 1966). 

<7\ , 

However a more customary formulation of heritability is h 2 » — j where <7^ is the phenotvoic variance of 
the individual tree. O* u 
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(1) Mass Selection 

Mass selection involves the selection of the best trees in forest stands 
and the collection of their open-pollinated seed. 

Predicted gain, AG =*ii^ (2) 
where k — i (one parent only is selected) 

erf = phenotypic variance of the forest stand 

2 i 2 < 2 2 
= CT w + *p T T °"G 

<t* = the environmental error variance of the small plot or 
jnicrosite. The total within-plot variance a\ contains a\. 
plus the remaining genetic variation depending on the family 
structure of the experiment. 



o-p = The error variance of large plots, or the failure of genotypes 
to behave in the same way between macrosites (commonly 
estimated in experiments as a family x replication inter- 
action), 

o-^ = The genotype x environmental or macrosite interaction, or 
error variance caused by the failure of genotypes to behave 
in relatively the same way in different environments. 

o% = Total genetic variance. 



(2) Mass Selection and Progeny Testing 

Mass selection and progeny testing corresponds most closely to the 
open-pollinated seedling seed orchard. This involves collecting open- 
pollinated seed from select trees and planting half-sib progenies in a test 
design which is later thinned to the best trees in the best families. 

JG = i t — - -r i2— r 13-— W> 

a t <r 2 <7 3 

The gain i 2 — — is due to selection of the best families and the element 

i 3 -^- corresponds to additional selection of the best trees within the best 

families. i lt i 2 and i 3 are selection intensities determined by proportions of 
individuals in the whole population, of families amongst those originally 
selected and of individuals within families that were saved. 



20 



= phenotypic variance in population, as equation (2) 
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fL + d + ^2£ + J al (Variance of half-sib family means) 
nre re e 

a 2 + a i + ff £ J jj. i- <x5 (within half-sib family variance) 

"? + <* + *i 

number of trees per plot 

number of replications/location 

number of locations 



(75 = dominance genetic variance 

The double prime on aj; indicates that this variance represents that in 
the population after selection, mating and thus genetic recombination 
have occurred. It can be assumed that cr-j = o£ , as genetic variance 
will be regenerated after recombination, oj. however represents the 
additive genetic variance left in the restricted, selected portion of the 
original population. It is less than a\ by an amount depending on 
heritability and intensity of selection. 

a 2 # = (j\ (1 — £/), where J5 = a\!<r { and v is a variable dependent 
on i r and is tabulated by Finney (1956) and shown in Appendix, Table 2. 

Alternatively open-pollinated seed from plus trees may be bulked, the 
offspring planted at close spacing, and thinned later on the basis of 
individual phenotype only. 



(4) 



Where: i 4 = combined intensity of i , and i 3 above; i.e. the proportion 
left after thinning. 



(3) Simple Recurrent Selection 

Simple recurrent selection is equivalent to a clonal seed orchard without 
progeny testing, or to a seed stand or seed production area where an 
existing stand has been thinned to the best trees and seed is collected from 
these. In both cases random interpoilination of selected parents produces 
seed. 



(5) 



where a\ = phenotypic variance as in equation (2). 
Note: k = 1 



21 



(4) Recurrent Selection for General Combining Ability 

This method involves phenotypic selection (as in (5) above) followed by 
either reselection amongst select individuals on their average test cross or 
polycross progeny performance (clonal seed orchard with progeny test 
reselection), or the selection of new select parents from the highest GCA 
families in the test. The latter procedure includes the control-pollinated 
seedling seed orchard, and the second stage clonal seed orchard where 
the new select parents are propagated vegetatively. 

(a) Clonal seed orchard with reselection by: 

(/) Control-pollinated North Carolina Design II test 

GCA can be estimated by making control-pollinated test crosses of each 
select clone to a standard series of testers (e.g. Roberds et al., 1967). 
Predicted gain is given by: 

2 ^ 2 

= ^-7- 2i,^ (6) 

where \ x and i : represent selection intensities of individual and progeny 
test selection. 

ctj = phenotypic variance as in equation (2). 

■» 2 ** *• 

2 __ g 7 j_ g P , g ge j_ gmf ■ °A' 

° 2 nre 1 re 1 e ' 4 1 4 

nre ' re 1 e ' 16 ' 4 

<rj f = male x female interaction 

-i»S 

a 2 

The gain is thus in two parts; i 4 — from phenotypic selection, and 

2 i 2 — — from progeny test reselection. It should be noted that the single 

prime on o\ refers to reduced genetic variance in the restricted portion 
of the original population that is being progeny-tested, and that the gain 
from progeny testing is multiplied by two as both future parents of seed 
orchard seed are being reselected (on the basis of their progeny per- 
formance). 

(ii) Polycross or open-pollinated test 

GCA can also be estimated by controlled pollination of all select 
clones using a standard pollen mix of some of them. If the number of 
clones in the pollen mix is large, the only component of equation (6) 
which will change will be cr 2 . Tbe component of a| is divided by 

1 Note that when considering genetic variances within and between families, cpistatic interactions are 
disregarded. 
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the number of tester clones; this component can be assumed to be zero 
if more than 10 testers are used. 

(hi) Clonal test for CCA 
GCA can also be estimated using vegetative propagules of the select 
parent such as rooted cuttings, as Libby (1964) indicated. Libby showed 
that the additional genetic gain from clonal selection is always con- 
siderably greater than that expected from full-sib or half-sib family 
selection. He compared the ratios of family selection to individual selection 
gains for half-sib, full-sib, and clonal "families". This is not the same as 
a progeny test situation where initial phenotypic selection is followed by 
subsequent reselection on the basis of a progeny test. Libby also assumed 
that topophytic or physiological age effects on the growth of the cuttinss 
were absent. Assuming that there is perfect correlation between gene 
effects of seedlings and cuttings: 



Ci <?<> 



(7) 



where a] = phenotypic variance as in equation (2) 



,2 = 



- ™ T — + T + (*V * 



nre 



re 



A) + 4 



A'* 



and <rp = additive, dominance and epistatic variance in select 
portion of population = <jg, 

0c = variance of "c" effects which are the effects due to 
vegetative propagation that result in non-equiva- 
lence of clonal variance, a^, and total genetic 
variance, a£. This can best be illustrated by the 
well-known dffference in morphology of young 
cuttings or grafts of pines from young seedlings" 
Propagation has caused an effect of increasing the 
apparent age of the cutting. This effect varies from 
clone to clone with a variance 

(b) Control-pollinated seedling seed orchard 

A further procedure advocated in tree breeding that is a variant of 
"Recurrent selection for GCA" is the control-pollinated seedling seed 
orchard. After phenotypic selection of plus trees, controlled polycrossing 
using a mixture of select tree pollen is carried out either on the parent 
tree (usually impracticable) or on its vegetative propagules. 

The polycross half-sib families are planted in a test that is later thinned 
on the basis of measurements to leave the best trees in the best families. 



jG-i^+i 



(8) 



where a u <r z and <x 3 are the same as in equation (3) (open-pollinated 
seedling orchard). 
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Note that the gain from phenotypic individual tree selection is twice that 
in the open-pollinated seedling orchard. 
The three components of the gain realised from the products of the 

2 

orchard [thus represent phenotypic selection it — , a half-sib family selec- 



tion for GCA u 



and within family selection i 3 ^ A * 



Each stage of 



°2 " cr 3 

selection has its own selection intensity, of which i, and i 3 are limited by 
the silvicultural problems of thinning the test area. 

Alternatively, as in the open-pollinated seedling orchard, polycross 
seed may be bulked, and the progenies thinned on the basis of phenotype 
only. 

In this case: 



(9) 



see equation 4 for a A ?.., cr 4 , etc. 



A further variant of this procedure involves instead of thinning the 
test, propagating the select individuals in a new "second stage" clonal 
orchard. The expression for expected gain would be the same as equation 
(8) if polycrossing were used and would be as follows where full-sib 
families were involved: 



2 12 12 

Ju = lj r i: r lj 

<T l <7 2 <T 5 



(10) 



In comparison with the control-pollinated orchard, this provides scope 
for increasing i> and i 3 at the expense of extending the time needed to 
produce improved seed. 

Expressions for a\ and aj shown for equation (3) will be modified to 
each include i <t£ and ^ <r£. respectively. The possibilities for intense 
family and within family selection are only limited by the number of 
families in the test and the number of trees per family. 



(5) Utilisation of Specific Combining Ability — Biclonal Orchard 

Perhaps the easiest way for tree breeders to make use of heterosis that 
has been a foundation for the hybrid varieties of plant breeding, is to 
exploit specific combining ability shown by individual cross combinations. 
The testing of considerable numbers of controlled crosses, arranged in a 
design to sample some of the possible combinations between the groups 
of phenotypically selected trees, could be followed by the setting up of 
several two-clone orchards. 



(H) 



where i x = intensity of individual phenotypic selection 
i 2 = intensity of selection among full-sib families 
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a\ = phenotypic variance as in equation (2) 

2 " nre * re ' e ^ 4 ' 2 

There are thus three elements of gain ; the first due to individual phenotypic 
selection is the same as in other procedures; the second is that due to 
selection of the best full-sib families, and the third element is again based 
on the full-sib family but is dependent on dominance variance. If this 
variance is large, the extra gain from this approach will be considerable. 



(6) Reciprocal Recurrent Selection 

This method involving parallel phenotypic selection in two populations 
with parallel test crossing of select individuals in each population to a 
select sample from the other population, allows the cumulative improve- 
ment and moulding' of two populations so that they have increasingly 
divergent gene frequencies and thus increased heterosis on crossing at 
subsequent stages in the programme. An expression for the genetic gain 
achieved by this system cannot be given. It will be made up in the same 
way as equation (11) except that non-additive gene effects are likely to 
be much greater when chromosome sets of widely differing origin are 
brought together. 
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(7) Selection and Mass Vegetative Propagation 

Phenotypic selection followed by bulk vegetative propagation of select 
tree planting stock can provide considerable genetic improvement. At the 
present stage of development in tree breeding its most important applica- 
tion is phenotypic selection in unimproved forest stands. 

A complicating factor with many species is that the process of propaga- 
tion maintains the adult characteristics of the ortet in the propagule. For 
instance, in many pines cuttings or grafts of older trees show the mor- 
phology of the upper crown of the parent tret. This effect is represented 
by an additional constant, c, in the prediction of genetic gain thus: 



(12) 



where a c = °a + a D + G i. 

= total genetic variance 
<j\ = epistatic variance 

c = average propagation effect 1 in using rooted cuttings rather 
than seedlings 

= phenotypic variance as in equation (2) 

i Different gene effects may be involved in the control of "adult" cutting versus "juvenile" seedling traits- 
Tbe present approach to the problem is to ignore this and assume that vegetative propagation of physio, 
logically "older" material merely produces an average modifying effect. 
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Depending on the size of the non-additive genetic variance (a^ and aj), 
gains will be greater than from simple recurrent selection. It is probable 
that for certain species the average propagation effect, c, may be greater 
than the genetic gain. It could be in a negative direction cancelling the 
effects of genetic improvement, though there are indications that at least 
in pine species the effects are supplementary to genetic gain for most 
morphological traits. 

(8) Selection, and Mass Vegetative Propagation of Best Tested Clones 

If clonal tests of rooted cuttings are established after phenotypic 
selection, and a further reseiection is made among the select clones on 
their test performance, an additional gain can be obtained. 



o> 4- cry) 

AG = c + u- 5 + u 



(13) 



where i 2 = selection intensity among phenotypically select clones 
ag. = reduced total genetic variance after phenotypic selection 



o% — 



the variance of "c ' effects (these may be permanent like 
physiological age effects in pines or transitory due to size 
and condition of cutting material). 



a 

nre 



p 
re 



It can be seen that the variance of clonal means (a?) is much reduced, 
firstly by the process of multiplication of each clone and secondly by the 
reduction in a*, .to <r£. With high broad sense heritabilities of, e.g. .80, 
this variance could be reduced to about one quarter of its original size 
with corresponding overriding importance of the "c" effect variance. 
However it is not known at present what size <j% is likely to be. Where 
broad sense heritabilities are high (over .80) and where a b is small, clonal 
testing is unlikely to be worthwhile. However where heritabilities are 
lower and/or is higher, then clonal testing will most probably be 
effective in substantially increasing genetic gain. Where tt is planned to 
only use limited numbers of the best clones for commercial propagation, 
clonal testing is an essential safeguard. 
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COMPARISON OF PREDICTED GAINS 
USING RADIATA PINE DATA 

Introduction 

It is possible to compare predicted genetic gains for any species from 
the different breeding methods discussed previously provided heritabilities 
(broad and narrow sense) have been estimated for the trait under com- 
parison. Actual variance data from Pinus radiata grown in New Zealand 
is now used to make comparisons of predicted gains under a range of 
selection intensities that are operationally practicable. 

Namkoong et ai (1966) made comparisons of genetic gains using 
simulated variances for traits of low and moderate narrow sense heritability 
(.05 and .20) and conditions of effective and ineffective testing, i.e. where 
the reducible variance in the progeny test is maximally or minimally 
reduced. These comparisons were made for five breeding methods: clonal 
orchard (without progeny test), clonal orchard with progeny test, open- 
pollinated seedling orchard, control-pollinated seedling orchard and 
two-clone orchard. They also considered the time taken to achieve gains ; 
and tabulated operational values per year, representing the total units of 
gain accrued throughout the breeding cycle divided by number of years 
in the cycle. To make gain from different methods mathematically com- 
parable, they restricted~all methods to a total selection intensity corres- 
ponding to 1 in 1000. 

In terms of total gain realised, they found that the control-pollinated 
seedling orchard was generally just the best under low and moderate 
heritability and under effective and ineffective testing, though the 
superiority over the gains by the clonal orchard with progeny testing was 
generally slight. When operational value per year was considered, the 
clonal orchard with progeny testing was just superior to the control- 
pollinated seedling orchard. The two-clone orchard, except under high 
heritability and highly effective testing, was generally slightly inferior to 
or about equal with the former two methods. Interestingly, total gains 
from the untested clonal orchards were much lower, but about the same 
as the open-pollinated seedling orchard. Under effective testing the open- 
pollinated orchard showed a considerably higher operational value per 
year than the untested clonal orchard. 

Two traits of radiata pine, stem straightness and breast height diameter, 
with contrasting heritabilities, will be used here to compare gains from a 
wider range of breeding methods, including those involving rooted 
cuttings. 

Methods 

Estimates of additive genetic variance, plot-to-plot and within-plot 
variance for diameter at breast height and bole straightness (subjectively 
rated 0 to 9) were available from an experiment involving 29 open- 
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poUinated families of 14-year-old P. radiata (M. H. Bannister, manuscript 
in preparation). There were no estimates for dominance variance (assumed 
zero for straightness and equal to additive variance for diameter) or for 
genotype x site interaction (assumed zero). The variance of "c" effects or 
propagation effects, namely the effect of vegetative propagation on the 
cutting's growth and morphology, was assumed to be one quarter ot 
within-plot variance for both traits. Heritability (narrow-sense) was .60 
for straishtness and .19 for diameter; heritabilities (broad-sense) were 
.60 and ~ 37 respectively. (Sec Appendix for details of variances and 
calculations.) 

Predicted genetic gains, expressed as a percentage of the population 
mean before selection, were calculated using these parameters for different 
breeding methods with varying intensities of selection. When methods 
involved phenotypic and subsequent genotypic selection (based on test 
performance), a phenotypic selection intensity corresponding to the best 
tree in 1,000 was used, and the intensities of progeny test reselection were 
varied from 1 in 2 to 1 in 100. * t 

It must be recognised that some of the extrapolations involved in 
prediction of improvement in diameter growth and to some extent straight- 
ness may break down in practice owing to the effects of intertree com- 
petition and to non-linearity of the subjective rating scale for straightness. 
This should not affect the ranking of gains for different methods and the 
effects of heritability on this ranking. 

Results and Discussion 

For the different breeding methods shown, genetic gains are expressed 
as percentages of the population means for straightness and diameter in 
Table 2*. The greatest aains in both highand low heritability traits (straight- 
ness and diameter) were predicted for method 6, the second stage clonal 
orchard (87% and 27° 0 ) and for method 10, clonal propagation and 
testing (76%and45%). However, usinaP. radiate as a basis of comparison, 
the stock from method 6 could only be planted 28 years after the initial 
phenotypic selection (see "Time to first improved planting stock", Table 
2) while that from method 10 could be planted after 1 i years. To provide 
sufficient unrelated families in which to make selections for the second 
sta°e orchard (Method 6) it requires a very large number of crosses with 
several different series of testers or a large number of individuals m the 
polycross pollen mix. 

Gains in straightness of 70% from Method 5, the clonal seed orchard 
with progeny testing, were about the same as for the control-pollinated 
seedling orchard, Method 7a and b (71% and 74%), and only a little lower 
than Method 10. The sain of 45% for diameter with clonal propagation 
and testing is much higher than for any other method mainly because this 
method utilises the considerable dominance variance assumed for this 

^An examination of test gains relative to total gains in Table 2 reveals 
that any form of testing results in relatively more additional gain with a 
low than with a high heritability trait. Also, where dominance or epistatic 

• On pages 30-31. 
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variance is appreciable (for diameter), clonal propagation is able to 
IS th total genetic variance, in contrast to other methods where only 
Jftwo-clone orchard is capable of utilising part of the dominance 
variance (reciprocal recurrent selection, not covered in this comparison, 

^^should beemphasised that the gains from clonal selection, propaga- 
tion and testing do not take into account the average propagation effect c, 
TA^^ complement or detract from genetic gam. However the 
donai test enables the net effect of gain and Propagation effort to be 
observed and the best performing clones can be identified and mass- 

•^r-KS* 0 *-"- «« "»v roseny 

(Method 5) depends very much on how progeny test reselection is applied. 
F^«-poLated test (Method 5b) was established at the same .time 
or earlier than the original orchard and then a new orchard established 
o the best clones later, it would be 22 years from initial selection on the 
P„Za time scale before the first progeny-tested p -^S 5t ^ r ^ S r 
available. If a control-pollinated test was used it would be 27 years or 
™« If however the first orchard was close-spaced, e.g. 6 x 6 feet and 
^ "pSTOWely on the basis of test data to 20 x 20 feet or wider 
progeny tested stock could be produced after Id years. This would be 
expensive if grafts rather than cuttings were used. 

Afferent progeny test methods (5a, b. c) do not materially affect gams 
from £ Jested orchards. Full-sib designs, e.g. Method 5a are more 
«,stW and more than one family per clone is necessary All control- 
S itS desiens including polycrossing are costly and delay test estab- 
Kent. Open-pollination is cheap, involves no delays but is less precise 
n Satins several combining ability (GCA). The use ot rooted cuttings 
to sc £n foVGCA is an interesting alternative to polycross and open- 
pollinated stock that allows early test establishment and tar fewer pian , 
Eer tested clone due to low within-cione variation, \\hen non-additne 
variance and topophysis or propagation effects are large, however, it 

w toc sxss — -» di " d ^ 7i) shows ab ° ut ■a 

d rains as the proseny tested clonal orchard; gains are composed 
SL S ly ^o^thi for the second stage clonal orchard because 
sSSLzl difficulties in thinning the test area limit the intensity of 
familv and within-family selection. ...... , u 

If however, polycross families are balked and indiv.dua phenotyp.c 
select^ Tonfy isVrried out in the bulked stand (Method 7b). the genetic 
gains for both high and low heritability traits are not seriously reduced. 
Such an orchard would be very much cheaper and eas.er to plant and 
than the control-pollinated seedling orchard. The mam requ.rernent 
for such an orchard would be a large number of clones selected (e g. .00) 
and the use of several pollen mixes with 10- or more pollens each to avoid 
the possibility of many individuals in the resulting orchards being related. 

TwcSlone orchards (Method 8) differ little in expected gams trom 
progeny-tested clonal orchards. For diameter where dominance variance 
was assumed equal to additive variance, genetic gam was predicted to be 
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23%, the same as for the clonal orchard and polycross test (Method 5b). 
Gains for straightness were lower than 5b as might be expected. 

Untested clonal orchards 1 show gains much lower for both straightness 
and diameter than tested orchards, e.g. 11% compared with 22% for 
diameter, underlining the value of testing and reselection. 

The open-pollinated seedling orchard, generally regarded as a low gain 
method, shows predicted gains as high as an untested clonal orchard (47 % 
and 13% for straightness and diameter). This method could be useful 
where a low cost, simple improvement programme is needed, or where 
serious difficulties are encountered with vegetative propagation, provided 
that the species flowers early in its life-cycle and that characters selected 
for can be assessed early. As for the control-pollinated orchard there are 
many practical advantages and little if any loss of gain in relying entirely 
on phenotypic selection (Method 3b) instead of between and within 
family selection for the open-pollinated seedling orchard. Open-pollinated 
seed is usually available in quantity from 200 or more select trees and when 
bulked permits large acreages of orchard, planted at close spacings and 
later thinned selectively, to be established very cheaply. Predicted gains 
of 50% for straightness and 12% for diameter show little difference from 
gains from Method 3a. 

The use of rooted cuttings of untested select trees (Method 9) shows 
gains similar to the untested orchard for straightness and as high as the 
tested clonal orchard for diameter. This method has the great merit that 
rapidly increasing quantities of improved stock are available within a 
year of initial selection. With 500 clones selected, six years would be 
necessary to build up to an annual production of 300,000 rooted cuttings 
in P. radlata (Thulin 1969). It should of course be emphasised that con- 
tinuing gains from this method and Method 10 require the provision of 
improved populations of seedlings for selection by methods 1 to .8. 

As interim methods of tree breeding, gains from seed stands and seed 
collection from selected trees in the forest show that really worthwhile 
improvement can be made, e.g. by restricting seed collection to the best 
tree in 20, a gain of 15% in straightness is predicted compared with 43% 
from an untested clonal orchard; or a seed stand selectively thinned to 
10% of original stocking could realise a 25% improvement in straightness. 
The important features of these methods are their low cost and early 
production of improved stock. 

Something frequently not fully appreciated by tree breeders is the 
logarithmic relationship between selection intensity (proportion saved) 
and genetic gain. The eains from the clonal orchard (Method 4) show that 
to increase gain in diameter from 11.0% to 12.9%, an increase of 1.9% 
(contrast possible gain of 24% shown in Table 1 with progeny testing), 
the work and thus cost of plus tree selection must be increased tenfold. 
The breeder should therefore know the relative costs of selection and of 
establishing and assessing progeny or clonal tests so that by adjusting 
the proportion of expenditure on each he can maximise his genetic gain 
for a given outlay. 

. Aq actual improvement of 63% in straishtnew and 8% in diameter w*i reaped for this method in a 
P. radiate control- pollinated progeny test aged 10 years, heiKht 50 feet. (Shelbourne 1969.) 
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These numerical comparisons of predicted gains are of limited value in 
a comparison of breeding methods. Other important factors are the length 
of the period before improved stock is available, the relative costs of 
different methods and stages of each method and the practical realities 
and problems of a particular situation. Generally the breeder is well 
advised to attempt maximum gains in any cycle, as it is the duration of 
these cycles that limit gains compared to those achieved by crop breeders. 
The high gains shown by methods of clonal breeding using rooted cuttings 
suggest that every effort should be made to solve propagation problems 
and implement these methods. In addition when initiating a programme 
with a new species, the limited gains that can be obtained immediately by 
mass selection, that is collection of open-pollinated seed from non- 
intensively selected trees should be exploited to the full whilst more 
sophisticated breeding methods are coming to fruition. 
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Appendix 



CALCULATION OF GENETIC GAINS 



Example — Diameter 

The following data (provided by M. H. Bannister) were derived from 
pollinated families of i\ radiaia grown near Nelson. 



J a- =.1483 /. oz = .5932 



» 2 = a 2 
t w 



= 3.000 
a 2 = 0.0310 

a 2 = 3.1793; al = 1.783 
Assume: a£ = .5932 
= 0 



"I 



i al = .4905 



heritabiltty (narrow sense) = .187 
heritability (broad sense) = .374 

x » 10.17 = population mean diameter. 



Assume that all progeny and clonal tests involved 10 replications of 10 tree plot 
one location, and that tr g - (genotype x environment interaction) is zero. 



1. Collection of Open-Pollinated Seed from Selected Trees 

AG = i i • -i 
a, X 

e.g. i = 2.06; AG - 2.06 (0.5) .5932 = 3 y 
* ' (1.783)(10.17) ' /Q 



2. Seed Stand 



^JG = i — 
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MS 



derived from 29 open- 



cations of 10 tree plots at 
m) is zero. "~ 



3. Open-Pollinated Seedling Seed Orchard 
(a) Thinned progeny test 



5932 . 

^ - = .187 = narrow sense heritabiUty 



2 

ff A' 



*a' 



*' from Finney (1957) = .9322 
: .4898 



2 2 
100 ' 10 



A' 



= .0300 + .0031 -r .1225 = .1556 



a, = .3945 



= 3.000 -f .0310 
a 3 = 1.741 
Assume n = 3.37; i : = 1.27; ij = 1.40 



/3 37 (0.5) .5932 , 1.27 (.1225) ^ 1.40(.4449) \ 

* G = \ n783 " .3945 1.741 ) 

t 

= (.5606 .3944 -f .3578) • 



10.17 



= 1.3124/10.17 = 12.9% 
(b) Thinned bulked progenies 

AC f'J°* ' 1 

/3. 37 (0.5) .5932 2.06 (.5674) \ 1 
= V T783 ' 1.776 J ' 10.17 

- (.5602 + .6581) • = 1.2183/10.17 = 12.0 J 

because " i a A" + i a A' 
= .4449 + .1225 
« .5674 

ff A " °A ? " = ,0258 

*J « * J - .0258 

- 3.1535 
* 4 = 1.776 . 
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4. Clonal Seed Orchard (without progeny test) 

<d i 

o { x 

e.g. i = 2.67 

AQ= 2-67 (.5932) = 8 ?% 
1.783 (10.17) /0 



5. Clonal Seed Orchard (with progeny testing) 
(a) 4 tester N.C. Design II control-pollinated test 



J G = fi 1 °i + 2i z i!l'Vl 

" loo- 10 4 4 

Where a^ f = J a% = .5932, 4 = .1483) 
2 = 3.00 -.0371 ,0310 ^ 122 ._ 

= .1923 
a 2 = .4385 

i 2 varies from 0.80 to 2.06. Assume i t = 3.37 

/ 3.37(.5932) ^ .80 (2)M225) \ J_ 
V 1.783 ' .4385 / ' x 

=* (1.1202 + .4470) • = 1.6170/10.17 = 15.9% 



(b) Polycross or open~poUitiQted test 

c| = ^ + !h + ^: 

2 100 10 4 

3.000 , .0310 , 17? , 
= 100 10" ^ 

= .1556 

02 « .3945 

e.g. i 2 = 0.80 

_ / 3.37(.5932) ^ .80 (2) f.!225) \ 1 
* U \ 1.783 ' .3945 / 10.17 

= (1.1202 + .4968) • jqjj = 1.6170/10.17 - 15.9% 
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(c) Clonal test for general combining ability 

\ C t Ql J X 

2 2 2 2 2 
<V = °t ~ * a A ~ a D 

= 1.9619 

a\ = 1.9619,4 = .4905 

100 10 v y 

_ 1.9619 , .OMO 
100 ^ 10 

= 1.492S 

<n = 1.213 



f .4893 + .4898 + .4905 



= SO AO = A37(.5932) ^_ .S0(.4S9S) \ 1 
* * V 1.783 1.213 J 10.17 



(1.1202 -r .3230) • Y^p; = 1.4432/10.17 = 1- 



6. Second Stage Clonal Seed Orchard 

As. 



\ ff[ ct> aj J X 



/ 3.37 (.5932) _ 2.06 (.2450) ^ 2.665 (.2966 A 1 

"A 1.7W .422S 1.653 / 10.17 

- (1.1202 + 1.1937 + .4782) • 

= 2.7921/10.17 - 27% 



7. Control-Pollinated Seedling Seed Orchard 
(a) Thinned progeny test 



qi (as for open-pollinated orchard) = .3945 
ff3 (as for open-pollinated orchard) = 1.741 
Assume: i r (at 1 in 1000) = 3.37 
i 2 (at 1 in 4) = 1.27 
i 3 (atlin 5) = 1.40 

A37(.5932) ^ 1.27 (.1225) ^ 1.40(.4449) \ 1 
JG = V T783 ' .3945 1.741 ) ' 10.17 

= (1.1202 + .3944 + .3578) • y-j^ = 1.87:4/10.17 = 18.- 

(Continued overleaf) 
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(b) Tlxinned bulked progenies 



As for open-pollinated orchard except initial phenotypic selection gain is doubled. 



-( 



3.37 (.5932) ^ 2.06 (.5674) \ 1 

J. 783 ' 1.776 J 10.17 

■ (U2Q2 + . 6581)^ 

•17.5% 



8. Two-Clone Orchard 

2 2 . i 2 . a 2 
°m ff w - * a A' ~ * °D 

= 1.9619 4- .2449 -r .4449 
= 2.6517 

2 2 2 2 

C2 "ioo T To^T + 4 

= .0265 -i- .0031 + .2449 -r .1483 
= .4228 
ff2 = .6502 

JG = ^1.1202 -p—^ . 65Q2 ) 101? 

= (1.1202 + .3014 + .1825) * = 1.6041/10.17 = 15.8% 



9. Clonal Selection and Bulk Propagation 
AG = 



where a G = <x A -f »5 4* of 

« .5932 + .5932 4- 0 
« L1864 



1 

10.17 



= c + 2.2404/10.17 = 22.0%-r c 

where c = average difference between seedlings and cuttings due to effect of propaga- 
tion. 
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ielection gain is doubled. 



10. Clonal Selection, Clonal Testing and Bulk Propagation 



^. = 0^(1-^0 
jS == -3 = .3732 



1.1864(1 - .373 (.9322)) 
> .7737 
.4905 

2 2 
ff w . °c i - i 2 

100 10 ° c 

,0310 ?737 4905 
100 10 



= 1.2869 
^ = 1.134 



= c + (2.2404 -r .8918) • 



1 



10.17 



= 3.1322/10.17 = 30.8% 



Variance Data for Bole Straightness 
The following data were derived from the same source as those on diameter: 



10.17 = 15.8% 



Assume: = 



a} = I a\ « 3.8418 


/. ai = 15.37 


2 2 . i 2 , 


afj = 20.262 


^ = 1.3880 




af = 25.466 




aj = 5.047 




2 -> _ *) 
°D = a I = U; a C = 


ia^ = 2.1843; 



heritability (narrow sense) = hcritabib'ty (broad sense) = .60; 

x = population mean straightness 
= 21.52 



lue to effect of propaga- 
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Appendix 



TABLE 2 

Extract from Table 1 from Finney, D. J,, The consequences of selection for a variate 
subject to errors of measurement. Revue Inst. Int. de Stat. 24: I ,-3 1956. 



p 




p 




Proportion 




Proportion 




of population 




of population 




selected 




selected 




.95 


.1903 


.25 


.7583 


.90 


.2879 


.20 


.7813 


.85 


.3595 


.15 


.8051 


.80 


.4169 


.10 


.8308 


.75 


.4653 


.075 


.8453 


.70 


.5071 


.05 


.8619 


.65 


.5442 


.025 


.8833 


.60 


.5777 


.02 


.8888 


.55 


.6083 


.0125 


.S989 


.50 


.6366 


.01 


.9031 


.45 


.6631 


.005 


.9142 


.40 


.6881 


.001 


.9322 


.35 


.7121 


.0005 


.94 (approx.) 


.30 


.7354 


.0001 


.95 (approx.) 











where a\> 



■ reduced additive variance remaining in a group of selected plus 
trees. 



a A = additive variance in unselected (original) population. 
p — -~ = narrow sense heritability. 



o~ = phenotypic variance. 

v' is tabulated above depending on the proportion of the population 
selected. 
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INTRODUCTION 



™. QOCC rMKir STRATEGY FOR THE NORTH CAROLINA STATE UNIVERSITY- 

,^ 1 S55SS to Movement program 

Steven 6. McKeand and Floyd E. Bridgwater' 
Attract -A strategy for the Cooperative's third-cycle Breeding program for 
lo b,olly »s A ,S U was «~»^^ S 
s „ ip-*-- each 

also be possible. Our strategy wium _™ " ..j.,,,,,, wi n consist of 1 60 selections pgz 
a, , different level of '7^ f ^^ , ^™^ pll |,ponl in . given geographic | ™ S 

^^££Ztt£^** S, « MW 't. nS rnos°, t P 

Sn n .s T .rr,rT..ch, primarily to ^^^SZ^SSL ? hT,°ht E S i! 
intensively selected and managed hierarchy ^Jbeth. «~ „ „ ovlde ^ 3 » 

S e,ec,=d group of trees A Ji hie „, chv „i„ g 2 | 

maximum short-term genetic gam TO * " oreserve genotypes with extreme 

be the genetic "^^^^^ as an 

breeding values for individual tra^ (not ° . ja cha in future 

insurance ™^°°J^ 
population of trees. 

The improved efficiency of this proposed breeding strategy alonfl , with the 
me ,m P ru comoared to our current program, will result m a 

reduct.on in populat. on «es ""JJJJ ooperators. The increase in selection 
substantial reduction in «J«^n*^ coope ^ ^ ^ 

SIS S^^^^ -~ — 

will be maintained by judicious management of all three h.erarch,es. 

swords : Elite breeding populations, genetic gain, mating designs. Pinus taeda L. 
testing designs 



5 
> 



CO 

o 

o 

GL 
CD 



As the North Carolina State University - »^ C w ^ 

progressed toward the third cycle °< members. A 

effective breeding strategy to ensure both short- and longi erm MMna ^tn. 

primary criterion for the breeding strategy was t ha it * flex We °^urre t ^ ^ 
Decisions made in the third-cycle program " 1 by the 

Cooperative's working area, and options ava table to the diverse 
methods employed in the next cycle. In addition the str ategy * * ^CSSSand objectives of 

objectives of the membership. While the members share m sup P° rt ^ level 

the Cooperative, there is a diversity in product ^/"^TI^ should be 
that must be considered. Additionally, as new ,nformat.on becomes ava.labie, 
flexible enough to incorporate the information. 



• Geneticist and Associate Professor. Cooperative Tree ^^^^^SSl 
Geneticist and Professor. USDA Forest Service. Box 8002. Department of Forestry, n 
State University. Raleigh. NC 27695-8002. USA. 
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Third-cycle breeding will commence around 1 997 when 30-40% of the third-cycle selections 
are available. A summary of the breeding strategy adopted by the Cooperative is given in this pap er . 
More details are available in the report prepared for Cooperative members (Anonymous 1992). 



MAINLINE POPULATION 



The size of the mainline breeding population is the driving force behind many breeding strategy 
decisions. The appropriate breeding population size depends upon the trade-offs of short-term versus 
long-term gain. Short-term gain is maximized by intensive selection and, hence, small effective 
population size. Long-term gains are maximized by milder selection and larger population sizes. 
Population size basically affects the risk of losing favorable alleles through genetic drift, and this is a 
more serious problem for those alleles that are at a low frequency. The loss of alleles affects the 
progress possible in breeding programs. The problem of choosing population size basically "reduces 
to one of deciding on the minimum effective population size that will avoid the major effects of drift, 
that is, that will retain most of the genetic potential in the population" (Rawlings 1970). 

Many simulation studies as well as long-term selection experiments indicate that effective 
population sizes (N E ) 50 are sufficient f or maintaining low frequency favorable alleles in a breeding 
population. Long-term gain should be possible for many generations provided that N £ ^50 is 
maintained (Namkoong et al. 1980). For favorable alleles at frequency ■» 0.10, a population size of 
N £ = 30±. gives a 95% probability that these alleles will be fixed at frequency = 1 .0 at reasonable 
selection differentials (i.e. 1/100) (Rawlings 1970). Alleles at very low frequency (.01) would require 
a much larger population size with similar selection intensity to have a high probability of fixation (Kang 
1979, Namkoong et al. 1988, Robertson 1960). Even when many of the simplifying assumptions 
(such as equal gene frequencies and only additive genetic effects) are relaxed (e.g. Baker and Curnow 
1969, Mahalovich 1990), the genetic gain from populations with N € < 100 are very comparable to 
gains from much larger populations. 

For alleles that are not under selection (but may be valuable for changing selection goals in 
future generations), the population size necessary to maintain these alleles with a high probability is 
higher than for those alleles under selection (Kang 1 979). Again, the appropriate N c depends upon the 
gene frequency and the desired probability of not losing the allele. 

Namkoong et al. (1988) also emphasized the benefits of maintaining a relatively large 
population size in the first generations of breeding. Low frequency alleles will be selected and wi« 
move into the intermediate frequencies where future selection will be more effective in smaller 
populations. Fortunately, the Cooperative managed very large breeding populations in the first two 
cycles. Each cooperator had access to an average of 600+ selections for a given geographic area • 
over 4000 selections throughout the Cooperative. 

Given the long-term nature of the Cooperative's breeding program and uncertainties abort 
future selection pressures (Kang and Nienstaedt 1987), a third-cycle breeding population with an,nit ~ 
effective population size of N c = 160 selections available for each cooperator (i.e. a recruitme 
population) is reasonable. The traditional approach of keeping the breeding population as large ast£ 
breeder can afford is appropriate, and 1 00 to 200 selections is sufficient to yield long-term gain, wn* 
this population size is much smaller than the 600^ selections currently available to each cooper*^ 
the reduction to 160 selections is warranted since 1) The mainline breeding populations 
composed of local selections only. For future generations, selections from adjoining or even 
geographic regions may be combined in breeding populations depending on the results from a | 
study designed to evaluate geographic variation among our plantation selection popuia 
(Anonymous 1987). And. 2) The 600.+ second-cycle selections will have been tested, s0 .^- 
estimates of breeding values for each tree will be available. Had we started with untested sele ^^|| 
1 60 trees with poorly estimated breeding values would be a bigger risk and probably unactept* J" 



$oxh parental and offspring selections will be used in the third-cycle breeding program. The top 160 
clones based on breeding value estimates from current tests, and some procedure to limit coancestry 
be utilized. 

The number of sublines and the number of entries per subline is primarily dictated by the size 
0 f the breeding population and the number of unrelated clones needed to establish orchards. With 1 60 
selections, and subline sizes of 4 clones, 40 sublines would be available for establishment of orchards 
0 r other types of production populations. There are several benefits from making the sublines as small 
as possible that include: 1 ) Inbreeding within a line will increase more rapidly, thus increasing the total 
genetic variance from which to make gain (Falconer 1989); 2) The 40 sublines allow for exploitation 
of the increased genetic variance (1 +F) among lines, since selections from only the best lines would 
be used for orchard establishment each generation; 3) Alleles can be fixed faster in small populations 
(Kang 1991); 4) If selection criteria change in future generations, small sublines will be more 
responsive to selection (Namkoong et al. 1988); 5) Breeding and testing each subline will move more 
rapidly with small sublines. Each can be managed as a discreet unit; 6) Small sublines can be 
combined into larger sublines in the future. If alternative population management strategies are 
deemed appropriate for future breeding cycles, the populations can be readily restructured; and 7) A 
minimum effective population size of 40 can be maintained even after many generations of breeding. 

Some of the potential disadvantages of maintaining small sublines include: 1 ) Inbreeding may 
be too rapid so that alleles (both favorable and unfavorable) are fixed too rapidly and selection is not 
effective; 2) Inbreeding will reduce fecundity so that seed production for breeding and operational 
needs is reduced; 3) Recombination is restricted to within lines in the breeding population; and 4) 
Sublines will be eliminated due to genetic sampling, drift, and/or selection. 

A sublining system is not only effective for managing inbreeding, but it is also very effective 
for managing a series of small, homogeneous populations of trees ordered along a continuous 
environmental gradient. There are ho discrete "breeding regions" in the Cooperative, but a continuous 
series of populations. Since environmental factors which affect genotypic performance are continuous, 
different selection criteria can be utilized in a continuous fashion in managing the sublines. For 
example, the importance of fusiform rust resistance gradually changes both N-S and E-W. The 
economic weight for rust resistance can be changed in a selection index for sublines occurring along 
rust gradients). Sublining also simplifies breeding, since each cooperator will manage a given 

set of sublines. The logistical problems of exchanging selections, pollen, seeds, and seedlings for 
breeding, testing, and selection will be eliminated. 

The mating design for the sublining system will be a complementary design that combines 
pollen-mix mating to estimate GCA's and control-pollinated matings for within-family selection, van 
Buijtenen and Burdon (1990) and Burdon and van Buijtenen (1990) emphasize the value of pollen mix 
designs combined with control-crossed designs. A pollen mix of 20j^ average clones for each of the 
eight Test Areas (Anonymous 1984) in the Cooperative will be used. An adjustment factor will be 
calculated for each pollen mix based on standard check seedlots. 

Control-pollinated families will be used only for within-family selection. With 4 parents in each 
subline, 6 crosses among them will be made. Seedlings from the control-pollinations will be planted 
in non-replicated full-sib 36-tree family blocks, and the best phenotypes will be selected. The crosses 
from which selections will be made will be determined by the results of the pollen mix tests. Crosses 
with the highest mid-parent value will be identified, and the best individual tree(s) within those crosses 
will be selected. 
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ELITE POPULATIONS 



The extra Qain to be obtained from breeding only the most highly selected trees for specific 
geographic regions (e.g. Cotterill et al. 1 988) was the main attraction to including elite populations intc 
our hierarchy of breeding populations. There is a large increase in selection intensity by utilizing only 
the top clones in a recruitment population rather than all 1 60 clones. The average current recruitment 
population is about 600 clones (Anonymous 1982). In the mainline population, 160 of these clones 
will be utilized for a selection intensity of i » 1 .231 . In an elite population, the selection intensity for 
40 out of 600 is 1.937 or 57% higher than for the mainline population. The selection intensity 
benefits derived from reducing the population size from 600+. will occur only once, since populations 
are not likely to be reduced in size in future generations. Another major advantage of breeding elite 
populations is rapid generation turnover. It is much quicker to breed only 40 entries rather than the 
160 entries in the mainline, thus gain per year will be higher (Mahalovich 1990). 

There are numerous options available for management of elite populations. Rather than 
developing one elite population breeding strategy for all cooperators, elite populations will be 
customized for each cooperator. Individual cooperators or groups of cooperators will decide how their 
elite populations will be managed. Populations can be structured using local or non-local genotypes 
to utilize only general combining ability effects or both general and specific combining ability effects. 
Since one of the main benefits of elite populations is that breeding cycles can be completed more 
rapidly than in the mainline, any strategy adopted must have fast generation cycling as a feature. 

One option we are evaluating is the use of inbreeding is to rapidly increase homozygosity and 
the frequency of favorable alleles as well as to increase additive genetic variance. Another advantage 
of very rapid inbreeding (e.g. selfing) is that the among-family selection intensity for selfed families is 
greater than for bi-parental crosses. With selfing, the best parent based on GCA performance is mated 
with itself and can be selected without reducing its breeding value by crossing it with another parent 

While selfing has potential to give the most rapid gain, there are concerns. Getting an adequate 
number of selfed seeds (200 jj is difficult or impossible for many clones. Based on past selfing trials, 
selfing will only be feasible for about 50% of the clones. If an elite population size of 40 clones is 
utilized, seeds from about 20 clones can probably be obtained at a reasonable cost. Selfing may also 
increase homozygosity so rapidly that selection may not have a chance to operate. Slower inbreeding 
such as full-sib mating or half-sib mating could be more effective in the long run. For the second and 
future cycles in the elite population, there are numerous options available for advancing to the next 
generation. Selfing can continue and S 2 offspring from the selected S/s can be produced. Progeny 
from outcrossed trees could be selfed to create new S/s. For milder levels of inbreeding that might 
allow selection to operate more effectively, crosses between full-sibs (9 = .25) or crosses among half- 
sibs (6 = .1 25} can be made. There are also numerous options for crosses among all clones within the 
elite population (e.g. among S,'s or among any selections within or among sublines would be possible! 

Additionally, there is no reason to consider the elite populations as "closed". There will be 
overlap in breeding values between trees in the mainline and elite populations (i.e. heritabilities art 
relatively low and selection is far from perfect). Each cycle, the best genotypes from the mainline 
population in the same geographic area or from other elite and mainline populations can be infused into 
the elite population if the breeding values meet selection criteria. Infusion of new genotypes will a 5 
enhance the genetic variation in the elite populations. 

If future breeders desire to maintain a sublining system within the elite population, related tre** 
f*om within each diallel can be maintained within lines. The full-sib and half-sib crosses could be 
*va V o f maintaining the genetic integrity of a subline. In fact, the full-sib crosses described above 
t>* considered as 2-tree sublines from within the 4-tree diallels. The setfs are nothing more than o 
sublines if S, # s are made. Flexibility for future options is a key advantage of this elite popula 
Management scheme. 
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GENETIC DIVERSITY ARCHIVES 

Th„ ™«d to conserve forest genetic resources has been discussed in a recent report by the 
.^^^S^wTlWI)- We have made no attempt to elaborate on the economy 
National Research iss ' recognize the need to manage the extremely 

biological, or eth.cal needs J£££*ll populations that we have developed 

diverse and nth gene poor ^J^J^^^^fSl^ and forest scientists today and 

-rting immediately if we are to 

insuta tna availabilitv of «hi« raaourca in tha tutura. 

c, ,.rm nraaarvation of ganatic diversity In our lobloKy Pine population!, an additional 
*JZ£S£Z EpSTSoTthi mainline br«din, progr«nl wil b. th. iudidou. us. of Gan.,ic 

M, r^o Jne ban^ Tc^ncerted effort will be made to preserve these selections in the genet.c 

m f„tur« breedina cycles rather than lowering the selection intensity in the mainline populat.cn 

u u\ -~ TK-« fl trflfl^ will have uniaue and extreme values for different traits le.g. very 

^onXneX le archives is warranted, all trees will be available for th.s effort. 

Genetic resources are also available from other loblolly pine breeding programs in South and 
from around the world. The Western Gulf Forest Tree '^7™!^ 

^^^^^^^ ttttlZS* event that 
desirable alleles are not available in breeding or archived populates. 

w 6 have taken a very conservative approach to managing and preserving the genetic "sources 
of our SrrpSi- "amkoong "™ ""l^ 

necessary for managing genetic diversity The mutaple ^^^^^^SZ^ much 
conservation strategy. "As a part of breeding populat.cn management. *™™™*™™" s identical to 

as a low-cost insurance of diversity- 

SUMMARY 

The breeding and testing workloadsfor the Cooperatives 
smaller then the workloads in the current breed.ng cycle. We "^^^£££4 crosse s in 
including polycrbsses. will be 50 to 1 30 per cooperator. or about ^'^^X^xor. 0 r 
the 2nd cycl* program. The testing program w... consist or 7 ttT 5- » ' " °< tes « S mp ,e aTpossib.e 
about 1/10 to US the size in the 2nd-cycle. Every e «°* to J^ populations 
in the field have been made. While the management and coord.nat.on of the arnerem p w 
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(especially elite populations with potentially different breeding strategies) will be complex, the actual 
breeding, testing, and selection activities by cooperators will be much simpler than the diallel system 
we currently use. 

The breeding strategy was developed without having made extensive genetic gain calculations 
for various breeding strategies. A program of high-priority research is planned to evaluate population 
structure and effective population size issues. Alternative theories of population management will be 
evaluated by computer simulation, and the most promising of these will be tested in model breeding 
populations. Model populations will be an integral part of our program's research effort since the 
information is required to guide our future breeding. 

A final issue concerns the compatibility of the third-cycle breeding plan with new technologies 
in forestry and genetics. If, for example, vegetative propagation of mature loblolly pine becomes 
feasible, breeders may want to reconsider breeding for specific combining ability as well as for general 
combining ability in the mainline program. The small sublines lend themselves very well to a paired 
subfining strategy or a reciprocal recurrent selection strategy to advance SCA and GCA (McKeand et 
al. 1 986). Another option available to any cooperator is to advance SCA and GCA in elite populations. 

Vegetative propagation could also change and enhance the testing and selection program. If 
individuals within families can be clonally Replicated, the within family herrtabilrty 4yould be much 
higher, and gain from within family selection would be increased. Additionally, to exploit SCA and 
GCA more fully in production populations, additional crossing among clones and planting of these 
crosses in replicated field trials may be desirable. 

New developments in molecular genetics can also be incorporated into the Cooperative's 
breeding strategy. If molecular marker aided selection (MAS) for quantitative traits becomes feasible 
(we are currently collaborating with our molecular genetics group on a large experiment to test the 
techniques), the line breeding and particularly the inbreeding in the sublines will be a major benefit. 
Marker aided selection depends upon the maintenance of linkage disequilibria between the marker loo 
and the quantitative trait loci (Lande and Thompson 1990). Line breeding as opposed to random 
mating among individuals is very effective at maintaining linkage disequilibria for multiple generations, 
and enhancing the efficiency of MAS. 

For other biotechnologies such as recombinant DNA to be applied to forest tree populations, 
breeders must have a better understanding of the genetic control of specific traits. The elite 
populations lend themselves to very intensive genetic experiments where breeders can intensively 
study a few, highly selected individuals to conduct the necessary physiology and biochemistry 
experiments. 
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Burdon, R. D., and van Buijtenen, J. P. 1990. Expected efficiencies of mating designs for reselection of parents 
Can. J. For. Res. 20: 1664-1671. 

Expected efficiencies, in terms of genetic gain from reselection of parents (backwards selection), were compared 
for hierarchical mating, factorial, partial factorial, modified half diallel, and partial diallel crossing designs and 
polycrosses; this was done in parallel with a separate study of expected efficiencies for advanced-generation (forwards) 
selection, assuming a fixed-resources model of 100 parents and 10 000 offspring. The present study considered a single- 
trait case, with variable numbers of crosses per parent, varying heritability (h 2 ) levels (A 2 = 0.1, 0.2, and 0.5), and 
varying ratios of specific combining ability to general combining ability variance (0, 0.5, and 1). Compared with the 
case of forwards selection, the relative efficiencies of the different designs were generally similar. Two notable excep- 
tions were the comparative inefficiency of small, disconnected factorial sets for backwards selection and the generally 
high (but not always maximal) efficiency of polycrosses for this purpose. 



Burdon R. D., et van Buijtenen, J. P. 1990. Expected efficiencies of mating designs for reselection of parents. Can 
J. For. Res. 20 : 1664-1671. 

Les performances attendues, en terme de gain genetique obtenu de la reselection des parents (selection recurrente) 
sont comparees pour des modeles de croisement hierarchiques, factoriels, factoriels partiels, semi dialleles et dialieles 
partiels modifies, et polycross. Ces travaux ont ete realises en parallel* avec une autre etude sur la performance estimee 
dans la selection de generations avancees (nouvelle selection), en assumant une situation avec un materiel fixe de 
100 parents et 10 000 descendants. L'etude porte sur un seul caractere stable avec des nombres variables de croisements 
par parents, avec des niveaux d'heritabilite (h 2 ) variables (A 2 = 0,1, 0,2 et 0,5) et differents rapports entre Inaptitude 
specifique au croisement et la variance de l'aptitude generate a ce dernier (0, 0,5 et 1). Comparativement a la nouvelle 
selection, les performances relatives des differents modeles etaient generalement semblables. II y avait deux exceptions 
importantes : la performance comparativement faible avec des petits groupes factoriels non relies dans le cas de la 
selection recurrente et la performance generalement elevee des croisements polycross dans ce contexte. 

[Traduit par la revue] 



Introduction 

Expected efficiencies of alternative mating designs, for 
giving genetic gains in advanced-generation (or forwards) 
selection, was analyzed by van Buijtenen and Burdon (1990). 
It was shown that provided designs were balanced in repre- 
sentation of parents and provided the parents represented 
a random sample of potential selections, a variety of designs 
giving full-sib family information differed very little in gain 
efficiency for a given number of crosses per parent genotype. 

As was pointed out, mating designs may need to serve 
several purposes that might all have to be considered in 
choosing designs. Relative efficiencies for the various pur- 
poses may therefore need to be considered jointly for pur- 
poses of decision making. For minimizing allele loss in 
breeding populations, Kang and Namkoong (1979, 1980) 
concluded that balance (equal representation of all parents) 
was the paramount consideration. For estimating genetic 
parameters, Namkoong and Roberds (1974) made a com- 
prehensive comparison of designs, although this function 
could often become less important for the next generation 
or so (Namkoong 1970) as genetic parameters would often 
become reasonably well known. For testing or reselection 
of parents (backwards selection), work has also been done 
(van Buijtenen 1976; Lindgren 1977; Pepper and Namkoong 
1 978), but it remains of interest to compare designs in situa- 
tions parallel to those in which they are compared for 
advanced-generation selection. Accordingly, the present 
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paper compares mating designs for reselection of parents 
in the same situations as covered by the companion paper 
by van Buijtenen and Burdon (1990). 

Genera] method 

As described by van Buijtenen and Burdon (1990), selec- 
tive efficiency is considered in terms of genetic gain for given 
numbers of crosses per parent. For the present study, how- 
ever, mass selection is not a relevant alternative against 
which other options might be compared, so expected gains 
are calculated as absolute values per unit selection intensity. 

Expected gains are formulated according to the modified 
model of Osborne (1957), as described by van Buijtenen and 
Burdon (1990). Briefly, expected genetic gain (AG) is 
expressed as 

[1] AG = / yfV; 

where / is the intensity of selection (reselection) in standard 
deviations and V f is the variance of the selection index. 
In subsequent derivations, the following coefficients are 
used: 

d - number of trees per plot 
/ = number of females 

k = number of crosses in which a parent is involved 

(in diallel) 
/ = number of locations 
m - number of males 
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p - number of parents in diallel 
q - total number of parents 
r = number of replications 
t = number of sets 

u - number of females/number of males in factorial 
jsiow the index, /, is for the single-trait case of the general 
form 

(2] / = £ h\ Az 
z 

where Az in this case is the effect of a seed parent in the 
progeny test classification, or of any group or subgroup 
within which the parents are nested and hi is the herit ability 
of the effect in question. Because the selection in this case 
is not being done within the offspring of the parents in ques- 
tion, within-family genetic sampling effects contribute 
nothing to the numerators of the relevant heritabilities 
(cf. van Buijtenen and Burdon 1990). 
For V/ we have 

13] V, = £(^)VAz 

Z 

where s 2 Az = variance of effects within category z. In turn 



[4] s 2 Az = ^ 



df, 



N z df, + 1 



SS, 
N 



when z denotes sets of parents. 

df z = degrees of freedom of clones (or subclasses) for 
category z 

MS 2 and SS Z = mean square and sum of squares, 

respectively, for category z . 
N z = number of individuals per class or subclass for 
category z 

N = total number of individuals in the progeny test 
classification 

However, s 2 Az equals 4SS Z /N when general combining 
ability (GCA) information on individual parents is being 
used, since GCA information reflects only half the additive 
genetic value of a parent. 

The progeny test situation 

The following situations are assumed (using the arbitrary 
convention that rare-parent classes are designated male and 
common-parent classes are female, although males can be 
equal in number to females): 

(/) hierarchical or nested design, with m males each mated 

to a different set of / females 
(//) factorial designs: 

(a) complete factorials, with m males crossed with 
each of / females within each of / disconnected 
sets 

(b) partial factorials, with /females each involved 
in k crosses with members of a group of m males 
within t disconnected sets 

for balance, / = urn, u > 1, / and m being integers 

(hi) diallel designs 

{a) modified half diallels, without reciprocals or 
selfs with p parents, within each of t discon- 
nected sets 

(b) partial diallels, with each of p parents involved 
in k crosses (k < (p - 1)) within each of / 
disconnected sets 



(/v) polycrosses, with each of the q parents pollinated with 
a mix from a large number of pollen parents 
Note that there are d trees represented in each individual 
cross. Limiting cases are the following: 

/ = 1 in hierarchical cross, single-pair cross 
p = 2 in modified half diallels, single-pair cross 
k = 1 (therefore p = 2) in partial diallels, single-pair 
cross 

t = 1, single connected sets 

k = p - 1, ordinary modified half diallel as special 
case of partial diallel 

k = um, factorial sets complete, instead of partial 
Assumptions made for present purposes include the follow- 
ing: (/) regularity of designs; (//') equal numbers of individ- 
uals (d) per cross within each design; (///') full randomiza- 
tion of individual offspring; (/v) all parents belonging to 
a single random population of possible selections, hence, 

(a) male and female groups have a common between-parent 
genetic variance (i.e., a 2 = = a\, where and a f 
are the variances among male and female half-sib families, 
respectively, and a\ = GCA variance, which is assumed to 
equal one-quarter the additive genetic variance (a 2 *)) and 

(b) any genetic variance among disconnected sets will derive 
from purely random genetic sampling error, so the two 
between-sets variance component (<r^ ls ) equals zero. 

From these assumptions, follow the analyses of variance 
already presented by Burdon and van Buijtenen (1990, 
Table 1). These analyses entail three variance components 
not yet mentioned here: a 2 (within (full-sib) cross 
variance), o\ (specific combining ability (SCA) variance), 
and o 2 < (variance within half-sib crosses (e.g., 
polycrosses)). The derivation of expected gains for specific 
designs is contained in the Appendix as follows: hierarchical 
crosses; factorials; diallels; polycrosses. 

Results 

Expected gains from different mating designs are shown 
in relation to crosses per parent, for the various combina- 
tions of heritability and a\/a\ y in Fig. 1. The impact of 
substantial SCA relative to GCA is clearly evident in the 
numbers of crosses needed per parent to approach maxi- 
mum gains (Fig. 1, Table 1). Where a\ = 0, the partial 
mating designs showed maximal gains as soon as k > 1 ; in 
other words, in all cases except single-pair mating. 

Polycrosses were consistently efficient in the presence of 
SCA, but were appreciably suboptimal with zero SCA and 
low heritability. Indeed, the numbers of crosses per parent 
required for control-cross designs to match the efficiency 
of polycrosses tended to be quite low, unless o\fo % and h 
were comparatively high (Table 2). 

Small disconnected factorial sets tended to be appreciably 
less efficient than the corresponding diallel sets (Fig. 1, 
Tables 1 and 2). However, where factorial sets are larger 
(i.e., m and /larger), the factorials compared much more 
favorably; a single partial factorial (not shown in Fig. 1) 
consistently approached to within 0.5% of the expected gain 
from a single partial diallel. This general feature of the com- 
parisons between factorials and diallels is further illustrated 
in Table 3 when k, the number of parents per set, and num- 
ber of sets, are traded off to give a constant 50 crosses per 
set. Where k is larger, the set sizes are correspondingly 
reduced and the factorials show poorly. 
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Fig. 1. Comparisons of expected gains per unit / from different designs, under different combinations of A 2 , a\/a l v in relation to 

jotal crosses per parent. Designs considered are hierarchical ( ), disconnected modified half diallels ( ), single partial diallel ( ), 

and disconnected balanced factorials ( — ). Asterisks denote polycrosses. Sets are assumed to be fully cross-referenced. 



The case of the single-pair cross was always substantially 
suboptimal, particularly as <j\ increased. Even so, with the 
low value of h 1 (0.1) and zero SCA, the expected gain from 
single-pair crosses, under the conditions assumed, was 79% 
of that from polycrossing. 

Hierarchical crosses tended to be decidedly inefficient, 
unless / was large and o\/o\ was small (Figs. 1 and 2). 

The significance of cross-referencing sets, which amounts 
to being able to include a term for sets in each of the rele- 
vant equations, is also illustrated in Table 4. Where sets are 
small, the contribution of such cross-referencing is clearly 
important. It is also slightly more important with lower 
heritability and with the higher ratios of o\/o% A further 
illustration of the impact of set size is given in Table 5, for 
factorials in which total crosses per set are fixed at differ- 
ing numbers but with no cross-referencing between sets. 

Figure 3 illustrates the significance of imbalance between 
numbers of males and females for different cases of fac- 



torial crosses and hierarchical crosses. The reduction in 
expected gain was greatest with the intermediate level of 
imbalance, but imbalance increases set size in relation to 
crosses per parent and thus reduces the practical scope for 
cross-referencing sets. For hierarchical crossing, however, 
increasing imbalance (/large) improved efficiency (Fig. 2), 
but this design was generally inefficient compared with the 
others, particularly when SCA was present. 



Discussion 

Validity and robustness of results 

A study of this sort must depend largely on the assump- 
tions made. Most of the assumptions have already been dis- 
cussed by van Buijtenen and Burdon (1990). However, some 
are likely to be more critical for backwards selection than 
for forwards selection and are accordingly considered again 
here. 
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Table 1. Approximate genetic gains, as a 
percent of those realisable with 25 crosses per 
parent, with 1, 2, 3, and 4 crosses per parent 
for the cases where a, = a g 





Total crosses 


per 




parent (i 


.e. ( k/2) 


Mating design 


1 


2 


3 


4 


Partial diallel 


84 


91 


94 


96 


Partial factorial 


84 


91 


94 


96 


Disconnected diallels 


79 


89 


93 


95 


Disconnected factorials 


76 


87 


91 


93 



Notc: Results arc those cases where applicable sets are 
assumed to be cross-referenced with respect to genotypic 
value. 



Table 2. Numbers of crosses per parent (k/2) required for 
various types of balanced mating design to exceed the efficiency 
of polycrossing, for various values of h and a s /a g 



Type of design 





Partial 
diallel 


Disconnected 
diallels 


Partial 
factorial 


Disconnected 
balanced 
factorials 


0 


1 


1 


1 


2 


0.5 


1.5 


2 


1.5 


2.5 


1 


2.5 


2.5 


2.5 


4 


0 


1 


1 


1 


3 


0.5 


2.5 


3 


3 


5.5 


1 


5 


5.5 


5.5 


7.5 


0 


1 


1 


1 


7 


0.5 


7 


7.5 


8.5 


13 



Regular mating designs, without missing crosses, have 
been assumed throughout. Where each parent is crossed with 
only a few others (i.e., k is small), the GCA estimates for 
the parents involved in missing crosses will be severely com- 
promised, with potential effects on the overall efficiency of 
selection. The loss of information from missing crosses could 
be far more severe than in forwards selection, which leaves 
scope for within-family selection, except of course in the 
missing crosses themselves. The comparative robustness of 
small diallels and small factorial sets has apparently not been 
documented. However, small factorials sets (m and/small) 
seem inappropriate on grounds of inherent inefficiency 
while it seems unlikely that large diallel and large factorial 
sets would differ materially in robustness of genetic gains. 

Where k is small it should be easier to approach the full 
number of intended crosses. In this connection, adopting 
small k values should make it easier to circumvent in diallel 
difficulties where particular parents are unsatisfactory pro- 
ducers of pollen or female flowers. Factorial crosses 
automatically accommodate such flowering constraints. 

For the partial designs it is assumed that the precision of 
estimating half-sib family means is as implied by the 
expected mean squares given in Table 1 of van Buijtenen 
and Burdon (1990). Some noteworthy results hinge on this 
assumption, so cross-checking by Monte Carlo runs is sug- 
gested as a follow-up study. 

For polycrosses it is assumed that pollen parents are 
numerous and that they make equal contributions to the of f- 



Table 3. Comparative expected gains from disconnected 
partial diallels and disconnected balanced factorials in 
which numbers of parents per set and total crosses per 
parent are traded off to give 50 crosses per set 



Crosses per 
parent 
(k/2) 


Disconnected 
diallels 


Disconnected 
balanced 

f o^trt ri n 1 C 
lo\,lUI lala 


i CLl CIllJ 

per set 


0.5° 


27.52 


27.52 


100.0 


1 


38.18 


37.89 


50.0 


1.5 


39.20 


38.72 


33.3 


2 


39.64 


38.96 


25.0 


2.5 


39.82 


38.93 


20.0 


3 


39.87 


38.77 


16.6 


3.5 


39.84 


38.52 


14.3 


4 


39.75 




12.5 


4.5 


39.63 




11.1 



Note: Sets are assumed not to be cross-referenced with each other, and 
h 1 ~ 0.2, a\/o\ = 0.5. Noninteger values of crossing parameters have 
been accepted. 

'Single-pair crosses. 



spring. If numbers in pollen mixes are not large, then 
unequal contributions will clearly become much more 
important. 

An absence of maternal effects has been assumed. Even 
if transient maternal effects exist, the breeder should use 
methods for correcting for such effects (cf. Burdon and 
Sweet 1976). 

Also assumed is that prior selection has not truncated the 
genetic variation among parents. This means that the a g 
values in the numerators in eqs. Al-Al 1 (Appendix) equal 
one-quarter of the base-population additive genetic variance. 
Some departures from this assumption will occur and will, 
in particular, increase the relative importance of <r g . 
Adjustments for truncation, effectively revising eqs. 12-16 
of van Buijtenen and Burdon (1990), should be possible. 

The weightings given to male and female information in 
formulating the efficiencies of hierarchical crosses and 
unbalanced factorials (eqs. A1-A6, Appendix) is not so 
automatic as for forwards selection. The basis adopted, 
namely proportions of males and females, may be arguable, 
but other plausible alternatives seem unlikely to affect the 
general picture. However, imbalance would presumably 
increase sampling variation about the expectations. 

A common time scale is implied for the different mating 
designs. However, any differences between designs in time 
required for completion will affect the comparative efficien- 
cies in practical terms (Cotterill 1986). 

Comparisons among designs 

A consistent and predictable feature is that where a s is 
appreciable relative to a\ 9 high efficiency depends on cross- 
ing each parent with a number of others. Apart from the 
general limitations of hierarchical mating, the one respect 
in which type of pair-cross design seems important is the 
relatively low efficiency of small factorial sets. This reflects 
the fact that males and females are not cross-referenced 
against each other, except in terms of expected sampling 
deviations from the population mean. 

A notable, and somewhat counter-intuitive, finding was 
the high efficiency, in the absence of SCA, of the partial 
designs where k > 1 but is still small. As mentioned earlier, 
this may warrant semiempirical checking. 
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Fig. 2. Detailed plotting of expected gains from hierarchical crossing in relation to total crosses per parent and /, for different com- 
binations of h and o\/o v denotes a\ = 0, denotes a\ = (1/2)ct*, denotes a\ = o\. 



Table 4. Expected proportion of potential genetic gain for 

varying h 1 and a\/o\ 
(A) Comparisons among sets sacrificed in disconnected modified 
half diallels 



Table 5. Comparisons of expected gain from 
disconnected balanced factorials when numbers 
of parents and total crosses per parent (Jt/2) are 
traded off to give 50, 24, and 12 crosses per set 









No. of parents per set (p) 




h 2 


a\/a\ 


3 


4 


5 10 


20 


50 


100 


0.1 


0 

0.5 
1 


0.788 
0.757 
0.734 


0.853 
0.840 
0.830 


0.887 0.949 0.978 
0.880 0.947 0.977 
0.874 0.946 0.977 


0.994 
0.994 
0.994 




0.2 


0 

0.5 
1 


0.804 
0.768 
0.742 


0.861 
0.847 
0.836 


0.893 0.951 0.979 
0.885 0.950 0.978 
0.879 0.949 0.978 


0.995 
0.995 
0.995 




0.5 


0 

0.5. 


0.815 
0.773 


0.867 
0.851 


0.897 0.953 0.979 
0.888 0.951 0.979 


0.995 
0.995 




(B) Comparisons among sets sacrificed in disconnected balanced 
complete factorials 








No. of parents per set (2m 


= 2/) 




h 2 




4 


6 


10 


20 


50 


100 


0.1 


0 

0.5 
1 


0.811 
0.799 
0.790 


0.892 0.944 
0.887 0.942 
0.883 0.940 


0.977 
0.976 
0.976 


0.994 
0.994 
0.994 


1 
1 


0.2 


0 

0.5 
1 


0.816 
0.803 
0.793 


0.896 0.946 
0.890 0.944 
0.885 0.942 


0.978 
0.977 
0.977 


0.995 
0.995 
0.995 


1 
1 


0.5 


0 

0.5 


0.820 
0.805 


0.899 0.947 
0.892 0.945 


0.978 
0.978 


0.995 
0.995 


1 

1 



Interestingly, it appears that the polycross can be 
appreciably less efficient than various pair-cross designs in 
the absence of SCA effects, particularly when heritability 
is low. Several factors contribute to this situation. Polycross 



Total crosses per set 

Crosses per 

parent 50 24 12 



1 


37.89 


37.02 


35.30 


1.5 


38.72 


37.35 


34.58 


2 


38.96 


37.07 




2.5 


38.93 






3 


38.77 






3.5 


38.52 







Note: Sets are assumed not to be cross-referenced with 
each other and h 2 = 0.2, o]/a\ - 0.5. Noninteger values of 
crossing parameters have been accepted. 



families will have an inherently low signal to noise ratio with 
low A 2 , unless d is very large. At the same time, under 
polycrossing, each parent was effectively represented by only 
half as many offspring as under balanced pair-crossing 
schemes. Overall, the advantages of polycrossing seem less 
important than is widely supposed, unless SCA is 
substantial. 

The range of possible regular designs was not exhaustive, 
but it should suffice to provide the requisite insights. 

Other issues in choice of design 

Depending on the situation, the breeder may want to be 
able to use the same trial for both forwards and backwards 
selection. Whatever the purpose in mind, it seems worth 
seeking a design (or combination of designs) that is efficient 
for both purposes. Many of the background considerations 
have already been covered by van Buijtenen and Burdon 
(1990), so only a few additional remarks are made here. 

The main divergence in requirements for forwards and 
backwards selection arises in the comparative efficiencies 
of polycrosses and single-pair crosses, the former being 
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suitable for backwards selection and the latter for forwards 
selection. Polycrossing is always near optimal for backwards 
selection, but unless a\/o\ is high, is in itself substantially 
suboptimal for" forwards selection (van Buijtenen and 
Burdon 1990). It also has the major disadvantage, for 
forwards selection, of sacrificing full pedigree. 

Overall, the presence of SCA is more important in choos- 
ing mating designs for backwards selection than for forwards 
selection. One other feature is that small disconnected 
factorial sets should be avoided for backwards selection, but 
larger partial factorial sets should be manageable. 

Conflicts in requirements for forwards and backwards 
selection- may become appreciable in field designs. For 
backwards selection individual randomization (single-tree 
plots) is the most efficient. For forwards selection, however, 
multitree plots may enhance the within-family component 
of selection. 

Acknowledgements 
Part of the work was done while R.D. Burdon was at 
USDA Forest Service Southeastern Experiment Station, 
Department of Genetics, North Carolina State University, 
while holding a David Henry Scholarship administered by 
New Zealand Forest Products Ltd. Thanks are extended to 
C.B. Low for handling the computing. Thanks are also 
extended to Dr. P. P. Cotterill, Dr. G. R. Johnson, and 
Dr. J. N. King for reading the draft. 
Burdon, R.D., and Sweet, G.B. 1976. The problem of inter- 
preting inherent differences in tree growth shortly after plant- 
ing. In Tree physiology and yield improvement. Edited by 



M.G.R. Cannell and F.T. Last. Academic Press, New York, 
pp. 483-502. 

Cotterill, P.P. 1986. Genetic gains expected from alternative 
breeding strategies including simple low cost options. Silvae 
Genet. 34: 212-223. 

Kang, H., and Namkoong, G. 1979. Limits of artificial selection 
under balanced mating systems. Silvae Genet. 28: 53-60. 

1980. Limits of artificial selection under unbalanced mating 

systems. Theor. Appl. Genet. 58: 181-191. 

Lindgren, D. 1977. Genetic gain by progeny testing as a func- 
tion of cost. In Third World Consultation on Forest Tree 
Breeding, Canberra, Australia. Food and Agriculture Organi- 
zation and IUFRO. Pap. 6/9. 

Namkoong, G. 1970. Optimum allocation of selection intensity 
in two stages of truncation selection. Biometrics, 26: 465-476. 

Namkoong, G., and Roberds, J. 1974. Choosing mating designs 
to efficiently estimate genetic variance components for trees. 
Silvae Genet. 23: 43-54. 

Osborne, R. 1957. The use of sire and dam family averages in 
increasing the efficiency of selective breeding under a hierarchical 
mating system. Heredity, 11: 93-116. 

Pepper, W.D., and Namkoong, G. 1978. Comparing efficiency 
of balanced mating designs for progeny testing. Silvae Genet. 
27: 161-169. 

van Buijtenen, J. P. 1976. Mating designs. In Proceedings, 
IUFRO Joint Meeting on Advanced-Generation Breeding, 14-18 
June 1976, Bordeaux, France. Institut national de la recherche 
agronomique, Laboratoire d'amelioration des coniferes, 
Bordeaux, France, pp. 11-27. 

van Buijtenen, J. P., and Burdon, R.D. 1990. Expected efficien- 
cies of mating designs for advanced-generation selection. Can. 
J. For. Res. 20: 1648-1663. 



Page 1 of 2 



From: Satish.Kumar@ensisjv.com 
Sent: Thursday, July 20, 2006 1 :40 PM 
To: Lambeth, Clem 
Subject: Re: Polymix breeding 

Dear Clem, 

It is great pleasure to hear from you. Following your leading paper on PMX/WPA, we recently 
completed a study to further explore this concept. We submitted a paper to the journal Tree Genetics and 
Genomes, which has been accepted subject to few minor changes. I am in the middle of preparing a 
revised version, which I should be able to send to you in about one week time. 

In the mean time, I have attached here a reprint of our earlier paper on estimating relatedness using 
markers. I will be keen to pursue some similar projects in collaboration with you. 

Kind regards, 
Satish Kumar 

Scientist (Tree Improvement) 
Ensis-Genetics 

(Registered office: New Zealand Forest Research Institute Limited) 
Private Bag 3020 
Rotorua 
New Zealand 

(See attached file: molbreedpaper.pdf) 

□"Lambeth, Clem" <clem.lambeth@weyerhaeuser.com> 



"Lambeth, Clem" To<satish.kumar@ensisjv.com> 

<cIemJambeth@weyerhaeuser.com> cc 

SubjectPolymix breeding 

19/07/2006 10:31 

Dear Dr. Kumar: 

I have recently heard that you and some colleagues are doing work in an area dear to my 
heart, polymix breeding followed by paternity testing. Do you have any published or 
unpublished works that you could pass along (electronic versions would be appreciated)? I will 
be more than happy to honor any wish not to cite the work if it cannot yet be cited though 
citable works would be most helpful. I also have some new results in this area that may soon 
be written up for publication and any related work will be invaluable for the review section of 
the paper. 

Many thanks for your any information you can provide. 



Sincerely, 



Clem Lambeth 

Weyerhaeuser Company 

Tree Improvement Program Manager 

Phone: 5016248510 

FAX: 5016248505 



Tree Genetics & Genomes (2006) 2: 53-60 
DOI 10.1007/sll295-005-0033-x 



ORIGINAL PAPER 



Nicholas Wheeler • Peggy Payne • Valerie Hipkins • 
Robert Saich • Stephen Kenny • Gerald Tuskan . 

Polymix breeding with paternity analysis in Populus: a test 

for differential reproductive success (DRS) among pollen donors 



Received: 24 October 2005 / Accepted* 15 December 2005 / Published online: 2 February 2006 
© Springer-Verlag 2006 



Abstract Polymix breeding with paternity analysis (PMX/ 
WPA) has been proposed as an alternative to traditional 
full-sib breeding and testing schemes. To fully capture the 
benefits of PMX/WPA, differential reproductive success 
(DRS) of pollen parents used in the polymix must be 
modest. DRS was evaluated in an operational test of PMX/ 
WPA for a hybrid poplar breeding program. A 16-parent 
pollen polymix (Populus nigra L.) was used to pollinate 
seven clones of Populus deltoides (Bartr. ex. Marshall) 
under greenhouse breeding conditions. Progeny were 
grown out briefly and randomly sampled (357) prior to 
out-planting in field trials. Twenty-eight simple sequence 
repeat (SSR) loci were evaluated and 15 were selected for 
genetic characterization in small populations of three 
Populus spp (P nigra, P deltoides, and P balsamifera 
spp trichocarpa Torn & Gray). Seven loci were ultimately 
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selected for paternity analysis of progeny. The average 
exclusion probability of the seven loci in P nigra was 
0.604; combined, the theoretical exclusion probability was 
0.9999. However, only 95% of sampled progeny were 
unambiguously assigned a single paternal parent. Missing 
data likely accounted for most of the ambiguity. DRS was 
statistically significant though not prohibitive for practical 
utility of PMX/WPA as a breeding system. Of the 112 
potential crosses in this study, 92 were represented. Eight 
of the 16 pollen parents contributed 83% of the progeny. 
Good pollen vigor, as measured by germination percent, 
did not ensure paternal success, but poor vigor was as- 
sociated with lack of paternal success. PMX/WPA appears 
to be logistically and economically attractive for hybrid 
poplar breeding and testing. 



Introduction 

Polymix breeding with paternity analysis (PMX/WPA) of 
progeny has been proposed as an alternative to traditional 
full-sib breeding and testing schemes [4]. In short, PMX/ 
WPA requires the application of pollen mixes of many male 
parents in controlled cross events, followed by paternity 
analysis of progeny using molecular markers. The approach 
offers good estimates of breeding values of parents and pro- 
duces progeny populations with superior gain potential for 
forward selection. Other advantages of this approach include 
simpler and easier (thus, less expensive) breeding and testing 
programs, greater flexibility for inbreeding control, and con- 
venient connections with other molecular genetic applications 
such as clonal fingerprinting and marker-aided selection [4]. 
These features are particularly relevant for trees, organisms 
which offer a number of obstacles to breeders (i.e., large, 
long-lived, and late-maturing habits requiring costly breeding 
and testing programs and facilities). 

To fully capture the benefits of PMX/WPA, differential 
reproductive success (DRS) of parents used in pollen 
polymixes must not be high (most seed sired by only a few 
males). DRS in conifers (gymnosperms), though observed, 
is typically rather small (reviewed in Lambeth et al. 2001) 
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and insignificant relative to virtually all of the other factors 
controlling paternal success, such as amount and timing of 
pollen shed, location and distance of pollen donor relative 
to seed trees, and weather effects. Nakamura and Wheeler 
[8] concluded modest DRS in Douglas fir (Pseudotsuga 
menziesii Franco.) was likely due to genetic incompat- 
ibilities between male and female gametophytes (post- 
zygotic factors) but Nikkanen et al. [9] have demonstrated 
significant differences among spruce clones for pollen 
viability traits such as germination percentage and pollen 
tube growth rate. Such differences could arise from poor 
pollen vigor or genetic affects and could contribute to pre- 
zygotic competitive effects. Competition among male 
. gametophytes has been well-studied in annual and non- 
woody angiosperms [7, 10, 11], but controlled studies with 
flowering forest trees are wanting. 

The study described here results from the application of 
PMX/WPA for the creation of select hybrid poplar clones. 
Specifically, this study addresses the null hypothesis that 
all pollen parents in a polymix of pollens are represented 
equally in the progeny pool resulting from controlled inter- 
specific crosses with an array of seed parents (poplars are 
dioecious). Additionally, this study reports on the selection, 
optimization, and characterization of an array of micro- 
satellite simple sequence repeat (SSR) markers in three 
species of poplar. 



Materials and methods 

Plants, polymix construction, and controlled crosses 

In the fall of 2003, multiple scions from each of seven 
select maternal P deltoides Bartr. ex. Marshall (eastern 
cottonwood) clones originating from the upper mid- 
western United States, were grafted onto potted rootstock 
{P. deltoides var siouxland). In January of 2004, the grafts 
were moved into a greenhouse to encourage pistillate inflores- 
cence emergence and receptivity. Similarly, in mid-January, 
branches bearing staminate inflorescences were collected 
from 16 select paternal P. nigra L. (black poplar) clones and 
placed in tap water baths in a greenhouse (15-18°C) to force 
pollen shed. All P nigra clones originated from northern 
Italy; however, staminate cuttings were collected from an 
established clonal orchard located in eastern Washington. 
Seven of the 16 males had half-sib relatedness to one another 
(two pairs and a trio of half-sibs). Pollen was collected in 
glass beakers, funneled through 100 mesh screens to remove 
pollen catkins and poured into 20 ml glass vials (1/3 full) 
which were stoppered with cotton plugs. Vials were placed in 
desiccators containing silica gel and maintained at 4°C for the 
duration of the study (~1 week). 

Upon completion of pollen collection, pollen germina- 
tion trials were conducted by first rehydrating pollen ali- 
quots for 2 h on dishes floating on a water bath in a 
hydration vessel and then, placing them on 1% agar con- 
taining 10% sucrose and minerals [1]. Pollen germination 
was scored after 18 h and reported as the proportion of 
pollen grains exhibiting pollen tubes equal to or greater than 



the diameter of the pollen grain. Individual clonal values 
were placed in categories (<10, 1Q-15, 15-20, and 20-25% 
germination). Beginning February 23, pollen mixes were 
created daily by combining equal weights of desiccated 
pollen from all 16 male clones. Prior to brush-applied 
pollination, the polymix was rehydrated for 2 h as described 
above. Individual inflorescences were pollinated once per 
day for 3 days and retained within protective paper bags. A 
single controlled-cross polymix family was created for each 
of the seven maternal P. deltoides clones. Seed capsules set 
and began to burst within ca. 30 days after pollination. 
Seeds were collected by family, weighed, and sown in flat 
trays. The rising crop (ca. 30 days after sowing) was 
transplanted into Ropac #3 trays (96 cells) in April, 2004, 
grown for ~30 additional days, and subsequently out- 
planted to a field test site near Pasco, WA. The hybrid 
poplar cross created here (P. deltoides * P. nigra) is one of 
the most common and valuable inter-specific hybrid taxa in 
the world and is recognized by the designation Populus x 
canadensis or P euromericana. The hybrid is made with 
ease when P deltoides serves as the female parent, though 
the reciprocal cross (P nigra as female parent) virtually 
always fails [6]. 

Molecular marker development 

SSR markers were selected for this study because they are 
plentiful, co-dominant, and highly polymorphic in poplar. 
Though literally hundreds of SSR markers have been 
identified in poplar ([17] http://www. ornl.gov/sci/ipgc/ 
ssr_resource.htm), few have been fully characterized and 
optimized for high throughput detection. One of the goals 
of this study was to identify and genetically characterize an 
array of SSR markers that were useful in three poplar 
species (P deltoides, R nigra, and P balsamifera Torr. & 
Gray), the latter being a species native to the Pacific 
Northwest (black cottonwood) and common to hybrid 
poplar breeding programs. 

DNA extraction DNA extraction was carried out on leaf 
tissue (progeny seedlings) using the DNeasy-96 Frozen 
Leaf Tissue Protocol following manufacturers instructions 
with tissue homogenization achieved via the Mixer Mill 
300 (Qiagen, Valencia, CA) and from bud tissue (parents 
and population samples) using the DNeasy Plant Mini Kit 
(Qiagen, Valencia, CA) following manufacturers' instruc- 
tions. DNA quality and approximate quantity were assessed 
by comparing all samples against 50 ng of Lambda DNA 
standard on 0.8% agarose gels stained with ethidium 
bromide (EtBr) under UV light. 

SSR amplification and electrophoresis A total of 28 
candidate SSR primer pairs were evaluated (Table 1). SSR 
screening was achieved by amplifying approximately 1.5- 
2.5 ng of template DNA from P nigra clones in a 10-^1 
final volume including lx polymerase chain reaction 
(PCR) buffer, 2.0 mM MgCl 2 , 0.4 \M of each dNTP, 
0.4 jjM of the forward and reverse primers, and 1 U of 
HotStarTaq DNA Polymerase (Qiagen, Valencia, CA). Am- 
plifications were performed using a MJ Research PT-100 
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thermal controller with the following touchdown conditions: 
15 min at 95°C, 94°C for 30 s, 55°C for 30 s, 72°C for 1 min 
(3*), 94°C for 30 s, 52°C for 30 s, 72°C for 1 min (3x), 94°C 
for 30 s, 50°C for 30 s, 72°C for 1 min (29*) followed by a 
final extension of 72°C for 7 min. The amplification product 
was then diluted to a ratio of 1:50 (amplification:ddH 2 0) and 
' 1 \d of dilute amplification product was added to 10 jjI of Hi- 
Di Formamide containing 1.2% GeneScan-500 [ROX] size 
standard. Samples were subsequently denatured at 95°C for 
2 min and placed immediately on ice for 3 min before the 
sample plate was loaded on an ABI Prism 3,100 Genetic 
Analyzer for detection of SSR product ABI software 
packages, GeneScan Analysis Software and Genotyper 
Software v 3.7, were used to visualize and evaluate alleles 
at each locus. 



Molecular marker characterization 

Each of the 28 candidate SSR primer pairs were evaluated 
for use in determining paternity in each of the three poplar 
species using the criteria of genetic interpretability (i.e., 
does the primer pair produce a single locus that segregates 
in Mendelian fashion, and are the resulting alleles 
interpretable), heterozygosity level, exclusion probability, 
and expected null allele frequency. Primer pairs were 
initially screened using full-sib, controlled-cross family 
progeny of each species (five progeny each in P deltoides 
and P nigra and six progeny in P. trichocarpa). Based on 
inability to amplify and/or the presence of uninterpretable 
peaks, 13 of the 28 candidate primer pairs were discarded 
from further study. The remaining 15 primer pairs 
produced 15 loci that were determined to be potentially 
useful for paternity analysis. These primer pairs/loci were 
selected based on: 

1. 16 progeny from each of the three Populus species (48 
individuals; unrelated to individuals in the polymix 
study), 

2. Eight progeny from each of six controlled-cross 
families, i.e., two families from each of the three 
species (48 individuals), and 

3. 12 parent trees of the six families noted above. 

Parental genotypes were confirmed by simple exclusion 
analysis of the paternal contribution to each progeny in 
each species family thus, providing a basis for the 
evaluation of each locus in the areas of interpretability 
and presence of null alleles. Data from each species were 
scored using Genotyper Software v 3.7. CERVUS 2.0 [5] 
was used to calculate heterozygosity, parental exclusion 
probabilities, and estimated null allele frequencies. Also, 
the number, size, and frequency of alleles were calculated 
for each locus. 



Paternity assignment and reproductive success 

In August 2004, one leaf was sampled from each of -50 
randomly selected progeny from each of seven polymix 



derived families and from those pollen parents not pre- 
viously genotyped in the marker characterization exercise. 
DNA extraction, SSR amplification, and electrophoresis 
followed protocols noted previously. Paternity assignment 
for each progeny was determined using programming and 
formal language (PFL), a computer program for estimating 
pollen flow using paternity exclusion and SSR markers 
(http://www.fsl.orst.edu/pnwtirc/publications/pnwtirc _pubs_ 
date.htm). PFL allows for exclusion calculated on the basis 
of one or more locus genotype mismatches. All crosses 
were evaluated at both one and two genotype mismatch 
levels. Missing data were accommodated in the program. 
Where two or more paternal parents could have potentially 
sired a progeny;, partial paternity was assigned (e.g., two 
potential pollen parents=0.5 each). A test for deviation from 
expected reproductive success across all females was 
performed using chi-square analysis. 



Results and discussion 

Pollen germination and seed set 

Date of pollen shed varied from February 3 to 17, with 
individual clones shedding pollen for up to 7 days, in 
general, pollen germination rates in this study were low, 
ranging from less than 10% to around 25%. Pollen 
germination rates of 10-55% are considered sufficient for 
adequate seed set in most controlled crosses [13]; it is 
unknown whether variation in germination percent might 
contribute to pollen competitiveness in polycrosses of 
poplar. In the few controlled crosses performed with single- 
clone pollen lots exhibiting low pollen germination percent 
(<10%) in this study, seed set was nil, and inflorescence 
development was abnormal. Four .of the 16 pollen lots 
tested (25%) had germination scores of <10%. There was 
no relationship between date of pollen shed and germina- 
tion percent (r=0.04). 

Seed set was generally very good. Estimated seed set for 
polymix crosses exceeded 3,000 for six of the seven female 
P deltoides clones (680-4,860 seeds), while 20 seed 
weights ranged from 0.0083 to 0.0150 g. The maternal 
clone with a few seeds also had the lowest number of 
inflorescences to pollinate. 



Molecular marker characterization 

Thirteen of the selected SSR primer pairs produced 
ambiguous peaks or failed to yield an amplified product 
(Table 1). Fifteen SSR loci from the remaining 15 primer 
pairs were characterized for the three poplar species 
studied. Most loci were highly polymorphic in all three 
species: the average number of alleles per locus was 8.6, 
10.7, and 11.1 for P nigra, P deltoides, and P balsamifera, 
respectively. Across all species, 368 alleles were detected 
at the 15 select loci (mean=24.5 alleles per locus; 
range= 11-41; Table 1). In general, relatively few alleles 
were shared among species, particularly between either 



56 



1.8 

81 

c .2 

§•! 

& 

s 1 

O. o 

o o 

to O. 
g O 

8 i 

T3 -g 
§ *g 
11 

3 -a f 

US 

«£ 

00 o *= 

a ic * 
« 

.22 o*^ 
« 1| 

M C 
^3 

S-8 £? 
S-S.S2 

K fl U 

^ 3 g, 

o 5 OT 

■« s s 



3 -a 

■El* 

Cfl 

«;=> g 
£ S 

■« S J 



2 § 



S 

■c 

Q, 

e 



-a 



E 



o 


o 


o 




oo 


o 




SO 




00 


© 










© 


o 


o 


o 


m 




^- 


<N 




© 


© 




© 


© 


© 1 1 1 1 


o 

<N 


r- 






m 




(N 
co 


fN 


o 






<N 

m 


s 


<N 


<N | I 1 i 








© 




r- 


Os 


00 


<N 


VD 


OS 




<N 


r- 






en 


r- 




oo 


\o 


r-» 


Os 


in 


Os 




© 


r» 


oo 






© 


CM 




<n 


m 


r- 


m 


OS 


«n 




oo 


oo 


© 


© 



<<<<<<<< 

oooooooo 




^ H ^ PC ^ ffi £ 



XXX 



x: 



m £ b 



2 


© 

CN 


m 
m 


r» 


Os 

s 


© 








ir> 




u 


U 


GG 


U 


U 


o 


a 


O 


O 


a 


o 


S 


s 


2 


S 


S 










Pu 







- VI ^ V) 



s 



fS (N m \D OO OO 



S 2 

Cm 0. 



S2S 

pu eu pl, 



© 
U ^ 

11 

Oh !> 





m 
© 


© 


SO 
00 




© 


© 


tN 


tN 


CN 






i—l 








00 


a 


RN 








[vol 


o 


O 


o 


o 


o 


Oh 



© 

^ - V) VI \C OO 



s s 

Oh 0. 



SSS2 

O. Ch 'ft. 0. 



57 



North American species (P deltoides and P balsamifera) 
and the European native P nigra (5.2 and 6.8% alleles in 
common, respectively). P. deltoides and P trichocarpa 
shared 13.6% of their alleles. Only 10 alleles (2.7%) were 
shared by all three species. Null alleles were detected in 
controlled crosses for six loci but were not considered 
shared between species because they likely do not repre- 
sent the same mutational event. The mean allele frequency 
of the most common allele, across all 15 loci, was 0.308, 
0.342, and 0.444 for P trichocarpa, P deltoides and P 
nigra, respectively. Allele sizes and frequencies of all 15 
loci evaluated are provided elsewhere (Appendix 1: Gene 
frequencies: http://www.ornl.gov/sci/ipgc/poplar polymix 
genefreg Appendixl.xls) Seven of the 15 loci characterized 
were ultimately deemed to be 'informative'' for all three 
species (Table 2) and were used for paternity testing in this 
study. For the seven selected loci, average observed 
heterozygosity was 0.681, 0.700, and 0.623 for P 
trichocarpa, P deltoides, and P nigra, respectively. 
Similarly, the average exclusion probability for any one 
locus, given that the maternal genotype was known, was 
0.754, 0.704, and 0.604. Although the Italian sample of P 
nigra was the least genetically diverse of the three species 
evaluated, it theoretically possessed sufficient diversity to 



distinguish among all paternal parents. In the final analysis 
of all progeny and parents, the paternal exclusion proba- 
bility for all seven loci combined was 0.9999. Measures of 
diversity (number of alleles, heterozygosity) observed for 
this breeding population of P nigra were slightly lower 
than those detected in a previous study characterizing a 
range- wide collection of 23 individuals using nine SSR loci 
[16]. 



Genotyping 

Multi-locus genotypes were obtained for 357 hybrid poplar 
progeny. Over half of the progeny (58.3%) failed to 
amplify a paternal allele at one of seven loci; however, 
fewer than 5% were missing paternal genotypes at two or 
more loci. Null alleles accounted for a very small 
percentage of the missing data, while simple failure to 
amplify some primers in specific crosses appeared to be the 
cause in most cases. The failure to amplify a product was 
most acute on the maternal side (P deltoides). For two 
markers (PMGC 2235 and PMGC 2675), five of the seven 
maternal parents displayed distorted allelic segregation due 
to a failure to amplify a product. For the remaining five 



Table 2 Properties of seven SSR loci selected for paternity studies in P deltoides (D) by P nigra (N) crosses (P nigra is paternal parent) 



Locus 


Species 


A ' 


k 


N 




#(E) 


Excl (1) 


Excl (2) 


Null freq 


PMGC 649 


T 


94-158 


19 


19 


0.895 


0.957 


0.759 


0.863 


0.019 




D 


100-152 


17 


17 


0.765 


0.959 


0.753 


0.859 


0.098 




N 


92-108 


7 


18 


0.722 


0.754 


0.336 


0.513 


0.013 


PMGC 14 


T 


186-210 


7 


19 


0.632 


0.804 


0.406 


0.586 


0.118 




D 


188-191 


2 


19 


0.579 


0.508 


0.122 


0.186 


0.079 




N 


197-220 


7 


20 


0.600 


0.771 


0.350 


0.526 


0.116 


PMGC 433 


T 


198-230 


14 


18 


0.722 


0.927 


0.667 


0.801 


0.110 




D 


193-221 


15 


19 


0.684 


0.909 


0.630 


0.774 


0.118 




N 


196-228 


8 


18 


0.556 


0.514 


0.147 


0.321 


0.069 


PMGC 2221 


T 


67-138 


21 


18 


0.667 


0.959 


0.759 


0.863 


0.167 




D 


73-287 


19 


19 


0.842 


0.950 


0.739 


0.849 


0.048 




N 


81-136 


15 


19 


0.526 


0.937 


0.700 


0.824 


0.265 


PMGC 2235 


T 


111-152 


16 


18 


0.833 


0.941 


0.708 


0.829 


0.046 




D 


114-161 


15 


18 


0.611 


0.927 


0.670 


0.803 


0.190 




N 


111-132 


9 


19 


0.895 


0.898 


0.591 


0.745 


0.011 


PMGC 2675 


T 


165-232 


14 


17 


0.353 


0.922 


0.656 


0.793 


0.441 




D 


161-220 


18 


19 


0.632 


0.923 


0.669 


0.801 


0.180 




N 


152-208 


10 


17 


0.412 


0.863 


0.521 


-0.688 


0.347 


PMGC 2885 


T 


291-316 


7 


18 


0.667 


0.778 


0.364 


0.543 


0.067 




D 


297-316 


10 


19 


0.789 


0.839 


0.483 


0.657 


0.013 




N 


288-322 


8 


20 


0.650 


0.817 


0.430 


0.608 


0.100 


AVERAGE OVER ALL SEVEN LOCI 


T 




14.0 


18.1 


0.681 


0.898 


0.617 


0.754 


0.138 




D 




13.8 


18.6 


0.700 


0.859 


0.613 


0.704 


0.104 




N 




9.3 


18.7 


0.623 


0.793 


0.439 


0.604 


0.132 



Values are given for P trichocarpa as well. A Range of allele size in bp, k number of alleles, N number of individuals typed, Ht 0 \ observed 
heterozygosity, expected heterozygosity, Excl(l) average exclusion probability without information on one parent, Excl(2). average 
exclusion probability given information on known (female) parent, Null freq estimated null allele frequency, Use: 1 indicates likely useful 
locus, 2 indicates possibly useful locus 

♦Indicates issues with genetic interpretability 
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loci, only three of 35 cases displayed segregation distortion 
of maternal alleles 0^0.05, X >3.84, 1 df). Marker PMGC 
2235 is known to occur on chromosome IV which is 
comprised almost entirely of maternally biased segregation 
ratios [17]. Interestingly, two of the 16 paternal parents 
(numbers 12 and 16) possessed triploid genotypes and a 
third parent (number 9), a tetraploid genotype, for marker 
PMGC 2675. Whether this was causal in the aforemen- 
tioned amplification problem is unknown. To properly 
account for paternal parents with three or four, alleles in the 
paternal assignment program, additional dummy genotypes 
were created for these three parents. It should be noted that 
polyploidy in poplars is not uncommon [3]. 



Paternal reproductive success 

The PFL program assigned unambiguous paternity to 86% 
(309) of the progeny based on a single mismatch (one 
mismatch excludes a parent) and 94.7% based on two 
mismatches. Of the remaining 19 progeny, 15 were 
assigned two potential paternal parents, three were 
assigned three paternal parents, and a single individual 
could not be matched. For numerical summary purposes, 
multiple paternal assignments were treated proportionately. 
That is, each of two potential parents received credit for 0.5 
progeny, and each of three potential parents received credit 
for 0.33 progeny. The theoretical expectation of exclusion 
probability exceeding 99% was above the observed value 
of 94.7%. Several factors likely contributed to this 
shortfall: genotyping error, including missing data, prog- 
eny that were undifferentiated genetically (half-sibs), and/ 
or pollen contamination. Indeed, five of the 15 cases 
sporting two potential pollen parents assigned the half-sib 



paternal parents number 7 and number 8. Contamination, 
though unlikely, in the controlled-cross conditions of this 
experiment, was still possible and may have accounted for 
the lone progeny that had no parental assignment. 

All 16 paternal parents in the polymix were represented 
in the progeny pool, although representation was highly 
variable (Table 3). Clearly, the null hypothesis that all 
parents would be represented equally was rejected 
(pO.OOl, x 2= 219.1, 15 df). Over all seven females, the 
expectation is that each paternal parent would be 
represented by 22.3 progeny (3.2 per individual female). 
Actual paternal contribution varied from 2.5 to 62.0 (1.0 to 
17.0%) across all female parents and ranged from 0-12% 
for individual female parents. The top eight paternal parents 
sired 83.3% of the progeny. In general, the top eight parents 
were well represented in each cross, though there were 
instances in which these males failed to sire a progeny for a 
given female (i.e., Male number 1, Cross H4). The reverse 
was also true, i.e., a number of parents that faired poorly 
overall did quite well with individual females. The expected 
sample size for each male by female combination was too 
low to consider performing interaction chi-square tests. 

It is noteworthy that paternal success appeared to be 
unrelated to pollen vigor, as measured by proportion of 
pollen grains germinating in a short term test (r 2 =0.057). 
Indeed, the least represented pollen parent (number 8) had 
the highest pollen germination score (>25%). Still, repro- 
ductive success was clearly low for all parents with 
germination scores of one. Given the observation that a few 
single parent controlled crosses with pollen parents that 
had germination scores of one resulted in no filled seed, it 
is surprising to see any representation in the polymix 
progeny for these parents. One explanation could be the 
mentor pollen effect [14]. The mentor pollen technique has 



Table 3 Number of progeny assigned to each of 16 paternal parents in a polymix for each of seven crosses (H2-H9), the total number and 
proportion of progeny assigned each parent, and a measure of pollen viability based on germination percentage (<10%=1, 10-15%-2, 15- 



20%=3, and 20-25%=4) 



Female (P deltoides) 



Male (P. nigra) 


H2 


H3 


H4 


H5 


H6 


H8 


H9 


Total 


Proportion of total 


Pollen viability 


1 


2.0 


5.5 


0.0 


8.0 


8.0 


4.5 


4.5 


32.5 


0.09 


2 


2 


1.0 


0.0 


0.0 


0.0 


0.3 


0.0 


2.0 


3.3 


0.01 


1 


3 


9.0 


7.5 


6.0 


5.0 


4.0 


5.0 


4.0 


40.5 


0.11 


3 


4 


2.3 


4.5 


4.0 


4.0 


7.8 


10.0 


3.5 


36.1 


0.10 


4 


5 


4.0 


3.0 


8.0 


4.0 


1.3 


6.0 


10.0 


36.3 


0.10 


4 


6 


12.0 


7.0 


11.0 


7.0 


8.5 


11.0 


5.5 


62.0 


0.17 


2 


7 


0.8 


1.0 


1.5 


1.5 


0.0 


0.0 


3.0 


7.8 


0.02 


1 


8 


0.5 


0.0 


0.5 


0.5 


0.0 


0.0 


1.0 


2.5 


0.01 


4 


9 


5.3 


1.0 


2.0 


6.5 


1.0 


0.5 


2.5 


18.8 


0.05 


4 


10 


7.0 


11.0 


5.0 


5.0 


10.0 


3.0 


6.0 


47.0 


0.13 


2 


11 


0.0 


2.5 


0.0 


0.5 


1.5 


2.0 


4.0 


10.5 


0.03 


1 


12 


2.0 


2.0 


0.0 


5.0 


2.0 


1.0 


2.0 


14.0 


0.04 


3 


13 


1.0 


2.0 


2.0 


1.0 


1.3 


1.0 


0.0 


8.3 


0.02 


4 


14 


2.0 


2.0 


7.0 


0.0 


3.3 


7.0 


3.0 


24.3 


0.07 


2 


15 


1.0 


0.0 


1.0 


1.0 


2.8 


0.0 


1.0 


6.8 


0.02 


1 


16 


0.0 


0.0 


4.0 


0.0 


0.0 


1.3 


1.0 


6.3 


0.02 


4 


Totals 


50.0 


49.0 


52.0 


49.0 


52.0 


52.0 


53.0 


357.0 


1.00 
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been used to make inter-specific crosses that otherwise 
cannot be made. The desired pollen is mixed with 
irradiated (killed) compatible pollen before pollination 
and fertilization is facilitated. The combination of pollens 
in this case may have actually aided pollen lots of low 
viability. 

Finally, all three paternal parents that possessed poly- 
ploid genotypes for marker PMGC 2675 were under- 
represented in the progeny pool (mean of 13 progeny per 
parent, overall) even though they all possessed high pollen 
germination scores. 

Though the experimental design and statistical rigor of 
the experiment are insufficient to differentiate among all 
the possible reasons for unequal paternal success, the 
results do guide speculation. The lack of association 
between vigor and paternal success implies pre-zygotic 
competition based on number of pollen grains competing is 
not a critical factor although issues such as pollen tube 
growth rate may still be relevant (assuming rate of growth 
is unrelated to germination success). Post-zygotic compe- 
tition in the form of genetic complementation between 
male and female parents could explain much of the pattern 
of paternal success observed. This may.be particularly 
relevant in this case where inter-specific crosses were 
made. The correspondence between DRS and paternal half- 
sibs among the pollen parents appears to support this 
hypothesis. For instance paternal half-sib pairs number 7 
and number 8, and number 2 and number 15 each exhibited 
very poor reproductive success for both individuals (mean 
of 1 .5% in each case), while each of the three pollen parents 
in the half-sib that included number 5, number 12, and 
number 14 contributed significantly (mean of 7.0% per 
parent) to the progeny pool. 

Applicability of PMX/WPA for hybrid 
poplar breeding 

PMX/WPA was. conceived for recurrent selection breeding 
programs in which pedigree control is vital to reduce or 
eliminate inbreeding in advanced generations while simul- 
taneously improving gain and decreasing breeding costs. 
For inter-specific F! hybrid breeding programs, such as the 
one described here, pedigree control for inbreeding is not 
directly relevant, though knowing the parentage of 
individual progeny is still valuable. Managing a pedigree 
against inbreeding would be relevant if reciprocal recurrent 
selection programs exist for each parental species to 
identify the best parents to enter the hybrid program. More 
importantly, PMX/WPA offers a method for increasing the 
number of parental combinations created, ostensibly 
improving the chances of finding the best combination at 
a reasonable cost. This is particularly true if the chosen 
alternative approach is a disconnected factorial design. A 
potential tradeoff, relative to traditional factorial designs, is 
that fewer progeny per cross are created with PMX/WPA, 
assuming that progeny test size remains the same for the 
two designs. It is necessary to understand the genetic 



variance structure of "among parents" vs "within parents" 
to determine if this trade off is significant. While it is too 
early to determine if gain was improved in this study, it is 
possible to evaluate the costs of deploying PMX/WPA 
relative to the traditional factorial mating design used in the 
hybrid poplar breeding in this program. For the combination 
of activities that included grafting of breeding materials, 
pollen handling, making controlled crosses, and seed sowing, 
the polymix approach cost was -30% of the traditional fac- 
torial approach ($9,716.00 vs $3 1 ,556.00). The added costs of 
genotyping all parents and a few selected progeny would be 
nominal (<$500.00), making this approach economically 
quite attractive in an operational setting. 



Summary 

Significant logistical and economic incentives exist for 
adopting PMX/WPA for the creation of hybrid poplars in 
the program described here. Though not particularly 
surprising, the skewed paternal contribution to the progeny 
pool is noteworthy, and suggests follow-up studies to 
determine if pollen competition induced by polymix 
breeding prevents otherwise good crosses from occurring 
while simultaneously encouraging fertilization from pollen 
sources with low vigor. Of the 1 12 potential crosses in this 
polymix study, 92 were represented in the progeny pool. 
This contrasts with the 49 crosses that would have occurred 
under the factorial study. Even the least successful pollen 
parent could be represented by -70 progeny in the 
operational field tests (0.7% success rate in sample * 
9,800 progeny in test). Still, it would be desirable to have 
the remaining 20 potential crosses represented if the 
highest possible predicted gains are to be achieved. 
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Abstract The conventional way to drive modifications in 
old forest tree seed orchards is to establish progeny trials 
involving each parent tree and then evaluate its contri^ 
bution to the performance of the progeny by estimating its 
general and specific combining ability (GCA and SCA). 
In this work, we successfully applied an alternative parent 
selection tactic based on paternity testing of superior 
offspring derived from a hybrid seed orchard established 
with a single Eucalyptus grandis seed parents and six E. 
urophylla pollen parents. A battery of 14 microsatellite 
markers was used to carry out parentage tests of 256 
progeny individuals including two independent samples 
of selected trees and one control unselected sample, all 
derived from 6-year-old forest stands in eastern Brazil. 
Paternity determination was carried out for all progeny 
individuals by a sequential paternity exclusion procedure. 
Exclusion was declared only when the obligatory paternal 
allele in the progeny tree was not present in the alleged 
parent tree for at least four independent markers to avoid 
false exclusions due to mutation or null alleles. After 
maternity checks to identify seed mixtures and selfed 
individuals, the paternity tests revealed that approximate- 
ly 29% of the offspring was sired by pollen parents 
outside the orchard. No selfed progeny were found in the 
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selected samples. Three pollen parents were found to have 
sired essentially all of the offspring in the samples of 
selected and non-selected progeny individuals. One of 
these three parents sired significantly more selected prog- 
eny than unselected ones (P<0.0002 in a Fisher exact 
test). Based on these results, low-reproductive-successful 
parents were culled from the orchard, and management 
procedures were adopted to minimize external pollen 
contamination. A significant difference (P<0.01) in mean 
annual increment was observed between forest stands 
produced with seed from the orchard before and after 
selection of parents and revitalization of the orchard. An 
average realized gain of 24.3% in volume growth was 
obtained from the selection of parents as measured in 
forest stands at age 2-4 years. The marker-assisted tree- 
breeding tactic presented herein efficiently identified top 
parents in a seed orchard and resulted in an improved seed 
variety. It should be applicable for rapidly improving the 
output quality of seed orchards, especially when an 
emergency demand for improved seed is faced by the 
breeder. 



Introduction 

Eucalypts are the most widely planted hardwood trees in 
the world. Estimates made in 2000 (FAO 2000) indicated 
that the eucalypt plantation area was globally greater than 
17.8 million hectares, with India being the largest planter 
with 8.0 million hectares — mostly in extensive low- 
productivity plantations — followed by Brazil with 3.0 
million hectares of mostly intensively cultivated clonal 
plantations of industrial forests reaching average produc- 
tivities of 45-60 m 3 /ha per year (Mora and Garcia 2000). 
Elite hybrid clones consisting of Eucalyptus grandis and 
E. urophylla are extensively used by the cellulose and 
paper industry because of its wood quality, rapid growth 
and high volumetric yield (Bertolucci et al. 1995). 

In Brazil, eucalypt genetic improvement programs 
were initiated in the early 1970s. The main purposes of 
the breeding programs were to match the most adapted 
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species to the existing environmental conditions and to 
provide sufficient plant stock of reasonable genetic 
quality (Brune and Zobel 1981). Several E. grandis, E. 
urophylla and E. saligna seed orchards were established 
on the basis of initial results of available provenance/ 
progeny trials. Furthermore, due to the early observation 
of hybrid superiority for growth coupled to resistance to 
eucalypt canker caused by Cryphonectria cubensis (Bru- 
ner) Hodges, hybrid seed orchards of E, grandis x E. 
urophylla, were also established, thereby generating the 
so-called Urograndis seed variety (Campinhos 1980; 
Vigneron 1991; Eldridge et al. 1993). These open- 
pollinated seed orchards often involved the deployment 
of one or a few self-incompatible or male-sterile E. 
grandis trees as female parents and a number of selected 
E. urophylla as male parents to ensure enough pollen 
pressure (Bertolucci et al. 1995). 

Clonal propagation of elite trees started in the 1980s 
followed the initial breeding/testing/selection efforts. This 
new technology, involving the, establishment of true 
clonal stands through vegetative propagation via rooted 
cuttings, was developed and applied on a wide scale, 
thereby providing significant productivity and uniformity 
gains (Campinhos 1980; Brandao et al. 1984) and con- 
sequently reducing the demand for seedlings in the 
establishment of new forests. However, recent increases 
in the industrial demand for wood have forced forest 
companies to rapidly expand planting areas and establish 
partnership programs with farmers. The only way to ful- 
fill this immediate challenge has been to use seeds from 
first-generation breeding orchards to complement the 
production of rooted cuttings. The development of im- 
proved seed varieties has thus become necessary again in 
the context of current breeding programs, both by 
establishing new orchards with trees of high general and 
specific combining ability (GCA and SCA, respectively) 
and by restructuring old ones. 

The conventional way to drive modifications in old 
seed orchards is to establish progeny trials to evaluate 
the contribution of the parent to the performance of its 
progeny by SCA estimation. This approach requires a 
long time to achieve its purpose (Namkoong et al. 1988). 
An alternative way to estimate the superiority of the 
parent trees of a seed orchard could be attained by a 
marker-assisted breeding tactic: carrying out DNA-based 
paternity tests of superior progeny trees in forest stands 
established with open-pollinated seeds from the orchard. 
Parents that display a low frequency of superior offspring 
could be culled from the orchard, thus practicing a 
backwards selection that would result in an improved seed 
variety involving exclusively trees of higher GCA and 
SCA (Ribeiro et al. 1998; Grattapaglia 2000). 

A similar approach, termed PMXAVPA (polymix 
breeding with parental analysis) was evaluated as an al- 
ternative solution to full-sib crosses in pine breeding 
programs (Lambeth et al. 2001). As pointed out by those 
authors, although PMX breeding is easy, provides good 
estimates of breeding values and allows the testing of a 
larger number of parental combinations, pedigree control 



is lost. Instead of using a single pollen for each cross, the 
concept then proposed by Lambeth et al. (2001) was to 
use PMX breeding involving many male parents followed 
by paternity analysis of progeny with microsatellites, thus 
allowing full pedigree control. 

Similar to the situation in humans (Hammond et al. 

1994) and domestic animals (Glowatzki-Mullis et al. 

1995) , the high degree of multi-all elism and the clear and 
simple codominant Mendelian inheritance of microsatel- 
lites provide an extremely powerful system for the unique 
identification of Eucalyptus individuals in parentage 
testing, particularly when individuals are expected to be 
related (Brondani et al. 1998). In the study reported here, 
we applied a marker-assisted breeding tactic -to an old E. 
grandis x E. urophylla seed orchard with the specific 
objectives of verifying the differential reproductive suc- 
cess of six male parents and potentially identifying those 
with higher SCA with the maternal tree to develop an 
improved seed variety. 



Materials and methods 

Seed orchard 

The seed orchard studied was established in 1982 in Aracruz 
(latitude 19°49", longitude 40°16") in Espirito Santo State, Brazil 
by Aracruz Celulose S.A. Six Eucalyptus urophylla plus trees 
identified as being canker-free were selected from provenance/ 
progeny trials of open-pollinated, half-sib families originally 
collected on the Bessi-Lau, Flores and Timor Islands in Indonesia. 
These six trees were clonally propagated by grafting and consti- 
tuted the male parents in the seed orchard. A single E. grandis plus 
tree selected from a progeny trial of open-pollinated, half-sib 
families originally collected in Atherton (Australia) was used as the 
female parent. This particular female parent was selected due to its 
good growth performance and putative self-incompatibility that 
was deduced from its very low seed set in controlled self- 
pollinations (Bertolucci et al. 1995). The seven parent trees were 
clonally replicated and planted in a 6x6 spacing in 267 hexagonal 
plots totaling 7.5 ha, with each of the £". grandis maternal tree 
surrounded by the six pollinator trees. The orchard was surrounded 
by a 300-m wide belt of native tropical forest to provide isolation 
from any pollen coming from nearby eucalypt plantations. Seeds 
were only collected from the E: grandis maternal tree. 

Selection of progeny individuals 

Sampling of progeny individuals was carried out in 6-year old 
stands established with seeds derived from the seed orchard studied. 
Two independent samples of 72 selected trees each were taken from 
two different forested areas — SSA selected sample A and SSB 
selected sample B. Individual trees with a circumference at breast 
height (CBH) above one standard deviation from the mean (mean = 
55±15 cm) were selected, reaching a selection intensity of ap- 
proximately 1:200. A control sample of 72 random non-selected 
(RNS) trees was also taken to allow a comparative analysis with the 
selected tree samples. To verify potential trees derived from selfing 
of the maternal tree, we also investigated a sample of ten trees with 
stunted growth and a sample of 30 trees that displayed symptoms of 
chlorotic leaves. Total genomic DNA was extracted from adult 
leaf tissue following the protocol described by Grattapaglia and 
Sederoff (1994). 
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Microsatellite genotyping 

Forty-seven microsatellite markers developed by Brondani et al. 
(1998, 2002) were screened for polymorphisms among the seven 
parent trees in order to identify a battery of markers with higher 
information content for the proposed study. Microsatellite marker 
amplification and detection were performed as described earlier 
(Brondani et al. 1998). The amplified products were separated on 
4% denaturing polyacrylamide gels, stained with silver nitrate 
(Bassam et al. 1991) and sized by comparison to a 10-bp DNA 
ladder standard (GibcoBRL, Gaithersburg, Md.) on a computer 
screen. Allele sizes were estimated using the software seqaid ii 
(Rhoads and Roufa 1990), taking into consideration the expected 
allelic series in basepairs for the locus. Multiplex loading in the 
same gel was carried out for up to three microsatellite loci 
simultaneously. 



Paternity testing 

Gels were scanned, and scoring was carried out manually on a 
computer screen. Multi-locus genotypes were determined for all 
seven parent trees and the five samples of progeny individuals 
(SSA, SSB, RNS, stunted and chlorotic). Maternity of the E. gran- 
dis tree was first checked, and seed contaminants identified. 
Paternity determination was carried out for all progeny individuals 
by a sequential paternity exclusion procedure. Exclusion was 
declared only when the obligatory paternal allele in the progeny 
tree was not present in the alleged parent tree for at least four 
independent markers to avoid false exclusions due to mutation in 
the microsatellite marker. Paternity was ultimately declared when 
the alleged parent tree shared the obligatory paternal allele at all of 
the microsatellite markers tested. Because the number of potential 
pollen parents was limited, paternity declaration was deterministic 
and did not involve the calculation of a paternity index using 
likelihood-ratio methods. Progeny that were not assigned to a 
specific pollen parent were classified either as selfs or as having 
been derived from outside pollen contamination. 



Paternity in selected versus non-selected trees 

For each pollen parent tree that effectively sired progeny individ- 
uals, a two-tailed Fisher exact probability test (Fisher 1934) was 
used to test the null hypothesis that the proportion of sired progeny 
individuals is not significantly different between the selected and 
non selected offspring samples. 



Improvement of productivity with parental selection 

The male parents displaying the lowest reproductive success were 
eliminated from the orchard in the aim of producing a new 
improved seed variety. Realized gain in mean annual increment 
(MAI) (m 3 /ha per year) was evaluated by means of a comparative 
inventory between stands established with seedlings derived from 
the original orchard (78 stands) and seedlings derived from the 
improved orchard following the culling of the worst parents (18 
stands). These forest stands were 2-4 years old and had an average 
size of 20 ha with 1,111 trees, per hectare; they were located at 
uniform sites in Bahia State. Student's /-probability test was used to 
test for significant differences in MAI between the old and the new 
variety. 



Results 

Out of the 47 microsatellite markers screened, a battery of 
14 markers was ultimately selected and used throughout 
the study (Table 1). On the basis of multi-locus mi- 
crosatellite data, paternity testing tables were set up for 
each progeny individual, and a sequential exclusion pro- 
cedure was applied (Table 2). In each progeny population 
sample, the analysis of 14 microsatellite markers allowed 
immediate identification of: (1) progeny from other seed 
sources (E. grandis mother excluded); (2) progeny de- 
rived from outside pollen contamination (all six fathers 
excluded); (3) selfed individuals (only maternal alleles); 
(4) progeny that belonged to one of the six pollen parent 
(Table 3). Approximately 29% of the offspring was sired 
by pollen parents located outside of the orchard, indicat- 
ing the ineffectiveness of the seed orchard isolation zone. 
No selfed progeny were found in the selected samples, 
which confirmed self-incompatibility of the seed donor 
parent. However, 6 of the 72 randomly taken trees (8.3%) 
and eight of the ten stunted trees (80%) were selfs. The 
sample of chlorotic trees analyzed indicated that these 
were sired either by outside pollen (Table 3) or by pollen 
tree no. 6 (Table 4). 

Pollen parent trees nos. 3 and 5 displayed indistin- 
guishable multi-locus genotypes at all 14 microsatellite 
markers (data not shown). Analyses were repeated three 



Table 1 Microsatellite markers used in the paternity testing of progeny individuals. The linkage group is as determined in Brondani et al. 
(2002) {nd not determined) 



Microsatellite Forward primer 5'-3' 



Reverse primer 5'-3' 



Observed 



Number of Linkage 



marker 


Embra 


6 


Embra 


10 


Embra 


11 


Embra 


16 


Embra 


21 


Embra 


22 


Embra 


27 


Embra 


30 


Embra 


37 


Embra 


40 


Embra 


49 


Embra 


52 


Embra 


53 


Embra 


131 



heterozygosity 


alleles 


group 


0.50 


5 


1 


0.83 


7 


10 


0.17 


4 


1 


0.67 


7 


1 


1.00 


7 


10 


0.67 


7 


11 


1.00 


6 


2 


1.00 


7 


8 


0.83 


9 


5 


1.00 


9 


10 


0.67 


9 


3 


0.67 


9 


7 


1.00 


7 


8 


1.00 


6 


nd 



AgAgAATTgCTCTTCATggA 

gTAAAgACATAgTgAAgACATTCC 

gCTTAgAATTTgCCTAAACC 

CAACgTTCCCCTTTCTTC 

AC A Aggg A A ACTTg ATCg 

gCACATgCACACACggTTg 

ATAACCACACCAATCTgCA 

TTAgTTgAATCCAACCATTg 

CACCTCTCCAAACTACACAA 

AAAgTATCTCCACgCTTCAT 

ATTATTggTTCATATTgAAAACC 

TAATCAgCATTAgCgAAAgA 

ATTAgCTTTTCTgTAACCCg 

ACTTAACATCTATACATATTTg 



gAAAAgTCTgCAAAgTCTgC 

AgACAgTACgTTCTCTAgCTCA 

gTAAAATCCATgggCAAg 

ATgTTAggCCAAACCCAg 

ggAACCgAACATAgCAAg 

AAggCCAgTggTCgTgAgTC 

TATAgCTCgAACgCTCAAC 

TATATAAggTgCAAATAATAAA-Cg 

CTCCTCTCTCTTCACCATTC 

TCCCAATCATgATCTTCAg 

AgATAgAgATTgAgTgAgACCC 

CgTATATgTTCAgAgTCAATCC 

gAATggACAAgCTCTgATg 

TgTCCTATCTggCTCA 
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Table 2 Example of paternity testing of a selected progeny 
individual with six microsatellite markers. Alleles were numbered 
from the smallest to the largest sized allele. At each locus, the 



maternal allele is indicated in italic and the obligatory paternal 
allele is indicated in bold and underlined. Pollen parent tree no. 6 
was declared the father of individual 46 



EMBRA locus 


Mother tree 


Progeny individual 46 


1 


2 


3 


16 


4/5 


2/5 


1/3 


1/1 


1/2 


37 


1/2 


2/6 


7/7 


6/9 


5/8 


6 


3/5 


111 


3/5 


?/4 


4/4 


131 


4/4 ' 


J/4 


1/6 


1/6 


1/5 


27 


■ 1/5 


7/6 


2/5 


2/4 


2/4 


10 


3/5 


5/7 


5/5 


2/5 


1/5 



4/4 
8/8 
3/4 
3/6 
5/5 
5/5 



2/4 
7/8 
4/4 
1/5 
4/4 
1/4 



2/3 
5/6 
T/4 
1/5 
7/6 
5/7 



Table 3 Parentage determination in the five progeny samples 



Progeny sample Number of Non-maternity Sired by outside Selfed Sired by orchard 

individuals (seed mixture) pollen parents pollen parents 

Selected-plus trees (SS A) 72 1 . 17 0 54 

Selected-plus trees (SSB) 72 4 23 0 45 

Non-selected trees (RNS) 72 14 22 6 30 

Chlorotic trees 30 1 11 0 18 

Stunted trees 10 0 0 8 2 

Total 256 20 73 14 149 



Table 4 Number of offspring sired by each pollen parent in the 
progeny samples from the seed orchard 



Progeny sample Pollen parents 





1 


2 


3 


4 


5 a 6 


SSA — Selected sample A 


32 


11 


0 


0 


11 


SSB — Selected sample B 


29 


1 


1 


0 


14 


Non-selected trees 


5 


8 


0 


0 


17 


Chlorotic trees 


0 


0 


0 


0 


18 


Stunted trees 


0 


2 


0 


0 


0 


Total 


66 


22 


1 


. 0 


60 



a Pollen parents nos. 3 and 5 are the same tree 



times with leaf samples collected from different clonal 
ramets; these confirmed the results. A probability of a 
random match between trees nos. 3 and 5 was estimated 
to be 1 in 89 billion using allele frequency estimates for 
these loci for E. grandis and E. urophylla (R. Brondani 
and M. Kirst, unpublished results). 

Pollen parents nos, 1, 2 and 6 sired 148 out of the 149 
offspring sired with orchard pollen. Pollen parents nos. 3 
and 5, actually the same tree, sired only a single offspring 
in SSB, and parent no. 4 did not sire any offspring 
(Table 4). For parents nos. 1 and 6, significant differences 
were found in the proportions of progeny individuals in 
both selected samples versus the non-selected sample of 
trees. However, while pollen tree no. 1 sired significantly 
more selected progeny than non-selected ones in both 
SSA and SSB (P=0.0002 and 0.00005 respectively), pol- 



len tree no. 6 sired significantly more offspring in the 
non-selected tree sample for both comparisons (P= 
0.00141 and 0.03365, respectively) (Table 5). For pollen 
tree no. 2, no significant difference was observed between 
the proportions of sired offspring in SSA versus RNS (P= 
0.58942), but significantly more RNS offspring were 
sired as compared to SSB (P=0.00221) (Table 5). A 
similar trend was observed when the test was carried out 
analyzing the outside pollen contribution — i.e. either the 
outside pollen contribution sired more non-selected off- 
spring in SSA (22 vs. 17, Table 3) (P=0.00028), or no 
difference was detected in SSB (P=0.09132) (Table 5). 

For practical breeding purposes, since parents nos. 1, 2 
and 6 displayed significant reproductive success and 
effectively contributed to the generation of superior hy- 
brid trees, they were kept in the orchard. Parents nos. 3, 4 
and 5 did not contribute to the generation of superior 
offspring and were thus eliminated. A significant differ- 
ence (P 0.01) in volume growth was observed between 
forest stands produced from orchard seeds before and 
after culling the worst parents (Table 6). The new seed 
variety (MAI = 44.5±5.8 m 3 /ha per year) resulted in forest 
stands that were on average 24.3% more productive in 
volume growth than the original seed variety (MAI = 
35.8±2.3 m 3 /ha per year). 



Table 5 Results (P-values) of p o]len source SSA— Selected sample A SSB— Selected sample B 

Fisher exact contingency tests — 

for binomial proportions of Pollen tree no. 1 0.00020 0.00005 

sired progeny individuals in the Pollen tree no, 2 0.58942 0.00221 

selected and non-selected off- Pollen tree no. 6 0.00141 0.03365 

spring samples for the different Outside pollen 0.00028 0.09132 

pollen sources 
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Table 6 Number of commercial stands (average 20 ha, with 1,11 1 
trees/ha) submitted to inventory. Mean annual increment (MAI) in 
cubic meters per hectare per year is given as the mean and standard 



deviation for seed varieties produced before and after elimination 
of parents with a low reproductive success for revitalization of the 
orchard 



Seed variety 



Number of stands 



MAI (mean) 



MAI (standard deviation) 



Old (before culling) 
New (after culling) 



78 
18 



35.8 
44.5 



2.3 
5.8 



Discussion 

A battery of 14 polymorphic markers was selected to 
carry out parentage tests of 256 progeny individuals. 
Maternity checks were carried on all individuals and, 
following the exclusion of seed mixtures and selfed in- 
dividuals, paternity tests were carried out. Markers were 
specifically selected to allow easy scoring in polyacryla- 
mide gels and to maximize allelic differences between the 
group of five pollen parents and between them and the 
maternal parent. This was possible as the five pollen par- 
ents were unrelated and the maternal parent belonged to a 
different species. We also specifically looked for markers 
for which all of the parents, maternal and pollen parents, 
were heterozygous so as to avoid the occurrence of 
null alleles that could result in false paternity exclusions 
(Moller 1995). A relatively large number of alleles were 
observed for the microsatellite markers used in the study, 
and for 6 out of the 14 markers all the parents were het- 
erozygous, thereby allowing a significant discrimination 
power and confidence in the parentage testing (Table 1). 

We also were interested in markers for which the 
maternal parent did not share alleles with any of the six 
pollen parents. Out of the 14 markers, three had this 
configuration. Based on these three markers the analysis 
of maternity versus non-maternity — i.e. seed mixture — 
was very fast and efficient. Progeny individuals that did 
not have any of the two maternal alleles were declared to 
result from the seed mixture. These loci were also very 
useful for determining the presence of selfed progeny 
individuals. The presence of homozygous loci in the 
maternal tree would also be efficient for the maternity 
checks, however null alleles could lead to false maternity 
exclusions. 

The analysis of pollen contamination in the orchard 
was based on multiple paternity exclusions when the six 
pollen parents were tested. In this procedure it was 
essential to have genotypes for several loci to declare 
non-paternity with absolute confidence with a minimum 
of three excluding loci to avoid false exclusions due to 
mutations (Gunn et aL 1997). Di-nucleotide repeat mi- 
crosatellites have been reported to have higher mutation 
rates following the step-wise mutation model (Valdes 
1993) than longer repeat motifs in humans (Brinkmann et 
al. 1999) and, more recently, in maize (Vigouroux et al. 
2002). Because no data are available yet on the frequency 
of null alleles for Eucalyptus microsatellites and, in fact, 
for the vast majority of plant species, our premise was that 
such a behavior would be the case for Eucalyptus as well. 
Although the most informative six markers for which all 



trees were, heterozygous could be sufficient, we used an 
over-abundance of marker loci in the paternity testing to 
avoid false exclusions due to the inheritance of mutated 
alleles from the pollen parent to progeny individuals. 

When we excluded the offspring resulting from pollen 
contamination and took into consideration the fact that 
pollen trees nos. 3 and 5 were clonal (i.e. identical), a total 
of 149 progeny individuals were tested against the five 
unique pollen parents, making a total of 745 paternity 
tests. In these 745 paternity tests, a putative occurrence of 
a null allele was observed at locus EMBRA6 where pollen 
parent no. 1 was apparently homozygous. However, 28 
out of a total of 66 offspring lacked the obligatory 
paternal allele. At all other loci offspring and pollen 
parent shared the paternal allele correctly. This result 
strongly suggests the occurrence of a null allele in pollen 
parent no. 1 that, as expected, was transmitted to ap- 
proximately half (28/66=42%) of its progeny. This locus 
was not considered in the paternity analysis for pollen tree 
no. 1. 

A total of 745x14= 10,430 allelic transmission anal- 
yses were carried out in the paternity testing. There were 
five putative occurrences of paternal mutations — i.e. 
where the obligatory paternal allele was different, usually 
by a single dinucleotide step, from the allele observed in 
the pollen parent. Again, at all other loci, a full consis- 
tency of obligatory paternal allele was observed. The 
resulting mutation frequency estimate of 4.8xl0" 4 is very 
similar to recent estimates in maize where the mutation 
rate per generation was estimated to be 7.7x1 0" 4 for 
microsatellites with dinucleotide repeat motifs (Vigour- 
oux et al. 2002). 

When the inconsistencies due to mutations and null 
alleles were not counted, the sequential procedure of 
paternity exclusion proved to be a very efficient and 
accurate method to quickly exclude non-fathers at several 
marker loci. By this approach, once a pollen parent was 
excluded in at least three loci, the remaining paternity 
tests were only carried out with the non-excluded fathers. 
Typically, however, all non-fathers were excluded in the 
analysis of the first five more informative loci. The 
proposed approach is therefore very efficient and should 
not demand the analysis of a large number of loci. It is 
obvious that as the number of mother trees and pollen 
parents increases and the relatedness of these trees also 
increases, confident paternity testing will require an 
increased number of polymorphic loci to reach conclu- 
sions as multiple males are often found to be genetically 
compatible with each offspring tested, even when the 
probability of excluding an unrelated male is high. In such 
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complex parentage testing situations, the deterministic 
approach used in this study is usually not possible, and 
likelihood based paternity inference methods become 
necessary (Meagher 1986; Marshall et al 1998; Chaix et 
al. 2003). 

A significant proportion of the analyzed offspring, 
approximately 29%, was sired by pollen parents located 
outside of the orchard. Campinhos et al. (1998) estimated 
a lower contamination rate of 14% in the same orchard 
based on isozyme markers, possibly due to limited in- 
formative polymorphism that would allow correct paren- 
tal discrimination or simply due to a sampling effect. 
Outside pollinator trees did not contribute significantly to 
the generation of superior offspring individuals (Table 5), 
although one cannot preclude the possibility of some 
occasional superior trees deriving from pollen outside the 
orchard. Heavy pollen contamination has been observed 
in a number of studies with wind-pollinated conifer spe- 
cies (El-Kassaby and Ritland 1986; Harju and Nikkanen 
1996; Pakkanen et al. 2000; Moriguchi et al. 2002), 
typically reducing the expected gains from seed orchards. 
For eucalypts, a number of studies have estimated the 
preferential mating system in natural populations and seed 
orchards using isozymes (e.g. Moran et al. 1989) and, 
recently, a complex pattern of mating was described in an 
E. regnans seed orchard in Australia. Gene dispersal was 
influenced by crop fecundity and orchard position of the 
mother trees with approximately 50% of effective pol- 
len gametes coming from males more than 40 m away 
from mother trees, indicating that insect pollinators are 
efficient promoters of cross-fertilization (Burczyk et al. 

2002) . Furthermore, in a recent outcrossing rate study 
carried out in an E. grandis orchard in Madagascar, a 
pollination rate from outside the seed orchard of 39.2% 
was estimated based on six microsatellite markers (Chaix 
et al, 2003). The results from our study also indicate that 
in our exotic conditions the 300-m-wide belt of native 
tropical forest maintained around the orchard to provide 
genetic isolation has not been effective in preventing 
insect-mediated gene flow from nearby eucalypt planta- 
tions. 

Pollen parents nos. 1, 2 and 6 displayed the highest 
male reproductive success, siring 148 out of the 149 
offspring analyzed. All three pollen parents were thus 
maintained in the seed orchard. Pollen trees nos. 3, 5 and 
4 did not father any offspring and were subsequently ex- 
cluded. The absence of offspring from these parents is 
probably due to asynchronous flowering in relation to the 
maternal tree as well as temporal variation in the flower 
crop, suggesting that physical proximity between eucalypt 
tree crowns not necessarily implies successful mating in a 
seed orchard. Sampling effects may also have contributed 
to a certain extent. Differential male reproductive success 
was also observed in a recent E. grandis seed orchard 
study, where only 199 out of 349 potential male trees in 
the seed orchard contributed to the pollination of' 440 
offspring and at a very variable siring rate (Chaix et al. 

2003) . 



Pollen tree no. 1 sired significantly more selected 
progeny than non-selected ones in both samples. For the 
other two pollen trees, either no significant difference was 
observed between the proportions of sired offspring in the 
selected versus non-selected samples or, in the case of 
pollen tree no. 6, significantly more non-selected trees 
were sired (Table 5). These results indicate a higher SCA 
and flowering synchrony of pollen tree no. 1 with the 
maternal tree, ensuing the immediate establishment of a 
bi-clonal seed orchard composed exclusively of pollen 
tree no. 1 and the E. grandis maternal tree. 

Although a bi-clonal seed orchard could have been 
established, involving only pollen tree no. 1 and the E. 
grandis seed parent, due to spatial distribution of the 
remaining pollen donors and seed parent ramets and with 
the objective of maintaining a slightly larger effective 
population size, it was decided to leave all three pollen 
parents and to cull only those that clearly did not 
contribute to the generation of offspring. It seems there- 
fore unreasonable to have accrued a 24% gain in volume 
growth leaving the same pollen donors, i.e. nos. 1, 2 and 
6. However, besides culling the worst parents, a signif- 
icant improvement was made in the silvicultural man- 
agement of the orchard so as to minimize external pollen 
contribution and to improve flowering of the remaining 
pollen donors. Intensive fertilization regimes were ap- 
plied to the orchard so as to stimulate abundant flowering, 
and a number of bee cages were introduced in nhe seed 
orchard so as to keep the bees from flying longer dis- 
tances to harvest pollen and nectar. These measures were 
likely successful in minimizing external pollen contribu- 
tion and at the same time increasing pollen contribution 
from the pollen donors in the orchard. Besides the re- 
vitalization of the orchard, it was also noted that the 
number of grafted ramets of pollen donors nos. 1, 2 and 6 
was not balanced due to tree death since the time that the 
seeds used to. establish the commercial plantation studied 
were harvested. The census number of ramets of pollen 
donor no. 1 was slightly larger that of pollen donor nos. 2 
and 6, possibly due to its more adapted growth and vigor. 
A potentially larger contribution of pollen parent no. 1 to 
the new seed crop from the orchard also contributed to the 
increased volume growth. 

As expected, no seifed progeny individuals were found 
among the selected trees. However, selfs were found in 
the random non-selected trees at a rate of 8.3%, which is 
in agreement with the selfing rates typically estimated in 
eucalypts (Griffin and Cotterill 1988; Gaiotto et al. 1997). 
This same selfing rate should therefore be expected even 
in a biclonal seed orchard. Not surprisingly, however, 
80% of the stunted trees were selfs, thereby confirming 
the inbreeding depression effect on growth typically 
observed in controlled selfing experiments in eucalypts 
(Griffin and Cotterill 1988; Hardner and Potts 1995). 
Finally, no selfs were found in the sample of chlorotic 
trees. The paternity testing results suggest a specific 
association between pollen parent no. 6 and the appear- 
ance of leaf chlorosis, although 1 1 out of the 30 chlorotic 
trees were sired by outside pollen parents. The establish- 
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ment of a more efficient barrier to outside pollen 
contamination should minimize the manifestation of this 
symptom. 

Our study demonstrates that a marker-assisted breed- 
ing tactic involving parentage testing with microsatellite 
markers can efficiently and rapidly be applied to existing 
seed orchards with the specific objectives of verifying the 
differential male reproductive success of pollen parents 
and identifying those that successfully generate superior 
offspring. It should be pointed out, however, that the 
tactic applied in this study represents an alternative 
measure, in situations where time is the critical issue, the 
breeder is faced with an emergency demand for improved 
seeds from existing orchards and no data are available on 
progeny testing of trees but rather only on commercial 
plantations derived from the seed orchard. This tactic is 
certainly no substitute for a well-conducted progeny trial 
followed by estimation of parental breeding values for 
backwards selection. We have shown, however, that this 
retrospective selection procedure coupled to improved 
management of the seed orchard did in fact result in an 
improved seed variety displaying a significant gain in 
volume growth over the average performance of the seed 
lots derived from the original unselected orchard. 

This short-term marker-assisted breeding tactic should 
be applicable to the improvement of old seed orchards 
currently supplying planting stock of other forest tree 
species. Although this demonstration experiment was car- 
ried out on a relatively small sample of parent trees, giv- 
en an adequate screening and selection of a battery of 
informative microsatellite markers, sufficient power of 
discrimination can easily be attained to resolve paternity 
and maternity of offspring derived from orchards con- 
sisting of several tens of unrelated parents. Relatedness 
among the orchard parents could be a problem when a 
limited number of markers are available. Lambeth et al. 
(2001), after genotyping a set of 45 parental trees at seven 
chloroplast and three nuclear microsatellite loci, found 
that more markers would be necessary for unambiguous 
paternal determinations of progeny from a complete 
pollen mix due largely to relatedness in the pine popu- 
lation studied. 

Finally, this tactic, as opposed to the conventional 
approach of estimating predicted gain based on controlled 
pollination and progeny trials, measures realized gain in 
an operational setting, which is a function not only of 
the SCA of pairs of parents but also of the effective 
reproductive success of mating between them, which is 
often a much more critical variable when large supplies of 
improved seeds are desired. 
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Population Management in the Age of Genomics 

Rowland D. Burdon 1 and Phillip L.Wilcox 1 

Abstract : Continued population improvement will remain the foundation for future 
genetic gain in forest trees. Any gains that can accrue directly from genetic 
engineering, and be captured through propagation technology, will be superimposed 
on those from population improvement. Efficient population management will need 
to be directed not only at pursuing recognised breeding goals, but also in preparing 
for addressing new, unpredictable breeding goals. Both these goals are served by 
maintaining genetic diversity. 

Genetic diversity depends initially on the population base. For characterising 
the base, genomics is a powerful adjunct, in conjunction with common-garden 
experiments. Within populations there are various possible measures of genetic 
diversity, involving polymorphisms for neutral markers and functional genes, and 
numbers of alleles and allele frequencies. Calibrating marker diversity against 
functional diversity is not simple, but will become easier with functional genomics in 
which the phenotypic impacts of nucleotide polymorphisms may eventually be 
identified. A need exists to counter the tendency for continued selection to erode 
diversity. 

An associated, but not identical issue is controlling levels of inbreeding, which 
is not wanted in commercial material and yet may serve as a breeding tool. A 
customary approach is to pursue intensive ongoing genetic gain in a pedigreed 
breeding population underpinned by very broadly based gene resources. Use of 
genomics can relax some requirements for population management, notably the need 
for pair-crossing in order to maintain pedigree. Related to that, genomics can serve 
to retrieve historical losses of pedigree information, even overturning the tenet of 
production populations being genetic dead-ends. It can also correct potentially costly 
misidentification of material. 

Genomic information has the widely hailed application of marker-assisted 
selection, although this has encountered various problems with forest trees. It is a 
promising tool in the quest not just for disease resistance but resistance that is 
durable. This requires a diversity of resistance mechanisms, yet the presence of 
some resistance genes can mask the presence or absence of others, unless the latter 
can be detected genomically. 

With future knowledge of functional polymorphisms the new understanding of 
the nature and origins of genetic variation should allow far more efficient 
management of populations in order to both exploit genetic diversity and preserve it 
for the longer term. 

Keywords : Population management, genetic diversity, genomics, DNA markers, 
pedigree 
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INTRODUCTION 

Genomic information on forest trees can be used for several purposes, apart from being part of 
the advancing front of fundamental knowledge. Of main interest for this paper is how the 
availability of genomic research tools and the resulting information can serve the underlying 
objectives of population management which, in turn, is all about capturing cumulative genetic 
gain while maintaining the genetic diversity that safeguards the future genetic gains. For 
capturing genetic gain, genomic information can be used for marker-assisted or marker-based 
selection which exploits linkage disequilibria between assayable polymorphisms and ones that 
underpin phenotypic variability. Such selection can be refined into gene-assisted or gene-based 
selection if it has become possible to identify actual nucleotide polymorphisms that govern 
quantitative effects on traits. Much has been written about the issues of refining and quantifying 
the efficiency of selection based on use of genomic information, but for this paper that is of 
interest mainly for its ramifications for appropriate structuring of populations. Genomic 
information also has value for safeguarding the capture of genetic gains by verifying genetic 
identity. Of prime interest for this paper are the various more direct applications of genomic 
information for appropriately structuring populations. 

The need for genetic diversity in forest tree breeding programs has never been seriously 
challenged, even if it has not always been honoured in how improvement programs are 
implemented. Long-term gain will depend on genetic diversity, which will tend to be run down 
by the process of intensive selection and associated maintenance of very finite breeding 
populations. Such limitations in long-term response to may not arise at all rapidly provided the 
same breeding goal continues to be pursued. However, breeding goals can change rapidly, 
owing to changes in market perceptions or biotic developments such as the emergence of 
significant fungal diseases or animal pests. Such changes can be substantial, and place heavy 
demands on available genetic diversity. While genetic engineering offers the prospect of 
conferring attributes that are not programmed by any of the available genetic diversity within 
species, the contributions of genetic engineering are widely seen as being properly superimposed 
upon a platform of genetic improvement resulting from classical breeding. For instance, 
engineering a new function from scratch can take a long time, and be a very uncertain venture 
compared with capturing the attribute from straightforward expression of existing genetic 
variation. On the other hand, genetic engineering technologies could also confer a narrow suite 
of attributes faster than existing breeding approaches can, although this may well impose added 
pressures of reduction in diversity, particularly in deployed populations. 

Even without a specific or pressing problem, obtaining a meaningful measure of genetic diversity 
is highly desirable. Diversity can derive both from diverse population origins and within- 
population variation. Genomics can be used to indicate origins, of material of mixed ancestry. 
Traditional taxonomic tools, including chemotaxonomic ones, can be effective in this area, but 
genomic information has major advantages through often requiring limited tissue samples, 
providing quick results, and being independent of the multifarious environmental influences. 
Information on population origins needs to be complemented by information from common- 
garden field experiments which yield information on a range of traits that include growth 
potential, environmental tolerances, biotic resistances, and wood properties. 
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MANAGEMENT ISSUES 

The issues will be addressed under the following categories: 

• Characterising the genetic base (sorting out taxonomy and crossability; identifying ancestral 
populations, migration history, and past adaptive pressures) 

• Quantifying within-population diversity 

• Managing inbreeding and relatedness 

• Selection using markers 

• Assuring durable disease resistance 

• The future (in the era of functional genomics). 

GENETIC BASE 

Knowledge of the genetic base that gave rise to one's population is always welcome, although 
not necessarily crucial or easily obtained. In the first instance, one may want to clear up the 
taxonomic status of one's material. Crossability, while an important but not absolute taxonomic 
criterion, is of immediate practical importance and something that needs to be addressed 
empirically. For addressing the general taxonomic issues, there are various classical herbarium- 
based criteria, which have come to be supplemented by phytochemical and cytological evidence, 
but genomic analysis is developing into a powerful new tool for the purpose. Even without a 
traditional taxonomic issue, it can be valuable, for instance, to know the geographic origins of a 
land-race population. For instance, the effects of sub-optimal origin may have been masked by 
either a shake-out of the neighbourhood inbreeding of natural populations, or by adaptive 
responses to selection pressures imposed by the exotic environment. Alternatively, a knowledge 
of the Quaternary migration history may warn one that a local population is likely to be 
significantly sub-optimal for commercial purposes. For addressing such issues, common-garden 
field experiments are likely to provide the best information, but genomic information may often 
provide a powerful adjunct, with the advantage of being far quicker to obtain. 

Attempts to characterize different populations or even interfertile species for their collective 
diversity, without the benefit of common-garden field experiments, can be perilous. Apparent 
affinities among such entities can depend very strongly on the set of genomic markers chosen. 
Besides, the inherently neutral nature of many genomic markers, whatever their power for 
certain purposes, means that there is no a priori reason why they should indicate values for the 
adaptive or commercial traits that are of direct interest to the breeder. In a case that will not be 
named, a genomic study showed a particular tree species to be undistinctive with respect to the 
markers used. From that, the suggestion was made that the species could forgotten about for 
purposes of genetic conservation. This flew in the face of common forestry knowledge that the 
species concerned is highly distinctive in respect of certain ecological tolerances and silvical 
properties. 

Indeed, indications are that neutral DNA markers can often be poor indicators of adaptive 
differences among populations within species (e.g. Karhu et al. 1996). 
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Where there is interest in preserving populations in a genetically pure state, genomics can serve 
as a new and powerful tool for detecting genetic contamination, quantifying it, and even culling 
individuals resulting from the contamination. This, however, will depend on identifying DNA 
markers that differentiate well between the populations in question. 



ISSUES OF QUANTIFTING WITHIN-POPULATION DIVERSITY AND 

IMPLICATIONS 

Apart from population origins being a source of diversity that can be both identified and 
characterized, there is the within-population diversity, which may often be the bigger general 
issue. There remains the issue of whether marker diversity parallels the functional diversity that 
is of real interest to the breeder. But there are also issues of what are the most meaningful 
measures of diversity for markers and functional traits respectively. 

Various possible measures of diversity exist, and they can have different implications according 
to the intensity of selection and the time frame of interest. For both markers and fimctional 
alleles there are numbers of alleles at the polymorphic loci, the expected heterozygosity over 
such loci which will often be dominated by intermediate allele frequencies, plus the numbers of 
loci involved. For functional alleles, the magnitudes of their effects are also of obvious 
importance. 

Numbers of alleles per locus and expected heterozygosity are two, largely complementary 
measures. Numbers of alleles are readily influenced by sampling effects. As such, they are very 
subject to founder effects or genetic drift, which can both occur incidental to intensive selection. 
They will also be affected by sample sizes in population surveys. 

Percentages of functioning genes that are polymorphic are still conjectural, at least for the typical 
outbreeding forest tree species. Nevertheless, this is an area where knowledge will doubtless 
make rapid advances (see later). 

Potential response to selection depends on functional diversity. The immediate functional 
diversity, and consequent short-term responses to selection, will tend to be governed by alleles of 
intermediate frequency - unless there is both very intensive selection and alleles of large effect 
for the trait(s) under selection. Longer-term selection responses, however, may depend more on 
alleles that are currently at low frequencies but will only contribute significantly to expressed 
variation when their frequencies increase. The ways in which low-frequency alleles contribute to 
functional diversity will be especially sensitive to the degree to which the effects of such alleles 
are recessive or dominant. 

In functional genes, many polymorphisms are silent, since they involve synonymous codons. 
Considerable discrimination is therefore possible between identity by function and identity by 
descent, which may offer important insights into the way in which different populations may 
have come to diverge. 
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Another issue is to what extent the functional genetic variation is governed by polymorphisms in 
the coding or promoter regions of genes. Some recent findings for plant species other than forest 
trees (e.g. Morgante & Salamini 2003, Paran & Zamir 2003) have indicated that large phenotypic 
variations can be governed by polymorphisms in the regulatory region rather than in the coding 
regions of the genes concerned. Such genetic variation will therefore not be evident in protein 
polymorphisms, although it could be detected via variation in transcript- protein- and/or 
metabolite concentrations in specific tissues. 

MANAGING INBREEDING AND RELATEDNESS 

Apart from running down genetic diversity, long-term maintenance of closed, finite breeding 
populations inevitably incurs some level of inbreeding. Inbreeding is not necessarily the same 
thing as loss of diversity, because deliberate inbreeding can be practised to both maintain 
diversity and preserve options for subsequent outcrossing. Deleterious effects arising from 
inbreeding, namely loss in fitness and/or reproductive capacity via the effects of lethal or semi- 
lethal loci, could be managed using information from markers linked to such loci (e.g. Kuang et 
al. 1998). However, costs involved in establishing populations for detection of marker-trait 
associations are likely to be considerable relative to the gains from such selection, particularly in 
breeding populations. Nonetheless, as coancestry builds in such population, the imperative is 
likely to increase for more proactively dealing with the problem. Options range from avoiding 
certain allelic combinations to actively purging lethal genes. 

For both managing inbreeding and maintaining genetic diversity, maintaining full pedigree 
information is extremely valuable. However, the traditional way of maintaining full pedigree has 
its costs, both direct and in achieving some of its aims. Making controlled pair-crosses, and 
keeping pair-crosses properly identified in replicated field experiments, are expensive operations, 
especially if the reproductive biology of a species makes pair-crossing difficult. In some cases, 
open pollination may even be preferred. Also, the size of fully pedigreed breeding populations 
may be constrained to levels at which loss of low-frequency alleles can become significant 
(Burdon 1997). 

Even where full pedigree is nominally kept, misidentification can occur in various ways. That 
can compromise genetic gain, the more so the more advanced a breeding program is. It could 
also compromise control of inbreeding, and even maintenance of genetic diversity. However, 
simple-sequence repeat (SSR), or microsatellite, markers are proving a very powerful tool for 
verifying identity of both individual clones and parentage of progenies. Indeed, the 
misidentification rates being thus detected are quite embarrassing - something that is admitted 
much more freely by word of mouth than in print. 

Genomic information, however, can allow retrieval of pedigree information. Lambeth et al 
(2001) has advocated the use of DNA markers in order to allow retrieval of full pedigree 
information after saving costs by making polycrosses instead of pair-crosses. Earlier, Burdon 
(1997) had proposed use of DNA markers be used to make good a complete lack of pedigree 
information if it was necessary to select from unpedigreed commercial stands. This eventuality 
runs counter to a traditional tenet of commercial stands being a genetic dead-end in population 
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management. It could, however, arise in the case of a biotic crisis whereby only very rare trees 
were resistant to a new disease or pathogen strain. Such resistant individuals might be products 
of very rare favourable de-novo mutations, and finding them would likely depend on the sheer 
numbers of genotypes that would only be available in large commercial plantings. Pedigree 
reconstruction would allow the breeder to avoid selecting resistant trees all from one or a very 
few pedigrees with the attendant risks of both inbreeding and loss of genetic diversity. Having a 
finite set of possible or likely parents, as may be the case in seed orchards, would be a big 
advantage. That said, the challenges of achieving complete pedigree reconstruction may be 
formidable, but even partial success may suffice for achieving the underlying objective. Such a 
scenario may seem fanciful, but biotic crises have affected some very important tree species. For 
the breeder to tackle a biotic crisis head-on would also require powerful propagation technology 
and/or the lack of alternative species that are strongly competitive commercially. 

Despite the obvious challenges in developing a sufficiently powerful set of markers for complete 
pedigree reconstruction, the sort of SSR marker diversity discovered in pines (Karhu et al. 2000) 
and and Douglas-fir (Slavov et al. 2004) augurs well. 

SELECTION USING GENETIC MARKERS 

The potential benefits of using genomics to detect and use quantitative trait loci (QTL) for 
marker-assisted or marker-based selection (collectively referred to as MAS) are well recognised. 
It is widely agreed that there are tree general areas where MAS could be especially beneficial: 
increasing selection intensity, earlier selection, and cheaper selection for expensive-to-measure 
traits. However, various methodological traps have come to light (e.g. Ball 2001), and much had 
already been written about both the potential advantages and the pitfalls of MAS (e.g. Strauss et 
al. 1992, Johnson et al. 2000). The magnitudes of QTL effects are subject to estimation errors, 
and some can be greatly overestimated, in what is widely known as selection bias. However, 
experience with annual (or other short-term) agricultural crops has indicated that where QTL of 
large effects are involved, as may happen with introgressing genes from wild relatives into 
domesticated crops, some major payoffs can accrue (cf Paran & Zamir 2003). 

Overall, MAS can be a selection tool for both breeding and production populations. While use 
of MAS as such is in no way central to population management, questions have been raised 
concerning the possibility that its intensive use would incur risks of otherwise avoidable loss of 
potentially useful alleles. On the other hand, as information on functional loci becomes 
increasingly available, tools to track specific alleles will become more and more effective, 
facilitating quantification and/or avoidance of allele loss during selection. 

With forest trees, which are typically outbreeders, there is usually not the linkage disequilibrium 
between non-genic marker alleles and QTL that is required for MAS. However, with fusion of 
well-differentiated populations, or with advanced hybrids between species (as in domesticated 
apple - Bus et al. 2000) - there may well be cases of the desired combination of large-effect QTL 
and strong linkage disequilibrium. Another exception is within known pedigrees, for which 
linkages have to be identified in each individual case. This area has been the focus of much 
effort over the past decade, driven in part by the development of within-family marker-assisted 
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selection. Such research has also resulted in the development of marker tools that can reveal 
much about population structure and existing inbreeding levels. However, these marker systems 
are largely based on polymorphisms in non-gene-associated regions, and may not adequately 
reflect the underlying functional variation. 

More recently, research has changed focus from pedigree-specific populations, to non-pedigreed 
association genetics. One of the reasons for this is the desire to develop the option of among- 
family selection. The marker systems concerned would be more likely to be based on expressed 
genes that have a suspected role in trait variation. The benefit of assaying polymorphisms within 
and associated with expressed genes, together with non-structured populations for ascertaining 
marker-trait associations, is that key parameters such as gene and nucleotide diversity, 
heterozygosity, inbreeding and coancestry relationships can be calculated (Neale & Savolinen 
2004). Therefore, the experimentation required for gene-assisted among- and within-family 
selection has the added benefit of revealing key information regarding those experimental 
populations. 

A key constraint to the implementation of markers as selection tools, is the substantive up-front 
costs required to detect marker-trait associations (Johnson et al. 2000, Wilcox et al. 2001). As 
such the more likely implementation is in deployed populations as opposed to breeding 
population advancement, because of earlier commercial returns. If marker-trait associations are 
used for selection in the breeding population, a build-up of coancestry may be hastened because 
of unequal selection rates in specific pedigrees. Such a build-up would need to be judiciously 
countered in order to maintain diversity. 



DISEASE RESISTANCE (DURABILITY) 

A related issue is ensuring durability of disease resistance against genetic shifts in the pathogens 
(e.g. Burdon 2001). While there is no certainty that lack of durability will be a problem in any 
particular case, measures to ensure durability seem prudent. Durability typically appears to 
depend strongly on a diversity of resistance mechanisms. For large-scale deployment of 
individual clones, pyramiding a multiplicity of resistance factors within individual clones is a 
likely need. For heterogeneous populations, it appears that population resistance is not dependent 
on a multiplicity of resistance factors all being present in every individual. Either way, however, 
knowledge of the various resistance factors, how they operate, and how to detect their presence, 
is extremely valuable. It is one of the problems with disease resistance that some of the more 
spectacular forms of resistance are vulnerable to pathotype shifts and yet tend to mask the 
presence of other resistance factors that may be needed in order to ensure durability. For 
addressing this general problem, genomic information when combined with infection studies in 
genetic experiments can be extremely powerful, and its use is progressing rapidly with annual 
crops. This is partly because there are often individual resistance genes of substantial effects in 
situations where there are superficial appearances of polygenic resistance. With the loblolly 
pine/fusiform rust pathosystem progress is being made in identifying genetic factors involved in 
resistance (Wilcox et al. 1997, Schmidt 2003). On the one hand, the short period of the pine host 
susceptibility during the rotation may make for a relatively simple pathosystem, and one in 
which the pine is not unduly vulnerable to rapid pathotype shifts (Schmidt 2003); on the other 
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hand, acquiring good data is relatively slow and difficult. In the case of leaf rusts affecting 
poplars more rapid research progress is apparently being made (e.g. Yin et al. 2004). Apart from 
poplars having come to be treated widely as model species, the biology of the rusts is conducive 
to much more rapid acquisition of data. 

With the search for resistance to blister rust on North American white pines much progress has 
been made by classical genetic studies. Various individual resistance genes have been identified, 
and the nature of their action (e.g. dominance or recessiveness) characterised (Sniezko et al. 
2004). This, however, has been a long, drawn-out process. It is to be hoped that, with genomic 
tools available, future progress can be much more rapid, in both research and its applications. 

For deployment of genetically resistant material the diversity of resistance mechanisms can, in 
principle, be used judiciously. Comprehensive pyramiding of resistance genes in individual 
genotypes cannot be expected to come quickly. But, even without such pyramiding, it should be 
possible to deploy mixes of tree genotypes that can be expected to show good population 
resistance. 



THE FUTURE 

Current population management is based largely upon various suppositions, which are often tacit 
rather than explicit, concerning the nature of genetic diversity and the ways in which it can be 
generated. The usual rule of thumb is for production populations to contain a (very) few tens of 
effectively unrelated individuals, breeding populations a few hundred, and gene resources a few 
thousands (e.g. White 1992). It tends to be assumed that mutation is a rare event at any 
individual locus, but that low-frequency alleles are potentially very important. Yet there are 
various long-term selection experiments, in which continued response has occurred for many 
generations in closed populations, which suggests that significant mutational or quasi-mutational 
events may not be rare. Since maintaining breeding populations and gene resources is expensive, 
and likely to become more so as progressive genetic gain increases the opportunity costs of 
maintaining such material, the question of requisite populations sizes and appropriate 
management systems is far from trivial. 

Questions of both population management and capture of genetic gain for operational use will be 
increasingly addressed by extending genomics into the study of functional genes. There are 
various ways in which this can be done. First, the genes need to be identified and located, and 
then their polymorphisms can be studied. Starting with just with information on tight linkages 
still leaves very long sequences to investigate and interpret. However, with libraries of 
expressed sequence tags (ESTs), belonging to genes whose function has been identified in model 
species, there are now much more powerful tools for locating specific genes. The quintessential 
model plant species has become Arabidopsis which, despite not being closely related to conifers, 
shows remarkably high levels of orthology with the full range of forest trees (e.g. Kirst et al. 
2003), even though much of the developmental physiology and morphology appear very 
different. However, with the availability of the full nuclear and organelle genomic sequences of 
poplar (Populus trichocarpa) (http://www.jgi.doe.gov/poplar) , it is likely that this will become 
the model species of choice for certain purposes. Among conifers, there are not only the 
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orthologies but also some remarkably high levels of synteny (e.g. between loblolly and Monterey 
pines). Thus, genes of common functions tend to be located in essentially the same positions on 
the same chromosomes of the respective species, so genomic information obtained from one 
species can be easily used for studying another. With the identity and roles of various genes 
identified, the developmental significance of polymorphisms can then be ascertained. Tactically, 
that can allow the breeder to screen for particular alleles, or to identify alleles that might be 
transferred between species by genetic transformation. 

The value of orthologies for probing the nature of polymorphisms has its limits, however, since 
the roles of particular genes in governing trait phenotypes can vary among plant species. 
Experimentally, the roles of different base-sequence alleles in governing phenotype can be 
probed or verified by using genetic transformation; however, some of the information may take a 
long time to come in with trees, particularly for traits that are expressed only at later ages and 
where no early indicator characteristics can be used. 

Various strategic implications will stem from the basic information on the levels and nature of 
polymorphism, on the dynamics and nature of mutation and on gene transposition. Knowing the 
level and nature of polymorphisms involving gene function, as distinct from recent descent, 
should give better insights into the significance of numbers, at different stages in the life cycle, in 
population management. 

We have various implicitly or explicitly pedigree-based measures of effective population size, 
which serve as indicators of changes in genetic diversity. These measures include Status 
Number (Lindgren et al. 1996), as well as more traditional measures (Caballero 1994). Although 
we have some intuitive appreciation of the value of different measures for different purposes, 
their exact practical significance is usually very uncertain. Hopefully, with knowledge of the 
level and nature of functional polymorphisms, and of mutational processes and transpositions, it 
will be possible to achieve not only better and more direct measures of genetic diversity but also 
much better interpretation of the measures of effective population size. 



CONCLUDING 

For many of those who are charged with population management, or see it as their responsibility, 
a strongly precautionary approach is indicated in the present state of genomic knowledge. At the 
same time, breeders may have to deal with people who want to see guarantees of returns from the 
direct costs and opportunity costs involved in population management, which can place them in a 
difficult situation. Better knowledge will therefore be extremely welcome. As yet we have no 
silver bullet, and we may never be able to majorly relax the requirements for sizes of the 
populations needed for continuing advance of breeding programmes. Nonetheless, the very rapid 
progress that is now occurring in functional genomic promises a much better basis for future 
decision-making. 
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Summary 



Polymix (PMX) breeding of interspecific hybrids of forest trees promises to greatly reduce 
the cost of controlled crosses, while increasing the number of hybrid genotypes that can be 
evaluated. However, severe pollen competition may make this approach impractical, as 
dominant pollen parents may not necessarily produce superior hybrid offspring. In this 
study, we used simple sequence repeat (SSR) marker fingerprinting and parentage 
analysis to determine the relative contribution of pollen parents in three 5-parent mixes (A, 
B and C) of £ grandis pollen used in four interspecific crosses with three E. urophylla seed 
parents. Sixteen SSR marker loci were tested, of which four were found to be informative 
enough to uniquely assign 247 of 248 tested F t polymix progeny to their respective seed 
and pollen parents. We found that the contribution of pollen parents to each polymix 
progeny set was highly biased in favour of one or two dominant pollen parents, while 
several pollen parents did not produce any hybrid progeny. The relative fitness of the pollen 
parents was remarkably constant in two crosses where the same £. grandis pollen mix was 
used to pollinate two different E. urophylla seed parents. These results suggest that severe 
pollen competition and/or differential hybrid viability occur in polymix crosses of E grandis 
and E. urophylla and that this phenomenon may be independent of the maternal genotype. 
It may also be necessary to re-evaluate polymix approaches to hybrid breeding in 
Eucalyptus and to investigate the possibility of using balanced pollen mixes of parents with 
equivalent fitness in interspecific crosses. 



Materials and Methods 



248 Leaf samples collected from 4 PMX families and parents 



CTAB DNA extraction 



(11 PARENTAL SCREEN 



(2) PROGENY SCREEN 



Amplification with all available SSR Amplification with informative SSR loci 
primers (Primer screen) chosen out of P rimer screen 



SSR analysis using Genescan 3.7 

and detection on an ABI 31 00 
sequencer (Applied Biosystems) 



Choose most informative SSR loci 
to use in progeny screen 



SSR analysis using Genescan 3.7 and 
ABI 3100 
(Applied Biosystems) 

+ 

Parentage analysis 



Statistical testing of pollen competitionVdifferential hybrid viability 
(x 2 test for deviation from 1:1:1:1:1 ratio) 



PMX Family Structure 
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Power of Detection 



Power to detect differentia) pollen fitness against poptJation 




Figure 1. Lower limits of detection of relative difference in pollen parent fitness (vertical 
green bar) determined by the range of family sizes tested in this study (horizontal green 
bar) 



Results and Discussion 

• Family sizes ranged from 37 to 82, which set a lower limit of detection of differential 
pollen parent fitness to a range of 1.2 to 0.7 (depending on family size) - (Figure 1) 

■ The four SSR loci selected were highly informative with individual allele frequencies 
mostly <0.15 and, in combination, allowed unique assignment of progeny to seed and 
pollen parents - (Figure 2) 

■ We observed a significant deviation from the expected 1:1:1:1:1 pollen parent ratio for 
all 4 families - (Figure 3) 

■ Some parents made no pollen contribution. This could be due to low pollen parent 
fitness (pollen competition or low viability of hybrid progeny produced by these pollen 
parents) 

■ We observed significant allelic transmission ratio distortion at two of the SSR loci (data 
not shown). This could be caused by selection on individual alleles due to linkage to 
pollen fertility or hybrid viability factors. 



Parentage Analysis Using SSR Marker Analysis 




Figure 2. Example of electropherograms obtained after SSR analysis of parents and 
polymix progeny on the ABI 3100 automated sequencer and Genescan software. The 
inheritance of parental alleles is indicated with vertical arrows. 
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Figure 3. Bar graphs indicating contribution of each pollen parent to polymix progeny as 
determined using SSR marker analysis. A significant tendency towards a dominant parent 
is visible. In family 2 & 4 the same PMX was used and similar results were observed, 
suggesting that the differential pollen fitness is independent of the maternal genotype. 
The P-values of the x 2 test at the 0.05 level of significance are indicated. 

Conclusions 

• Four informative SSR loci were sufficient to uniquely assign 247 of 248 polymix 
hybrid progeny to 3 seed parent and 15 pollen parents 

• Large differences were observed in the relative fitness of pollen parents in pollen 
mixes as measured by the relative contribution of each pollen parent to polymix 
progeny 

• This could be problematic for PMX breeding since some genotyplc combinations 
may be severely underrepresented 

• A possible solution would be to make "balanced" pollen mixes with parents of 
similar pollen fitness 
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